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On a Self-Exciting Process in Magneto-Hydrodynamics. 


By 
Hitoshi Taxrucui and Yasuo SHimazu. 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


In connection with the convection current model proposed by E.C. Bullard, the 
problem here considered is whether the magneto-hydrodynamical equations have a 
solution representing the self-exciting process which is appropriate to the model. In 
Section 8, the problem is reduced to solve the equation (3.9) with conditions (3.3) and 


(3.4). 


In Section 4, an approximate solution of the equation is obtained. 


In Section 5, 


this approximate solution is shown to be good enough for studying the present problem. 
The results obtained in this paper show that such a self-exciting dynamo is possible 
by which the earth’s main magnetic field is produced and maintained. 


1. The main part of the earth’s magnetic 
field and its secular variations are known to 
have nothing to do with things outside the 
earth. Here lies the most difficult point of 
the problem concerning the earth’s magnetic 
field. What is going on within the earth in 
order that its magnetic field can undergo 
sensible changes in a hundred years or so 
and in order that the field itself can, never- 
theless, be maintained? Recentiy, a simple 
and plausible explanation of the non-dipole 
geomagnetic field and its secular variations 
was proposed by W.M. Etsasser (1946, 1947) 
and E.C.Buiuarp (1948, 1949). The varia- 
tions are attributed to hydrodynamical eddies 
occurring in the earth’s core and to changes 
in electric current systems in the outer part 
of the core caused by the eddies. In the 
course of these studies, it became apparent 
that an improbably strong inducing main 
field in the earth’s core has to be assumed. 
This difficulty is removed, however, if one 
notices that, in the presence of primary dipole 
field and a differential rotation of conductive 
fluid, a secondary field much in excess of 
the primary dipole field will be produced by 
electro-magnetic induction. After examining 
all possibilities, Bullard came to conclude that 
the differential rotation is caused by the ther- 
mal convection of the fluid within the earth’s 


core. 


The success above referred to leads us to 
the idea that the whole magnetic field of the 
earth is ultimately produced and maintained 
by a mechanism of induction. ELsassmr dis- 
cusses the relation between the field and 
motion in a very general way and concludes 
that there may exist chains of relations in 
which a field H; interacts with a motion Vi 
to produce a field Hy which, in turn, interacts 
with another motion V2 to reproduce Hj. If 
it is possible to show the existence of sucha 
“Self-exciting dynamo’’ in the earth’s core, 
the main field might also be explained. All 
that is necessary is a weak initial field to 
start with, and the field will then grow auto- 
matically until the energy dissipated becomes 
equal to that supplied. Thus the problem 
is reduced to examine whether MaxweE..’s 
equations of electrodynamics possess solutions 
representing such a self-exciting dynamo. 

On the other hand, T.G. Cowxtne (1934) 
showed that the self-exciting process is im- 
possible if the motion of the fluid is rotation- 
ally symmetrical and is confined in meridian 
planes. In addition, Ensasser showed that 
the same is also true for a purely toroidal 
flow. Many believe that these theorems are 
foretastes of a stronger theorem that will 
prove the process is altogether impossible. 
This point is to be studied thoroughly. 

In this connection, BuLLARD’s convection 


74 


current model, which does not violate the 
above theorems, is to be reconsidered. In fact, 
Bu.LLarp has proposed a process by which 
the self-exciting dynamo may be possible. 
It is this convection current model. that we 
are going to study in the present paper. If 
it is possible to prove the existence of a self- 
exciting process by this model, we solve at 
the same time the following problems. 


(1) It is shown that Maxwe.u’s equations 
have a dynamo-solution. The questions raised 
by the studies by Cowxine and ELsAssER is 
put to an end. 

(2) A simple and natural explanation is 
given for the existence of the earth’s main 
magnetic field. Exsass—erR and BuLLARpD’s 
theory on the non-dipole field of the earth 
and its secular variations are afforded wider 
backgrounds. 


(3) New factors, such as the electric con- 
ductivity of the earth’s core or the magnitude 
of the convection current velocity there, are 
furnished for studying the internal constitu- 
tion of the earth. 


(4) The results obtained may also be ap- 
plied to the magnetic fields of the sun or 
of other celestial bodies in general. 


2. In order to investigate the maintenance 
of a magnetic field by electro-magnetic 
induction, the following equations must be 
solved under certain boundary conditions 
(AtFvin, 1950): 


curl H=47nJ=4nn(E+V x A), (alo) 
OH 

curl E= aie =(0). (2.2) 

diy H = 0, (2.3) 

div B = 4npc’e-?, (2.4) 


Customary notations are adopted. From (2.1), 
(2.2) and (2.3), we get 


y?H+4ne« curl(V x H)=0. (225) 


The convection current model proposed by 
BuLLarRD is shown schematically in Fig. 1. 
The expressions for the (7, 0, ¢) components 
of H, curlH and 2H are as follows. 
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For T,*° S-type magnetic field ing the magnetic fields coexist under certain 
0, boundary conditions disregarding their pro- 
| Baap. cesses of production. In Fig. 2, are shown 
oe, ee sin 2 the magnetic and velocity fields schematically. 
} — 2 a Pi 
; sin 00¢ : Ss 7 2c 
2 P2 
a( ) 
2 Peps fps SE JE a 
Pa r Ts 00 ? 
—6rT5 (28 24 Ps) , 
; : a(Spr2d P.*) 
curl H —( 2 AT2** .9¢,.9)\SIN 
iy iB ) 00 7” S Vi S Vig 
OG st Oh 2 Pal 
ow ie ) singdd ”’ C) © 
‘ Oe 
CT.2" S aT. S a(Si26 P) 
Fut —(# s—+6r— ) ee es g 
dr’ dr sin 60¢ p= 
‘ ar ) Fig. 2. Schematic representations. of the 
+( ” ) = a8 : , magnetic and velocity fields. : 
Pees sin 2A 2 eS. (2.8). 3. The system for which the solution is 


required is shown in Fig. 3. Since the equa- 
tions (2.1)- (2.5) do not change if we go from 
a system at rest to a uniformly rotating 


The expressions for the (7, 9, ¢) components 
of velocity fields are as follows. 


For T\°-type velocity field system of reference, we shall use the coordi- 
0, nates that are fixed to the rotating mantle. 
0 Therefore, V in (2.1) or in (2.5) denotes the 
y : oP. velocity relative to the rotating mantle. 
Vi(r)r aa ee Se ee ae (2 
For, S,?°-type velocity field 
—67rV2?"(7) cos 2¢ P# > 


» dV 27" -\ 0(cos 2¢ P”) 
V —(r ant +37V 2” \ eo 9 


Cn) 
—( 9 rend O66 oh ae10): 


In’ (2:6}— (2.10), Si2(7), -.3.-- , V2"(r) are un- 
determined functions of 7 only. The induc- 
tion process from S;° to T:° field by T,° mo- 
tion shown in Fig. 1 is that which was studied 
by Exsasser and BuLLARD in order to obtain Mantle 
a strong magnetic field in the earth’s core. 
Since this is an irreversible process, the self- 
exciting dynamo required cannot be construct- 
ed by the S,° and T° fields only. Fig. 1 
shows the chains of induction in their pro- : 
cesses of production. The condition for the The magnetic. field in the mantle which is 
- maintenance of the fields is satisfied by mak- an insulator ‘is 


Fig. 3,. The system of the core and 
the mantle. 
core: convective fluid with conductivity x. 
mantle: uniformly rotating insulator. 
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swn=A"(2),, 


T2(r) = T2?(1)= T2?5(7)=0, (3.1) 
where A’ is a constant to be determined. 
The boundary conditions are that H,, Hy, H», 
E, and E, are continuous and J,=0 at r=a. 
By the assumption of ‘‘viscous’’ liquid, V,, 
V, and Vs must vanish at r=a. These 
boundary conditions will be satisfied by put- 
ting 
oe a) et 


Vi(n)=V2(n)= (3.2) 
aT Wits (3.3) 
~ S.%r)+ - = =0 (3.4) 


at 7=a. 


In our present discussions, the velocity V 


in (2.5) is assumed to be given as a linear 
combination of (2.9) and (2.10) with (3.2). The 
induction mechanism proposed by BULLARD 
may be shown to be self-exciting, if we can 
make (2.3) and (2.5) satisfied by a properly 
chosen linear combination of (2.6) -—(2.8) under 
the boundary conditions (3.3) and (3.4). 


H in (2.5)=H in (2.6) with (3.4) 
+H in (2.7) with (3.3) 
+H in (2.8) with (3.3) 


=H,+H.+H;+Hi, say. (3.5) 


Since (2.6)-(2.8) are given so as to satisfy 


(2.3) independently, H in (3.5) will satisfy (2.3) 
also. Boundary conditions imposed on H are 
statisfied by Hi. Hz, Hs and H; in (3.5) inde- 
pendently, so we shall call them as coordinates 
(vector) functions. We shall now determine 
the radial functions S,°(7), in (2.6) - (2.8) 
by the following ‘‘orthogonal’’ conditions : 


ha 


eee ee 


2a 


\\r 2H+4ne curl(V x A)]- 
00 


et 


ae 


xr sin édrdod¢=0, (3.6) 


where H is to be taken as that in (3.5). This 
is a generalization of the so-called GALERKIN’S 
method (FrAzER, DuNcAN and Co.Luar, 1938). 
Executing integrations with respect to @ and 
od, we get 


ta, ae gate Pea 
‘s "(S8 S.0+ Te S0+Te" S°+Tes S8]-r°dr=0, 


1 AS oS fo oS ae jes ees 
{1S Te+Te Fo+Te T2+ Tes Te) Pdr=0, 


etc., (Sans 
where ° 

hy ae CAC an ae hohe ea f @Si° , 4 d§S,o dS\° 
S,° S°—— 1 P1 C ere 2 4 1 is be 1 

hip: zp ogatten a Set teary Ue pee oe ar “+25;°) , 

Pee 
T:° o°=0, 

Rie 5 axe bc eet EC SAS d daS,° 
2S S$.9 (dng)-l= pre Jr2) aS 26.77 2c, Some 
Ti295 S,° (Ane) CID nV Te 8 SEC De (AV Ta) (7 SE + 25,0), 
on ae 
SF hoe) =— 38 50 ere 

Pera 24 @T,° aT.» 

0 1) ses J 2 2 
T2 T2 5 (7 dr: ++ dr 2 ATS ’ 

om 1 288x6 /1 , 
Ti 8 T9 (dan) = — 2B Spe SVs 6 Tt 2 (AVE rT 8}rTe , 
, 96x 2,/0\, dlp dS 

So T3*5 (Ane) =— are {seve So < 4.28,9 ike Lge Vense [byTs9s 
ek rea 
T? T2°8 (dng) = —~ = Goes ("V8 T) + AV WT S\ rT 05 
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Sas = 2x6 LON ( , T3298 dT37°S 
EP pea BEB ON (4 PT 6 AT ares ig BEE. rms yap 


In (3.8), deo: in the expression for T2585 So : 
respectively. From (3.7) and (3.8), we get 


for example, means 0 and 1 for 7,2" and T275 


SS $A nerve Tato, 
| 1g SEE tne wv s) 
26 a [ap HVT + 2 (PVs) Te" |=0 , 
= SAP meee rhe: 98% 6 aeVior TS 
ee tte | op OVE TOES (CVT S]=0, 
8X6 [PT 4 ATE 6X6 4 yay 
— tne Fe {ar vate Z (SE + 250) + 21 2 v.88 ]} <0. 


These are the simultaneous differential equa- 
tions for the radial functions. 
Considering (3.2), we shall put tentatively 


Vi(r)=2 Vi. Eacelk 
V2?"(r) = Vee 1 (Z)] 


in (2.9) and (2.10), and take their sum as V 
in (2.5) and (3.6). In (3.10), Vi.o and V2.2, are 
undetermined constants. 


(3.10) 


4. Taking (3.3) and (3.4) into consideration, 


Thus our problem is reduced to solve (3.9) 
with (3.10) under the boundary conditions (3.3) 
and (3.4). It is easily seen that this is an 
eigen-value problem for Vi.o and V2... If 
we get real eigen-values for Vi.o and V2.2, 
by solving (3.9), we can thus prove that the 
process proposed by Buuuarp is self-exciting. 
It is, however, almost prohibitive to find the 
eigen-value by trial and error method. In 
the next section, we shall try to find them 
in an approximate method. 


we shall put 


sornmalt—3(2)} +a (EI $(Z)e a(S) LAME (Eft 
; T° (7) | Bo , Bi ( toy e 

a{1-(Z)] eae “jolie Ig) fen Ic") cos, (4.2) 

. F38(7)) Dy! \ Dy! D,,} 

in (2.6)-(2.8). Inserting (4.1) and (4.2) into (3.5), we get 

A, = Aoih.ot+ AiMi.it::++:: tA, Aint: t- >: ; 

HA, = Bofe.o+ Bifza+ +++: + BablanF****>° ; 
Ms ee ots (4.3) 


etc., 


Table I. Equations for the third approximation. 
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34 24 =e Te ie 48, 432 ' 
35-45” 35 35 5 35-55 35-55 36-65-13" 
2 28 pas = a4 26 72 48— 432 - 
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8-19 ae _ 26 248 8-102 ,| 9-482 | 72-38 _, . 
26-45-11” Ce 25-3 25-33 35-55-13" |(35)?-11-13" |(35)?-11-13 
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35-65 36-554 |. 35-65-13" 35 35 5 25-11 
. ret - 48 432 216 _ 176 _ 104 16 
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0 = a os al 
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where H.o, for example, denotes H which is 
obtained by inserting si(n=1-+(2) into 


4 

wee ia Sr ea =O,” Susy 2, 2. 22) 
in (4.3) satisfy the boundary. conditions im- 
posed on H independently, so we may take 
them as the coordinate vector functions in 
our approximate solutions. The results ob- 
tained by putting all the A;;’s with jy n 
equal to zero will be called the n-th appro- 
ximations. We shall now determine the con- 
stants Ao, Ai, .... by the following orthogonal 
conditions (loc. cit.). 


2 


\ | pH dna Vx A))- Hi; 
v0 0 


rt) 


x sin ddrdéd¢ =0, (4.4) 
where H is to be taken as that in (3.5) with 
(4.1) -(4.3). Inserting (3.10) and (4.1)-(4.3) 
into (4.4) and taking (3.7) and (3.8) into con- 
sideration, we get the equations of the follow- 
ing types: 

24 


8 2 
—6Ao— 3A = pee oR yDo 


WD 12 ae 
— 95 PPi— ey Vet eeee — 
etc. 


All of these equations for the third approxi- 


mation are shown in Table I. In Table I 
x=4nraVi.0, y=4nkaV2.2¢ (4.6) 
respectively. 


In order that these equations are compatible, 
the determinant formed by the coefficients 
must be equal to zero. Putting 


a=, (4.7) 


and omitting the coefficients for Ai, A», Bi, 
B:, Ci, Cs, D; and D,, we get the -first ap- 
proximation for y. Similarly we get the 
second and the third approximations. The 
results of the computations are shown in 
Table II. 


Table II. Eigen-values obtained by the 


approximate method. 


a eh ee ee 
a 1st 2nd 3rd 3rd 
0 169.36 269.7 Bye ra 0 
1 163.54 265.6 267 25.0 
10 82.322 94.80 | 111 105 
100 68 .664 72.05 92.8 870 
foe) 68.507 faloral 9253 co 


In order to express ‘‘magnetic REeyYNoLpD’s 
number’’ in terms of the maximum velocities 


of 71° and S.2*-types, Y= ey and X=Ga 
are shown in Table II. From Table II, we 
see that our method of approximation gives 
slowly increasing and converging real eigen- 
values. For the real eigen-value Y=323.2 in 
the case when a=0, all the A’s and D’s be- 


come zero, 1.e., we get no dipole field. 


5. In the last section, we have got real 
eigen-values for y and x, and it seems plau- 
sible that the process proposed by BuULLARD 
is self-exciting and is capable of maintaining 
the magnetic fields. As this is an important 
conclusion, we shall now make a numerical 
check for it. It is easily seen from Table I 
that as x00, T° 0. We shall now show 


y=92.3x > (our third approximation) is a 


very good approximation for the eigen-value 
for this case. 

Changing the independent and dependent 
variables by the following equations 


r=aé, 
be ie a eer (e) 
T2°5 (7) © { T28(E), 


Si (7) = Si’ (€) 


V;° (7) il Vi.9Mi° tS, 
= (5.1) 
Vestry 2 OU Vege Vo" (Ee); 


we get from (3.9) 
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9x 24 ¢ 
p EOP ge OE VTS 


T2 aT! 2 ie 2 e5 417.20’7T\ 2¢ 372° 26 
gO +6 Ge = 3% ge Set SEV ET H+ GEV IETE 
Gis dT,” 4. 2¢ 0 3 mee A 19 3% 5 
feo ge 7 aw V22°T2) + EVs ET» |-2x VET: 
2 dS ol 2c 0 
g OO 460 2a V8ET +3 Pi aeVa ae Ge +25!) + gel ageeve 35,1 | 
:. ae (5.2) 
Putting further 
To (E) = % T° €), (5.3) 
and making x— ©, we get, taking (3.10) and (5.1) into consideration, 
soe) = —3(e G6), 
P22 SSS 0 ? 
2 mes dT. j 
Ts (6)=79[£—£) Ge 47-11) TI. (5.4) 


OTe VIZe TY | 2850 Tr o0% ie, : viens 4 mo|— 
pete GE +e Ge tg reeled pe +c-upTe|=0-. 65) 


Boundary conditions (3.3) and (3.4) will be satisfied by putting 


T(E) =0, (5.6) 
and ae) = TEE + 10 Fe + 13 = 0 5.7) 


at &=1. (It is to be noted that our equation (5.5) with (5.6) and (5.7) is analogous to that 
in the problem of stability of viscous fluid motion (PEKERIs, 1948)). Expanding 7.°(&) in 
(5.5) in power series of &, we get the following two independent solutions 


Solution I for which 7.°=1, aT oie i ty 


dé AE Fare 
. aT! _ _ PTY __€TY 
Solution II for which Vo 1, T2= dE Fi aaa =0 (5.8) 


at €=0. The most general solution of (5.5) will be given by 


T.. (E) = A x Solution I + B x Solution II, (5.9) 
where A and B are undetermined constants. Integrating (5.5) with (5.8) numerically, we can 
determine the solution I and II, and using the numerical values thus obtained, we get two 


values of a from (5.6) and (5.7). The ratio 


a from “ 7) 


R= (5.10) 


Ql = 

a = 
a ene (5.6) 
¥ 


will become 1 for the exact eigen-value 
Similarly the determinant 


_| T:'(E=1) for Sol. I 7." (E=1) for Sol. I | | 
Ay Q (€=1) for Sol. I Q (€=1) for Sol. IL | (5.11) 


eS 


On a Self-Exciting Process in Magneto-Hydrodynamics. 9 


will become zero for the exact eigen-value Y. 
For Y=0, we get Solution I=1, Solution II 
=€, R=0 and 4=+18. This means that 
Y=0 is not the eigen-value for (5.5) with 
(5.6) and (5.7). In the same way, we get 
R=0.83, 4=+12 and R=1.88, 4=—63 for 
Y=92.3 (our third approximation) and Y=185, 
respectively. In view of these results, the 
exact eigen-value is considered to be at about 
Y=100. Our third approximation Y=92.3 is 
a good approximation to this exact value. 
Our approximate method for obtaining the 
eigen-value is, thus, shown to be fairly effec- 
tive. 


6. The eigen-values X and Y in Table II 
are of the orders of 107. These results are 
to be compared with those obtained by 
Exsasser and Buuuarp in their discussions on 
the secular variations of the earth’s magnetic 
field, i.e., c<=3-10-8e.m.u., V~10-2, a=3- 108, 
4nkaV~10?. Thus, using the results obtained 
in the present paper, we may construct a 
self-exciting dynamo by which the earth’s 
main magnetic fieid is produced and main- 
tained. The self-exciting dynamo is con- 
sidered to have the similar structure as sup- 
posed by Exsasser and Bu.uarp in their 
studies on the secular variations of the earth’s 
magnetic field. As to this self-exciting dy- 
namo, further studies will be made in forth- 
coming papers. 
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Density Distribution and Concentration of Heavy Materials 
in the Mantle of the Earth. 


By 
Yasuo SHIMAZU. 
Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


The fundamental assumptions underlying K.E. BULLEN’s study of the density distribu- 
tion in the mantle of the earth are the existence of a density discontinuity at the 
~ depth of 400km. and chemical homogeneities of the materials on both sides of this 
depth. The assumption of the continuous distribution of one homogeneous material 
throughout the mantle is known to lead us to an unreasonable result as to the physical 


ll 


state of the core. 


Recently, H. BROWN who studied the chemical compositions of the earth’s interior 
referring to those of meteorites, came to believe that the concentration of metal 
phase increases linearly with depth in the mantle. 

An earth model which has no density discontinuity in the mantle and which has a 
continuous concentration increase of metallic material towards the centre is possible. 


Such a model is investigated in this paper. 


The main results obtained are as follows: 
(1) The density at the bottom of the mantle may be between 6.289 and 7.970. 
(2) The total contents of metal phase in the mantle may be between 11.4 and 


22.9 percent. 


(3) The atomic ratio of iron to silicate in the whole earth will be 1.559. This 
value is not affected by the mode of distribution of iron phase between the mantle 


and the core. 


1. Fundamental Idea. 

We have knowledges about the velocities 
of seismic waves within the earth. That the 
density distribution within the earth can be 
derived therefrom on the assumption of hydro- 
static equilibrium was shown by K.E. BULLEN 
(1936), to whom we owe what is called BULLEN’s 
' density distribution. Whatever the mode of 
distribution may be, it must be consistent 
with the total mass of the earth on the one 
hand, and with the total moment of inertia 
about the axis of rotation on the other, the 
values of both of which have been known from 
astronomical observations. As we shall see 
later, BuLLEN’s method reduces itself to solve 
an integro-differential equation, and as a 
result, a unique solution can be obtained 
which fulfil the above two requirements, if 
a homogeneous earth is assumed. But, the 


propagation of seismic waves reveals the 
existence of a velocity discontinuity at 2900 
km. depth from the surface, associated with 
a density discontinuity. If there is such a 
discontinuity, the density distribution through- 
out the earth cannot be determined uniquely 
unless the density jump at 2900km. is 
known. But Buxtuen found that whatever 
reasonable jump in density at 2900 km. is 
assumed, no continuous density distribution 
in the mantle from the surface down to 
2900km. depth can yield the conceivable 
result as to the state of the core. Thus, he 
came to assume the existence of a density 
discontinuity at 500km. depth. Later, he 
(ButLEN, 1940) was led to consider that the 
density changes rapidly within a transient 
layer between 400-900km. depth in the 
mantle, and that the hydrostatic equilibrium 
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is not satisfied within this layer. Besides, he 
assumed chemical homogeneities on both 
sides of this transient layer. Thus, finally, 
he could obtain a consistent density distribu- 
tion, adjusting the density variation within 
the 400-900km. transient layer. But this 
supposed density variation in the transient 
layer cannot be said very certain, for it is 
the so-called second order change, only the 
gradient in density being varied. Although 
the 400 km. discontinuity is often claimed by 
several seismologists, BuLLEN’s distribution is 
after all nothing but one of the important 
possible solutions. 

Recently, H. Brown and C. PLATTERSON 
(1947, 1948) studied the chemical compositions 
within the earth as inferred from those of 
meteorites. In contrast to similar discussions 
already made by many geochemists, which 
are purely chemical and qualitative, Brown’s 
is based on physical chemistry in a much 
more quantitative way. This has been made 
possible, mainly owing to the recent develop- 
ments in the technique of quantitative micro- 
analysis. It has been shown by Brown that 
if the observed distribution of elements in 
meteorites is regarded to correspond to that 
in equilibrium state, then the temperature 
and pressure necessary to give rise to that 
distribution of elements can be estimated. 
The data so far indicate strongly that all 
meteorites have their origin in a planet similar 
to the earth in general physicochemical charac- 
teristics. In what proportion Fe and Ni are 
distributed between silicate and metal phases 
are controlled by pressure only. Knowing the 
relation between the equilibrium constants 
and pressure, the distributions of Fe and Ni 
within the earth can be estimated. In such 
a logical manner, Brown obtained the radial 
distribution of metal content in the mantle 
of the earth. These results of Brown indicate 
an almost linear increase with depth of the 
metal content by weight which is mainly due 
to Fe. Such an increase of heavy metal 
towards the centre of the earth must also 
have influence upon the density distribution 
within it. Instead of the BULLEN’s model with 
a discontinuous mantle, a model having such 


a concentiation increase may be also possible. 
The object of the present paper is to find 
such a solution. 

2. Fundamental Equations. 

If V,; and V, are the velocities of longi- 
tudinal and transversal seismic waves respec- 
tively, they are expressed by 


V; + (K+5n) vi=# QD) 


2S 

= 3 
where p is the density, » the rigidity and K 
the bulk modulus defined by 


1 eles 


—=— 2 Zea: 
From (2.1), we obtain 
K 4 
wire Vio @ (2.3). 


If hydrostatic equilibrium is assumed within 
the earth, then 

where 7 is the distance of any point from the 
earth’s centre, p the pressure, g the attraction 
force per unit mass at this point. Also 


“3 a (k?: constant of gravitation) (2.5) 


(2.4), 


if the earth be taken as spherically stratified. 
M(r) is the mass enclosed within a sphere of 
radius 7 which is concentric with the earth 
and 

do __00 Op , Oe A6 

dr Op gr 006 Or 
In the above expression, the first term re- 
presents the density increase caused by hydro- 
static compression. 6 is an independent vari- 
able representing the concentration of heavy 
materials. Physically, it will be convenient 
to consider 6 as the contents ratio of heavy 
material by volume. 

Inserting (2.1)}-(2.6) into (2.7), we obtain 


d Rb? r 
=— Lae | may—4n\" ortdr |+ 98 06 


+ (2.7). 


The density distribution with respect to depth 
can be obtained by integrating (2.8) numeri- 
cally. (o/K) is a known function of 7 by 
(2.3), if V; and V, are known. If we omit 
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the second term of the right hand side of 
(2.8), the BuLien’s fundamental equation is 
obtained. 

Following BuLLEN, we start the integration 


from 35km. depth towards the centre. The 
conditions given are 

radius of the earth: a=6338 km. 

total mass: M(a)=58.97 x 1076 gr. (2.9). 


total moment of inertia: 
i (e)—(9nl4s 10 rersem 2: 


Let the initial value of density at the start- 
ing point be 

0(a)=3.32 (2.10). 
The layer above 35km. depth is removed 
because the condition of hydrostatic equili- 
brium may not be satisfied in this layer. 
When Buen put the density p(a@) as 3.32, 
there were some criticisms on this particular 
value, for instance, by F. Birow (1939). But 
BuLLEN (1940) pointed out that there were 
many geophysical arguments in favour of 
this value, and showed that even if the initial 
density was in error by 0.1, this would affect 
the densities below this depth by no more 
than 0.5 percent. 

Now, we shall define a quantity 

2= 4s (2.11). 
In case of a uniform sphere, z=0.40. From 
(2.9), we see that z=0.3342 for the whole 
“earth and this value indicates a strong increase 
of density towards the centre. Also 


z<0.40 (2.12) 


must be satisfied in the core, because the 
increase of density towards the centre will 
also be expected in this part. This additional 
condition as to the state of the core is useful 
in determining the density distribution in the 
mantle part. But the condition z< 0.40 gives 
only one of the two extremes of the density 
distribution. The condition for the other 
extreme will be discussed in section 6. 

3. Numerical Integration. 

We shall transform the independent variable 
from 7 to 


ae me. 
a=1 7 (3) 


The mantle corresponds to 

x=0~0.452 (GAY 
Moreover, we shall put (00/06) (0d/0x) as a 
linear function of 2, such as C)+C, x, where 
CG and -C, are unknown constants. The 
density can be obtained from (2.8) as an 
expression which includes the two unknown 
parameters Cy) and C,. These parameters 
must be so determined that the condition 
(2.12) is satisfied. 

To get the value of (e/K) by (2.3), we shall 
use the distribution of V, and V, published 
Dyas JEFFREYS (1939). For integrating (2.8) 
numerically, the method of expansion by 
power series will be adopted. The mantle 
part is divided into six shells in order to. limit 
the errors in the integration within one percent. 
The solutions for each shell are connected. 
In a shell in which the variation of (p/K) is 
large, one step of integration is taken small. 
The six shells are shown in Table I. Omitting 
the details of numerical integrations, we shall 
show only the densities at various depths in 
Table I. 


Table I. Densities at Various Depths. 
xv depth (km.) 0 

0.00 | 33 8.82 

0.06 | 413 3.6684 6.3100C)+-0.187160C, 
0.09 | 603 3.8054 9.5962C)+.0.236420C, 
0.12 | 794 3.920 +-12.9208C') 4-0. 288659C, 
0.28 1808 4.460 +31. 5900C) + 0.642648C, 
0.44 2882 = 4, 968 +. 49. 7999C') +. 1.009977C, 
0.452} 2900 5.009 +-51.0976C)+1.019388C, 


Integrating the values given in Table I again, 
we get ‘ 


mass of the mantle 
= (36.179 + 185.526C)+3.84261C;) 
lO orca (eZ) 
moment of intertia of the mantle 
=(63.729 +265.870Co+5.40926C1) 
x 10" ot. Cm (3.0). 


Subtracting (3.2) and (3.3) form (2.9), we get 


mass of the core 
M,=(22.791—185.526C)—3.84612C;) 
; 1078 sr S:4); 
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moment of inertia of the core 
I, =(15.414—265.870C)—5.40926C1) 
x 104 gr. cm?. (3.5). 
If the ere is chemically homogeneous as 
BuLuEN first considered it is, then 


G;—G—0; 
From this, it follows 
2= (eon 


for the core, indicating that the homogeneous 
mantle model is unreasonable. Next we shall 
assume that 6, the content ratio by volume 
of heavy material, is a linear function of 2. 
Since 6=0 at the surface, we may put 

(3.6). 
If the increase of bulk density by concentra- 
tion of heavy materials does not depend on 
depth, or the difference in bulk moduli of 
the mother and metallic phases is ignored, then 


d=kxz («: constant) 


Oa Oy Const. PAG 


Bound to the condition-(2.12), we must have 
Co = 0.025051 (3.7). 


Since the total mass of metallic material in 
the mantle M,, and the total mass of the 
mantle M are 


My, = 185.526 x 107°C) =4.648 x 106 (aaa) 
M = (36.179+4.648) x 10?°= 40.827 x 1076 (3.9), 


respectively, the total content of metallic 
material by weight in the mantle must be 


M, 


ap = 0.1138 (3.10). 


4. Concentration of Fe. 

If (00/06) is independent of depth, the critical 
values for the content ratio can be determined 
without signifying what is the concentrating 
material. But if (00/00) is a function of depth, 
we must specify the material. Here is the 
reason that the bulk-moduli of materials must 
be taken into consideration. Now, we assume 
that the mantle is composed of silicate and 
Fe phases as follows, 


volume 
acu bo leahaee ratio 
Silicate Pine) 3.82 | 1.085 x 10” 1-9 
Fe phase 7.87 | 1.683 x 102 6 


cere 


As Brown has shown, nearly 90 percent by 
weight is due to the Fe metal phase in the 
mantle. Thus our assumption that the metal 
phase is exclusively composed of Fe may not 
be in a serious error. Then the density of 
the mantle material is given by 


0=0re0+0si(1—94), (4.1) 
OF ee 
or a0 = 0re— Psi > 


where pr. and ps; depend on the depth through 
pressure. On the other hand, in the theory 
of finite elastic strain (Birca, 1938; SHImMaAzu 
1949, 1950), we have 
p= poze). (P=Kee 

where e¢ is the strain caused by external 
hydrostatic pressure P and 0), Ko are density 
and bulk modulus at the state of no external 
pressure respectively. Eliminating e from 
the above equations and inserting the result 
into (4.1), we obtained (09/06) which is shown 
in Fig. 1. Since the pressure distribution 
within the earth is not sensitive to the 
density distribution, the transformation of the 


and TAKEUCHI, 


ap/as 
6 
S 
4 
=) 
cS 1.0 15x10 12 
: do 
Fig. 1. 59 38 function of P. 
(Ore —Psi~P) relation into the (o7.— 0si~depth) 


relation can easily be made. Because (Ore— Osi) 
is approximately linear with respect to P 
which is also approximately linear with respect 
to depth, we obtain 


00 
Og 4-986 +-3.654 « 


(4.3). 
From this and (3.6), we get ; 


Co+C «=(4.586-+3.654 2) «, 
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Table II. - Densities at Various Depths. 
er a ee A a ae Sls Aine | i oe a 
x Depth (km i 2 : Metal phase 
pense “ oe, (3) &) content by weight 

0.00 33 3.320 3.320 8.320 3.320 0.000 

0.06 413 3.668 3.827 4.026 3.643 0.041 

0.09 603 3.805 4.046 4.368 4.186 0.059 

(ge 794 3.920 4,244 4.670 4.480 0.076 

0.28 1808 | 4.460 5.251 6.291 5.187 0.151 

0.44 2882 | 4.968 6.251 7.857 5.645 0.201 

0.452 2900 I 5.009 | 6.289 ie OkO) 5.676 0.204 
a 


Co=4,586 «, C,=3.654 « (4.4). 


We know that the effect of depth upon 
(00/00) is rather small. This is the reason 
that the difference in bulk-moduli of the two 
phases is important. From (3.4), (3.5), (4.4) 
we obtain 


& > 0.005376 (4.5). 


The density distribution thus obtained is 
shown in Table II and Fig. 2, in which 

(1) continuous and homogeneous model, 

(2) critical values for the model with con- 

centration variation, (z=0.40), 

(3) same as (2) for z=0, (See Section 6) 
(4) BULLEN’s distribution (BULLEN, 1947). 
The metal phase content by weight given in 
Table II is shown in Fig. 3 by a thick curve. 
The total metal phase content is 11.4 percent 
in the whole mantle. The metal phase 
content estimated by Brown is shown also in 

Fig. 3 for comparison. 

Now we shall study the distribution of the 
metal phase in detail. That (K/p) has rather 
_-abnormal values between 400-900 km. depth 
is an observational fact. Since we have 
assumed that the volume content of Fe phase 
is a linear function of depth in this paper, 
the abnormal value of (K/o) is attributed to 
that in the elastic constant K. According to 
Brown, while the metallic Fe phase content 
increases almost linearly with depth, Fe in 
silicate phase as chemical compounds (com- 
bined iron in Fig. 3) shows a remarkable 
abnormality at 400-500km. depth. This is 
- also true for Mg content which is another 
important constituent in the silicate phase. 


T a 
4° esl 


al +2 3 


Fig. 2. Density distribution 
within the earth. 


3x 10°km 


Metal phase content by weight. 


Fig. 3. 
@ Metal phase by Brown. 
© Combined iron in silicate phase 
by BROWN. 
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Although the above circumstance has been 
ignored in our simplified model, it is certain 
that the density distribution is affected by 
abnormal behavior in chemical compositions 
at 400-500 km. depth. Previously, J.D. BERNAL 
(1936) tried to explain the discontinuity at 400 
km. depth by polymorphic transition of con- 
stituent materials. It is remarkable that there 
is such a chemical irregularity at 400-500 km. 
depth besides the physical one which is known 
from seismology. It may also have some 
connection with the discontinuity in electric 
conductivity found at 400km. depth (RIKITAKE, 
1950, 1951). But on account of the difficulty 
of discussions about the relation between 
concentration of metallic phase and elastic 
constants, we shall not try to discuss the 
problem in detail in this paper any more. 

5. Atomic Ratio. 

In order to discuss the distribution of ele- 
ments, it is necessary to express the results 
obtained in section 4 in terms of the atomic 
ratio, which is defined by M/A, where M is 
the content ratio by weight and A the atomic 
weight (molar weight in the compounds). 
Regarding the core as composed wholly of 
metal phase, we get from (3.2), (3.4): 

mass of Si-phase in the mantle 

= 36.18 x 1078 gr., 


mass of Fe-phase in the mantle 
=(185.526C)+3.84261C;) x 1078 =4.65- 107% gr. 
mass of Fe-phase in the core 
=M, =18.14 x 10 gr. 
Then 


My, = 36.18 x 10 gr., 
My, =4.65 x 10° (in the mantle) 
+18.14x 10° (in the core) 
27.19 <10%6 pr. 
in the whole earth. On the other hand, A= 
142 (mean value for rocks), A=56 (Fe). Then 
the atomic ratio in the two phases becomes 


Nre _ Mre /Msi _ 
Ne = A | A = 1-599 


aS. at (5:1). 
Considering the increase in A with depth by 
the gravitational concentration (Ramsry, 1948), 
Nr. / Ns will be slightly larger than the value 


given above. Goxpscumipr obtained N,,/Ng; 


0.89 (ELsassER, 1950) for meteorities, Unsold 
2.7 (ELsAssER, 1950) for the solar atmosphere 
by spectroscopy and Brown 1.7 for the earth 
by physicochemical considerations (1949). As 
My, is determined by the first term of (2.7) in 
our discussion, My, is obtained uniquely for 
the whole earth. Thus the result (5.1) is un- 
changed even if the distribution of My, be- 
tween the mantle and core is varied. 

6. The Density Distribution for Another 

Extreme Model. 

Another extreme case as to the state of the 
core will be given by the condition z=0. 
This amounts to mean 7J,=0. Corresponding 
to this case, we get from (4.4) 

«<0.01244 (6.1), 
and 
Ox =0-452 = 7.970, 
M, > 121082% 10" or: 


Mre __ 
Mre+ Msi rake 
These results are given in the column (3) of 


Table II. 


for the mantle. 


Acknowledgements 


The writer wishes to thank Prof. C. Tsuv- 
Bol of the Geophysical Institute, Tokyo Uni- 
versity, for his continued encouragements. 
The present study has been suggested by Dr. 
H. Taxrucui of the same Institute, to whom 
the writer’s thanks are due. He also wishes 
to express his thanks to Prof. E. Minami of 
the Chemical Institute, Tokyo University, for 
his kind advices. 


References 
BIRCH, F.: 

1938 ‘The effect of pressure upon the elastic 
parameters of isotropic solids, according 
to Murnaghan’s theory of finite strain’’. 
Jour. App. Phys., 9, 279. 

“The variation of seismic velocities within 
a simplified earth model, in accordance 
with theory of finite strain’. Bull Seism. 
Soc. Am., 29, 463. 
BROWN, H. and PLATTERSON, C.: 
1948 ‘‘The composition of meteoritic matter. 
III. Phase equilibria, genetic relations 


1939 


Density Distribution and Concentration of Heavy Materials in the Mantle of the Earth. 17 


and planet structure’’. 
85. 
“Tables of relative abundances of nuclear 
species’”’. Rev. Mod. Phys., 21, 625. 
ISON, IMSS 

1936 ‘The variation of density and the elliptici- 
ties of strata of equal density within the 
earth’’. M.N.R.A.S., Geophys. Suppl., 3, 
395. 
“The problem of the earth’s density 
variation’. Bull. Seism. Soe. Am., 30, 
235. ; 
“Tntroduction to the theory of seis- 


Jour. Geol., 56, 


1949 


1940 


1947 


mology’’. Cambridge Uniy. Press, 218. 
ELSASSER, W.M.: 
1950 <‘‘Harth’s interior and geomagnetism’’. 


Rev. Mod. Phys., 22, 1. 
JEFFREYS, H.: 
1939 “The times of P, S and SKS, and the 


velocities of P and 98”. 
Geophys. Suppl., 4, 498. 
RAMSEY, W.H. : 


M.N.R.A.S., 


1948 “On the constitution of the terrestrial 

planets”. M.N.R.A.S., 107, 406. 
RIKITAKE, T.: 

1950 <‘‘Hlectromagnetie induction within the 
earth and its relation to the electrical 
state of the earth’s interior’. Bull. 
Earthq. Res. Inst., 28, 45. 

1951 ‘Electromagnetic induction within the 


earth and its relation to the electrical 
state of the earth’s interior’. Bull. 
Earthq. Res. Inst., 29, 219, 268. 
SHIMAZU, Y.and TAKEUCHI, H.: 
1949 “Effect of tension upon the elastic para- 
meters of isotropic solids’. Geophys. 
Notes, Tokyo Univ., 2, No. 23. 


ys hig f ij 
ve, ane eo ws =. 7 : 
. } ‘ ~ Ste Ne 
£ ot Oe » 9 be See 
oA pa 
nt “ee Dinan cinnget- dedi 
; Mie Fel ys (htter 2 Soy AE 


‘7 7 ‘, 
at at Ne ‘ Ms 1 yr Lia mahi =) 
: 7 


ae Tri Aces at jb, ui ak Z i f 
, PA a 
an wre er a) om htt ed mee a 
Pietaten® Qreeisbiewes. oye ae cae 
¥ ae Sigs r Ve ot ria ee a. Nfs oie “hy ~ - 
a ~@ a ? aw 
x : Le e: os 
cae ab We aie eel Lee RE - 
a ~ le a : a. > 
— ots $F: «aa! gh Leds Tg at Se oo he wag ——— a 
aap th eS? gros teers Jy oh i> tue =~ 
BR MSW dlnctit 88 wae te OL 
. © a ‘ 
Cee RO SEE ak eet See 


ie . a Sp spew oe | 


JOURNAL OF PHYSICS OF THE EARTH, VOL. ENG e Lobe. 


Griuneisen’s Parameter and the Depth of Isostatic Compensation. 


By 
Haruo Mixt. 
Geophysical Institute, Faculty of Science, Kyoto University, Kyoto. 


Abstract 


From the considerations on GRUNEISEN’s parameter, it has been found that the 
condition of hydrostatic equilibrium is satisfied beneath the depth of 90km. from the 
earth’s surface, although it is not satisfied above that level. This depth can be inter- 
preted as representing that of the isostatic compensation. The value of 90km. was 
deduced by using the velocities of seismic waves in Japan, so that it can claim no 
general yalidity. Nevertheless, this line of approach may afford a new method for 
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estimating the depth of compensation from data other than gravimetric. 


The GRUNEISEN’s parameter of a material is 
defined as 


Pe —ovUlog ¥,,)/ lor VW) i525. ose .2 8: (1), 
where 

Bryne ee {(INz/4nV)/ (Up? +2u5-3) 8 . .(2), 

Nrz....LOSCHMIDT number, 

Vanes. volume per one mol, 


.velocities of longitudinal and 
transversal seismic waves. 


We assume that the fundamental idea of 
GRUNEISEN can be applied to the materials of 
which the earth’s crust and mantle are com- 
posed. We divide the earth’s crust and mantle 
into a number of strata and assume that 
neither the crystal structure nor composition of 
the material changes within any one of these 
strata and that the GrineIsen’s parameter 

is therefore constant in each stratum. Then 
from (1) and (2), we get 

(Wp +20.-)i[(Vp?-+20.- = (Vil Varo"? (3), 
where the suffixes 1 and 2 correspond to the 
ypper and the bottcm surfaces of the stratum 


under consideration respectively. 
The condition of hydrostatic equilibrium is 


written by 


Vp Ueecns 


OD ase Bepede, oP aT 
9°= "ar odo dr ' OT.ar** 


2. .(4), 


in which 


(5), 


(kK: Volume modulus). . 


can be evaluated from seismic data and 


from GRUNEISEN’sS equation of state. 
Combining the above three relations, (4), 
(5) and (6), we get the following equation ; 


40+) B AO) 


2 Gi TG 1 1 elie 


oO re (GA 
TT Che AA fen ee 
where 
[Es Ace specific heat at constant volume 
per one mol, 
Ae, atomic weight, 
and 
Mowe a parameter depending upon crystal 
structure and defined by 
AA 
re a) 


If the a four conditions, viz., 


(1) ei) that is, 
creases with increasing depth, 

(2) the crystal structure and composition 
of material do not change in any 


stratum, 
(3) GriNEISEN’s parameter has no extra- 


ordinary value, 


the temperature in- 
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(4) the condition of hydrostatic equilibrium 
is satisfied, 

are satisfied in the interval of integration 

between 1 and 2, then it follows from (7) that 

the following inequality must hold in that 


interval : 
50+) [>| dal To 


The both v) of the Ae: inequality can 
be evaluated from observed velocities of seis- 
mic waves propagated through the earth’s 
interior. The results of calculations based on 


actual seismological data in Japan (Honpa, 
1942) are shown in Fig. 1. 


CGS 
4x10" 
1 V2 
————— 1-= 
3x10 Pp Te-h Ve 
h 
/ 
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 G=1.0 
Noo 
10), eee, errr ro) Snreoerel Conan cuen 
/ 
beh, 
t=20 
0 00 200 300 400 


Depth in km 
Ui bens) Abe 
The points in the figure correspond to the 
middle points between the upper and lower 


limits of integration indicated as suffixes 1 and 
2. We see in the figure that 


- oF (1-7 ) 


decreases at first steeply: with increasing 
depth down to 90km., beyond which the 
change is rather gentle. We see also that the 
inequality (8) is not satisfied before this depth 
of 90km. is reached. This shows that one 
or more of the conditions prescribed above 


is not satisfied by the material in this part of 
the earth. 


We may put the condition (1) out of ques- 
tion, since its validity is already an esta- 
blished fact. The existence of polymorphic 
transition under the influence of pressure and 
temperature has been ascertained experi- 
mentally. But considering the results of high 
pressure experiments, it is difficult to conceive 
that this transition occurs from the earth’s 
surface to the depth of 90km. (about 30,000 
atm.) continuously throughout a considerable 
extent. The same can be said with respect 
to the composition of the material. There- 
fore the condition (2) is also taken as valid. 
GRUNEISEN’S parameter is known to have the 
value of 1~2 under ordinary conditions, and 
is almost constant against temperature varia- 
tion. Unfortunately no experimental result 
has been available for the effect of pressure 
on it. That it has an extraordinary value 
throughout from the surface to the depth of 
90 km. is, however, difficult to conceive, con- 
sidering that the equation (8) is satisfied with 
its ordinary value 1~2 beneath the depth of 
90km., where the pressure is higher than 
that in the portion where the equation (8) is 
not satisfied. Therefore the condition (3) 
must be regarded to be satisfied. The condi- 
tions (1), (2), and (3) being taken as satisfied, 
we have come to conclude that it is the 
condition of hydrostatic equilibrium that is 
not satisfied in the part shallower than the 
depth of 90 km. 

There are two different opinions as to the 
model of isostasy, the one is Prarr’s and the 
other Atry’s. Many efforts have been made 
in vain in order to see which of these opinions 
is better for explaining plumb-line deviations 
and gravity anomalies. Our results stated 
above give a new method for determining the 
depth of compensation from data other than 
gravimetric. 

The writer wishes to express his hearty 
thanks to Prof. K. Sassa and Prof. E. Nisat- 
muRA‘ for their instructions, encouragements 
and many valuable discussions. 


Reference 
H. Honpa. : 
1942 “Seismic Waves” (in Japanese), Tokyo. 
pp. 80. 
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Three Types of the Distributions of Temperature and 
Density in the Interior of the Earth. 


By 
Haruo Mik1. 
Geophysical Institute, Faculty of Science, Kyoto University, Kyoto. 


Abstract 
Based on the considerations on GRUNEISEN’s parameter as in the preceding paper 
it bes been shown that three types of the distributions of density and temperature in 
the interior of the earth can be considered according to the following three suppositions 
(a) GRUNEISEN’s parameter has an extraordinary value throughout the earth’s 


mantle. 
(b) In the stratum between the depths of 400 km. and 800 km., GRUNEISEN’s para- 


meter has an extraordinary value. 

(c) GRUNEISEN’s parameter has an ordinary value in the part shallower than the 
depth of 400 km. and an extraordinary value in the deeper parts. 

In (a), the temperature decreases with increasing depth between the depths of 400 
km. and 800 km. In (e), the physical states or chemical compositions of material change 
at the depth of 400 km. dividing the earth’s mantle into two parts. 


In the preceding paper (Miki, 1952), the nu- CGS 
merical values of the both sides of the equation 4 io” 


V2 | K i V2 | 
saa | armas (1-y;) 

were calculated down to the depth of 400km., 
using seismological data obtained in Japan, 
and the results were shown graphically in a 3,403 
figure. In the present paper, the calculations 
will be extended to deeper parts of the earth’s 
mantle. Necessary materials for this purpose 
were taken from BuLLEN’s ‘‘Introduction to 
the Theory of Seismology”’ (1947). Poisson’s 
ratio was also calculated from the same 
material. The results of calculations are 
graphically shown in Figs. 1 and 2. We see 
in Fig. 1 that between the depths of 400 km. 
and 800km., the conditions stated in the 
preceding paper are not satisfied. 103 

Since the condition of hydrostatic equili- 
brium must undoubtedly be satisfied in these 
depths as usually accepted, we are led to 
conclude that one or more of the following 
three conditions must be true between the 
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Fig. 2. Porisson’s ratio at various depths. 


(2) The crystal structure, or the composi- 

tion of material, or both, vary. 

(3) GRUNEISEN’s parameter has an extra- 

ordinary value. 

These anomalies of physical states corres- 
pond to the gradual increase of Potsson’s 
ratio from the depth of 400 km. to 600 km. 
and steep decrease from the depth of 600 km. 
to 800km. as shown in Fig. 2, although at 
present, we know no reason by which to 
explain this correspondence. It is obvious, 
however, that this anomaly which begins at 
the depth of 400km. and extends to the 
deeper parts of the earth is closely related 
to the 20° discontinuity which appear in the 
gradient of the travel-time curve of seismic 
waves. 

Since we have no decisive data to examine 
the truth of the three conditions stated above, 
let us consider the following three cases as to 
the value of GRUNEISEN’s parameter in the 
depth of the earth. 

(a) GRUNEISEN’s parameter has an ordinary 

_value throughout the earth’s mantle. 

(b) It has an extraordinary value between 

the depths of 400 km. and 800 km. 

(c) It has an ordinary value in the part 

shallower than 400 km. and an extra- 
ordinary value in the deeper part. 
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It is likely that the successive polymorphic 
transitions and the variation in composition 
of the material do exist in the interior of the 
earth (NISHITAKE, 1951). But for the sake of 
simplicity, let us assume that these transitions 
and variation occur only at the depths of 
400 km. and 800km. in case of (b), and at 
the depth of 400km. alone in case of (c). 
Under these assumptions, we get the distribu- 
tions of density and temperature from the 
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Fig. 3. Density at various depths. 
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equations (3) and (4) in the preceding paper, 
corresponding to the above three cases. The 
results of calculations are graphically shown 
in Figs. 3 and 4. In each figure, we take 
y¢=5.0 as an extraordinary value of Grijnet- 
SEN’S parameter and y¢=2.0 as an ordinary 
one. The amount of increase in density at 
the depth of 400km. and the amount of 
decrease at 800 km. must be adjusted so that 
the mass and the moment of inertia of the 
earth calculated from the density distribution 
obtained by the above procedure agree with 
the observed ones. 

From the two figures, we notice the follow- 
ing facts. 

In case of (a), between the depths of 400 
km. and 800 km., the temperature decreases 
with increasing depth, and the total mass of 
the earth’s mantle comes out too small. 

In case of (b), there is a stratum having 
an extraordinary value of GRUNEISEN’S para- 
meter. 

The case of (c) implies that the physical 
states or the chemical compositions of material 


change at the depth of 400 km., dividing the 
earth’s mantle into two parts. 

These results seem to be very suggestive 
although successive polymorphic transition 
and the change in composition were neglected 
in the computations. Unfortunately we cannot 
proceed further, because of lack of available 
materials. 

The writer wishes to express his hearty 
thanks to Prof. K. Sassa and Prof. E. Nisut- 
muRA for their instructions, encouragements 
and many valuable discussions. 
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Rheological Properties of the Earth’s Crust and Alluvial Layers 
in Relation to Propagation of Seismic Waves. 


By 
Akira KusorTera. 
Geophysical Institute, Kyoto University, Kyoto. 


Abstract 


The earth’s crust is generally divided into three layers, the upper two of which 
are called granitic and intermediate layers, while the deepest is believed to be composed 
of peridotite. The visco-elastic properties of the three layers have been studied in the 
present paper from the velocities and amplitude changes of seismic waves which are 
propagated through them, Principal results obtained are as follows. 

1. The attenuation coefficient and relaxation time in each of the three layers were 
computed from the amplitude changes of initial P-waves in case of near earthquakes 
with the results as follows: 


Attenuation coefficient Relaxation time 


Granitic layer Ka=7.9 10-7 em=1 0.45 sec. 
Intermediate layer K*=1.3x10-* 3.5 
Periodotite layer Ky,=2.6x10-8 16.5 


2. The velocity of P in the granitic layer is about 5.4-5.75 km./sec. from natural 
earthquakes (Py seis.), while it is 6.0 km./see. from explosion experiments (Py, expl.). 
This difference can be accounted for by the dispersion of seismic waves in the visco- 
elastic medium, since the periods of the waves in the two cases are decidedly different. 
On the other hand, searecely any difference can be found between P-wave velocities 
from natural earthquakes and artificial explosions in the intermediate and peridotite 
layers. 

8. The viscosities of the layers were calculated by the methods recently developed 
and were compared with the values obtained by other authors. 

4. Similar studies were made of the results of artificial explosion experiments on 
alluvial layers, such as have recently been attempted at various places in Japan by 
members of our institute. It was found that earthquake damages on the surface of 
the alluvial layers are closely connected with the rheological properties of the layers. 


1. Visco-Elastic Body S=GF+L7 oO pts SEN an 
a) Models of the Visco-Elastic Body. Ot 
For the’ purpose of investigating the rheo- 2) Maxwet.’s models 

logical characteristics of the earth’s crust, it 3S G 

will be appropriate to start with some synoptic - ap at F —F Sh eae dea) 


reviews on various models of the visco-elastic 


body in general. Typical classical models where S: ples 

are those proposed by VoierT, Ketvin and F: strain 
Maxwe.u. The stress-strain relations accord- G: elastic constant 
ing to these models are; 7: viscosity. 


1) Vo1er’s (KELVIN’s) model 7 is called the internal friction in the model 
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of Vorer and the firmoviscosity in that of 
MAXWELL. 

GutTensperG and ScaLecutwecG (1930) and 
Sezawa and Kanai (1940) calculated the 
viscosity of the earth’s crust using VoicT’s 
model. On the other hand, Scoweypar (1912), 
Ventna Mernesz (1934), Hasken (1936) and 
others obtained the viscosity of the earth’s 
mantle using Maxwe.w’s model. 

In place of these classical models which 
have but low approximations, those which 
have higher approximations have been pro- 
posed recently by Frenkett and OBRAsTZov, 
Desye, Osipa (1941), Topa (1947), and others. 
Instead of one definite value of the relaxation 
time in the classical models, several different 
values can be conceived according to these 
new models. 

In the present paper, Topa’s model will be 
adopted because it is easy to deal with. The 
relation between stress and strain tensors 
according to this model is 


Sip=2uPiz+a 01; 4 fal anegel a) <176)'—) (ES) 
where 4 and y» are following operators : 
0 aro 
ries nape 2 ae 
3 0 
Kg+tt Caw It+tay 
: 0 Hate (ls) 
ree ot 
— a 
1 @=)) 
and oOyj= - 
¢ {0 (Aj), 
Siaeeestress 
EF strain 
*: instantaneous compressibility 
ks: statical compressibility 
ct: relaxation time for shearing 
t’: relaxation time for volume change 
G: rigidity 
4: viscosity. 


If the material in question is solid, we may 
take 


b) Damping and Dispersion of Visco-Elastic 
Waves. 


According to Topa’s model, the velocity of 
propagation of the longitudinal wave Vp is 
given by : 
1224-6 wt? 
Oke 0 lee 


V,?= 


where p is the density and w= Sh, T being 
the period of the wave. For w=0, the velo- 


city is given by 


1 
a= (a — Vp, shelisel eletereie GE); 
while for w=, it is given by 
zh 4G) ; 
Vp= aa Se San V fica. cctete (Lye 


Generally for «© >w>0, Vp has of course 
a value between V,. and V,». Ke, G, Tt, 0 
are constants, if the medium is homogeneous. 
As is seen from (1.6), the velocity of propa- 
gation of the longitudinal wave varies accord- 
ing to its frequency, or in other words, the 
medium is dispersive. 

The velocity of propagation of transverse 
waves Vs is given by 


which also depends on @. 

When a wave is propagated spherically in 
an isotropic homogeneous medium, the am- 
plitude of the wave decreases with’y according 
to the following relation, 


yrAe eo .. e, (1.9), 
where 2: distance from the origin of the 
wave 
y: amplitude of the wave 
A: a constant 
k: attenuation coefficient 


At): arbitrary function of*time ¢. 
If the factor e~* in the above expression is 
due to the viscosity of the medium in which 
the wave is propagated, then according to 
Topa’s model, we can obtain the following 
approximate relation between k and the visco- 
elastic constant, 


_ 2° 1a 
_ kas v4) oe (1.10). 
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2. Decrease of Amplitudes of Initial Motions 
of Near Earthquakes. 


The change in amplitudes of initial motions 
of near earthquakes was studied by H. 
Honpa (1931), in connection with the North 
Idu and Ito earthquakes. The same materials 
will be used in this paper, for the purpose of 
studying the damping of intial motions of 
near earthquakes. 


As the hypocentres of these earthquakes 
were very small in depth, one can suppose 
that the horizontal component H of the initial 
motion observed at a station is approximately 
a function of the epicentral distance 4 as 
well as of its asimuth with respect to the 
epicentre. Thus we may put 


As vA eee ae (2.1): 


Now supposing that the function F can be 
represented as the product of f(4), a function 
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Fig. 1. Amplitude-distance and time-distance 


curves for the North Idu earthquakes. 
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Fig. 2. Amplitude-distance and time-distance 


eurves for the Ito earthquake. 


only of 4, and g(¢), only of ¢, we put as 
follows : 
Ne LON AG NG. CARB ob oe (Zz. 


where C is a certain constant. 

In view of determining the functional form 
of f(4) at first, assume that the seismic body 
wave is propagated spherically through the 
earth’s interior. The amplitudes of the initial 
motion of horizontal component H as well as 
of vertical component Z are proportional to 


ae for a given azimuth. For a given 


distance from the epicentre, the amplitudes 
vary according to the azimuth. A (or Z) 
observed at a station is thus represented by 


H(or Z)=C5e SHO Ae 0 cos (253) 


where C is a constant, which may have dif- 
ferent values in different pairs of quadrants 
into which the area around the epicentre is 
divided according to the direction of the 
initial motions. ¢ is the azimuth of the station 
measured from the epicentre with respect to 
a nodal line which divide the area into 
quadrants. kis the attenuation coefficient of 
the medium, and 4, H, ¢ are values which 
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can be obtained from the data concerning 
particular earthquakes. k can therefore be 
calculated. ; 

Fig. 1 and 2 were obtained by taking 


log sno in ordinate and J in abscissa. 


The points in these figures are distributed ap- 
proximately on three straight lines, of which 
the slopes represent the value of & each. 

It is likely that the three straight lines 
correspond respectively to the three branches 
of the time-distance curves of P,, P* and P,. 
So each of the attenuation coefficients may 
be considered to correspond respectively to 
the grantic, intermediate, and _ peridotite 
layers. From Figs. 1 and 2, the attenuation 
coefficients for the three layers are 
Ky=7.9 x 10-7 em. -1 
K*=VSx10—* erm. 1 
Kn=2.6 x 10-8 em.-1 


Granitic layer 
Intermediate layer 
Peridotite layer 


In the cases of the Fukui and Tottori earth- 
quakes, also, the similar relations were found. 
By assuming V,y and Vz as follows: 


Vp km./sec. | Vg km./sec. 
Granitic layer 5.5 Se2 
Intermediate layer 6.15 4.0 
Peridotite layer 7.5 4.5 


the relaxation times of the three layers were 
obtained from the equation (1.10) as follows: 


Table I. 
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Granitic layer tT 7==0.46 see. 
Intermediate layer t*=3.5 
Peridotite layer Tm =16.5 


3. Dispersion of Earthquake P-Waves in the 
Earth’s Crust. 

According to the reports concerning with 
artificial explosion experiments which have 
been carried out at various places of the world, 
the propagation velocity of the longitudinal 
wave through the granitic layer is about 
6.0 km/sec. This value is sensibly greater 
than 5.5 km./sec., which is the value obtained 
from the analyses of the time-distance curves 
for near earthquakes. 

GUTENBERG (1951) and Roru# (1950) gave 
the velocities of P-waves as listed in Tables I 
and II. Besides, GuTENBERG stated, ‘‘It seems 
more likely that the velocity of longitudinal 
waves below the sediments is about 6 km./sec., 
increasing to 6.5 km./sec. or more at a depth 
of 10km., and possibly decreases below a 
depth of between 10 and 15km. Such a 
decrease in velocity is to be expected in rocks 
with an appreciable content of quartz, since 
in laboratory experiments a decrease of elastic 
constants of quartz with increasing tempera- 
ture has been found approaching the tempera- 
ture at which transformation from alpha- to 
beta-quartz occurs.’’ 

But, assuming that the earth’s crust behaves 
as a visco-elastic body for waves propagated 


Velocities of Longitudinal Waves in Continental Areas 


from Artificial Explosion (GuTENBERG). 


senses 


Region Reference Py, km./see. P* km./sec. Py km./sec. 
Northern Germany WIECHERT 5.98 6.7 _ 
Northern Germany BROCK AMP 5.9 — 8.18 
Helgoland WILMORE 5.95 6.5 8.18 
Helgoland SCHULTZ and FORTSCH 5.34 6.55 S23 
Black Forest REICH et al 5.9-6.0 6.54 8.15 
Black Forest ROTHE et al 5.97 6.45 So 
Canadian Shield HopGson 6.15 — 8, 
New England LEET 6.0 6.7-7.1 8.15 
Near Washington D.C. TUVE 6.0-6.7 = a= 
Southern California Woop and RICHTER 6.0 — — 
Near Corona, Calif. GUTENBERG 5.7-6.0 6.5-7 8.1-8.2 
Near Corona, Calif. TUVE and TATEL (Has! 6.8 — 
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Table II. Velocities of Longitudinal and Transversal Waves from 
Natural Earthquakes and Artificial Explosions (RorH®). 


Natural Earthquakes 


Region Reference P, km./sece. S, km./sece. 
Japan MATUZAWA 5.0-5.5 ~ og 
Garfagnana ROSINI 5.25 3.06 
Mississipi JOLIAT 5.4 3:3 
New Zealand HAVES 5.4 3.4 
England JEFFREYS 5.4 3.3 
Dinarde Tables MoHOoROVICIG 5.54-5.59 — 
Isére ROTHE 5.55 = 
Bodensee HILLER 5.55-5.63 3.d1-3.02 
California GUTENBERG 5.57 3.23 
Tables JEFFREYS 5.57 3.36 
SW Soltau WILLMORE 5.6 = 
Southern Germany GUTENBERG 5.6 3.2 
Central California BYERLY 5.6 3.26 
Northern Tyrol GRAFE 5.7% 3.5 
Cansiglio CALOI bai; 3.86 
Alpes Carniques CALOI 5.7 3.3 
Brabant GEES 5.7 3.46 
Belgium SoMVILLE 5.95 3.5 
Houtes, French Alps ROTHE 6.0 3.5 

Seismic Explosion 

Helgoland N.O.L. 5.9 a 
Haslach RoTHE 5.97 3.42? 
Haslach REICH 6.02 = 
California GUTENBERG 6.0 = 
New Zealand LEET 6.0 3.5 
Eastern U.S.A. TUVE 6.0-6.17 = 
Northern Germany WIECHERT 6.0 — 
Soltau, Profile E WILLMORE 6.0 = 
Soltau, Profile W WILLMORE 6.3 wr 


through it, the dispersion phenomena may 
occur. From the equation (1.6), the velocity 
of the visco-elastic body wave is a function 
of the period and visco-elastic constants of 
the medium. 

The visco-elastic constants, tc, G, t and 0 
of each layer of the earth’s crust are not 
available. The periods of waves caused by 
explosions and natural earthquakes differ as 
follows : 


Explosions about 0.1 sec. 
Natural earthquakes about 1-2 sec. 


Using the equation (1.6), we get Fig. 3, 
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which shows the relation between - = and 


1+0?7? 
w, taking t as parameters. From Fig. 3, 


if T=0.1sec., whether the value of t is 0.45, 
or 3.5, or 16.5sec. It is by this reason that the 
propagation velocity of longitudinal waves 
through the earth’s crust caused by explosions 
is equal to V»., while that by natural earth- 
quakes is variable according to period. 
Now, if we write the longitudinal wave 
velocities in the artificial explosions and in 
ear earthquakes by Pexp:. and P seis. respec- 
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—- wlsec™!) 01 


Fig. 3. 


where M= 


tively, we get from the equations (1.6) and 
(1), 


“Thy 4G. ij 
Pay ee . a . Se. sere “- - -(3.D, 
1 4G wo??? } 
P seism.= Oks Sy ol eartes €3-2). 


For calculating the values of Pseis. for the 
granitic layer, we take Pexp).=6.0 km/sec and 


/ c eaiflonjcee. \Thetcalculated) yalues: ot 


P seis. are shown in Fig. 4 for t,=0.1-1.0 sec. 


If it is preferable to take /< =3.2-4.0 km/sec., 


P seis.. for t7=0.45 sec. are as follows: 
P seis. =5.85~5.75 km/see for T'=1 see 
=5.71~5.39 km/sec for 7 =2 sec. 
These results agree with the observed facts. 
On the other hand, in the intermediate and 
peridotite layers, the relaxation times are 
t*=3.5sec. and r,=16.5sec. respectively. 


Onn 
ST 1 for T=1-2sec., 


Since from Fig. 3, - 
we see that 

P seim. == WAS = 
The above considerations may explain why 
there is hardly any difference between the 
longitudinal wave velocities caused by explo- 
sions and earthquakes in the case of the inter- 


mediate and peridotite layers. 


as 1 
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Fig. 4. Peis, curve for T=1sec. and T=2sec. 


1.0 


4. Viscosity and the Relaxation Time of the 
Earth’s Crust. 


The idea of viscosity is somewhat different 
according to the different models of visco- 
elasticity. As stated in section 1, there are 
many such models. Here, in this paper, the 
approximate values of viscosity are calculated 
adopting the models due to Vorct, MAXWELL 
and Topa. 
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We can get the relation between the inter- 
nal friction or the firmoviscosity and k, with 


special reference to the P-waves. We have 
pall DAY AP a er 4 
Rk Gee Vp peu PRCA NO a EC (4.1) 
“‘h «1 A+2p 
k=5 Wp tte (4.2), 


where (A’+’) is the internal friction, and » 
is the firmo-viscosity. The viscosities for 
different models may be obtained. With the 
values of 


| 
Layer Aerie (km, See.) (em--») e | ie 
Granitic| 5:5 | 3.2 |7.9x10-72.8| 1 
Bee | Gelb 4.0 1.8x10-7 3.0| 1 
ee 4.5 (2.6 10-8 3.4; 1 


we get the calculated values of viscosity as 
follows : 


Layer | 24! +2u! (pois.)| v (pois.) | 7 (pois.) 

E | | 
Granitic 1.9x10 | 9.8x104 | 18x 10u 
Intermediate) 4.6x10" | 7.1 10" | 1.7 10 
Peridotite | 19x10" | 4.910 | 1.1108 


GUTENBERG and ScHLECHTWEG (1929, 1930) ob- 
tained the following values: 


Mantle 5x102 poise 
Surface 1.5x109 7» 
Sediments 107» 


while Sezaba and Kanar (1940, 1941) obtained 
the following values: 


Depth (km.) 
0~70 10° ~—s poise 
70~100 100~10" 
1000~3000 10° ” 


On the other hand, Scaweypar and others 
obtained the value of viscosity of the earth’s 
mantle, using Maxwetu’s model which is 
about 10”? poise. 

Next, we shall try to determine the relaxa- 
tion*times for the granitic layer, by using 
pe values: Of Pexcin Pi tig, and Siccis..! The 


VAlUe Of Peasy S ceewand P exp!. are obtained 
as follows. If we assume that the period of 
Peis, is equal to that of S seis, we get from 
sections 1 and 3, 


Pte ete 


Piss 3°90 
1 4G wo 
FOIE pate fy nea SI nee ond Hee 
seis. Oks 2 3 p 1+u27? -- (4.3). 
G wt? 
S seis. = SF EE 


From the equation (4.3), the relaxation times 
of the granitic layers at various places can 
be determined. But, unfortunately, the values 
Of Pseis., Sseiss and Pexp, at the same place, 
are not avilable. In this paper, by assum- 
ing Pexp1.=6.0 km./sec., ct, is determined. 
These values come out to be 0.36~0.645 sec. 
and this is shown in Table III. These 
values agree with those determined from 
damping. 


Table III. The Relaxation Times calculated 
for Various Places in the World. 
nee eee Ce 
Japan B-5.5.| 3.25 36-.485 
Garfagnana 5.26 | “8.060 || 489 
Mississipi 5.4 im os-3 .465 
New Zealand 5.4 | 3.4 .465 
England B.A 8.3 .465 
California 5.57 3.27 .53 
Tables 5.57 3.36 -53 
Southern Germany 5.6 3.2 -552 
Central California 5.6 3.26 .53 
Northern Tyrol 5.7 Bhi -635 
Alpes Carniques 5.7 | 3.3 64 
Brabant 5.7 8.46 .645 
Belgium 5.95 3.5 1.46 


5. Some Rheological Properties of Alluvial 
Layers and their Relation to Earthquake 
Damages on the Layers. 

Similar studies as above can be made for 
the cases of seismic prospecting on the earth’s 
surface, such as have been carried out at 
various places by the members of our Insti- 
tute, in view of determining relations be- 
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tween the physical properties of the earth’s 
surface with earthquake damages experienced 
on it, (Basa, Hatanaka and Fusimoro, 1949 ; 
Hatanaka, 1949; Kamerant, 1949; Koporera 
and YAMANAKA, 1949; Kusotera, 1950; 1951 ; 
Momosz, 1950; Yamapa and Kawanam1, 1949 ; 
Yoxoo and Basa, 1949; YosHrkawa, 1949, 
1950). 

According to these studies, the bed rock 
(Vp=2.0~3.0 km./sec.) is from 10m. to 100 m. 
in depth and is overlaid by an alluvial layer 
(Vp=0.6~1.8 km./sec.). Surface waves which 
are of RAYLEIGH type (Kusotmra, 1951) are 
propagated along the surface. The propaga- 
tion velocities (Vz).of these waves are about 
0.1~0.5 km./sec.. The ratio of Vp/Veg was 
found to be closely connected with earthquake 
damages on the surface. 

However, it is more important to determine 
the relation between other properties of the 
earth’s surface, specially Rk or rt and earth- 
quake damages, 

The values of kp and kg at Wakayama are 
as follows (Kusorrra, 1950): 


A 
fo) 


CONRAN WrOH | 


| Ap (in 10-4cm.-1) || ke (in 10-4 em.-1) 


Nowak 
acorns 
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Here, kp and kz are the attenuation coeffici- 
ents for P- and Rayleigh-type waves. kp» and 
ky are roughly proportional to the earthquake 
damages or Vp/Vz. These relations are illus- 
trated in Figs. 5, 6, 7 and 8. 

The relaxation times obtained from the 
equation (1.10) by assuming Ve=+Vs, are as 
follows: 


| t (10-3 sec.) 
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The relation between r and V»/Vz is illustrated 
in Fig. 9. From Figs. 5 and 9, it is seen that — 
the relaxation times of the material of alluvial 
layers are closely connected with the earth- 
quake damages experienced on them. 

In general, when zc decreases, the medium 
becomes more fluidal. It is therefore likely 
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Fig. 5. The relation between damage per- 
centage of dwelling houses and at- 
tenuation coefficents (in 10-*). 
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Eig. 6. The relation between damage per- 
centage of dwelling house and Vpl Vr. 
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Fig. 7, The relatien between kz (in 10-*) 
and V,/Vr. 
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Fig. 9. The relation between < and V)/Vp. 


that earthquake damages are large, where 
the alluvial layer is fluidal. The popular 
statement that “‘the earthquake damages are 
large, when the alluvial layer is weak’’, can 
reasonably be understood in this manner. 
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Changes in the Equilibrium Position of a Sensitive Torsion 
Balance Placed at Inuyama in Relation to Flood Flows 
in the Nearby Kisogawa River. 


By 
Naomi Mryasp, 
Geographic Survey Institute, Chiba, 
and 


Ikuei Muramatsu, 
Gifu University, Gifu. 


Abstract 


A sensitive torsion balance was installed at a place of Inuyama town and continuous 
observations have been made of the change in its equilibrium position. The instrument 
was placed at a point on the bank of the Kisogawa river, which flows through the 
town, with the intention to investigate the change in the gravitational potential around 
the instrument caused by the rise and fall of the level of the river water, especially 


at the times of flood flows. 


Caleulations have shown that the observed deflections fall much short of those 
expected from the visible changes in the level of the river water, thus suggesting there 
are notable changes in the underground water conditions associated with the floods. 


The equilibrium position of a torsion balance 
is determined by the state of the gravitational 
field in which it is placed. So, if a change 
occurs in the field itself, a corresponding 
change in the equilibrium position of the 
balance is expected to be observed. A very 
sensitive torsion balance was constructed and 
it was installed in Inuyama town, Aichi Pre- 
fecture, at a place on the bank of the Kisogawa 
river. The equilibrium position of the balance 
has been continuously recorded photographi- 
“cally (Muramatsu, 1950). It was our intention 
to detect the change in the gravitational field 
around the instrument caused by the water 
mass of the occasional flood flows of the 
river, and from the observations to get infor- 
mations, if possible, regarding the conditions 
of the underground water in these cases. 

Some portions of the actual records obtained 
are shown in Fig. 2, together with the heights 
of the water level of the river. As may be 
seen from Fig. 2, the deflections of the torsion 
balance are evidently associated with the rise 
and fall of the water level of the Kisogawa 


river as expected. At the times of the flood 
flows of the river, on March 7, April 6 and 
June 20 of 1950, the deflections of the torsion 
balance were especially notable. 


We Scale for Map 
0 30 60 90 120M. 


Topographic map of the vicinity of 
the observation room. 


Big. is 


The torsion balance used in our observations 
is of the type of inverted T, so that its de- 
flections are caused by the change in the 
curvature terms of the gravitational equi- 
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potential surface. 

For a torsion balance of the type here 
employed, the deflection which is caused by 
a change 3U in the gravitational potential is 
given by 


s9-2 il ee OU es: oO” 6U sin Or 
TL2{ Oy Ox? 


where J is the moment of inertia of the sus- 
pended beam of the balance around its central 
axis, and t is the,restitutive torsional moment 
of the suspension wire given by 

n ura 


ae 


a is the angle between the direction of # and 
the suspended beam. 
The gravitational potential dU due to an 


elementary mass pdv situated at a distance 


ry from the instrument is given by 
sua Bee 


Now the torsion ee was installed on 
the bank of the Kisogawa river in such a 
way that the suspended beam of the balance 
is in the direction 45° inclined to the general 
trend of the course of the river in the neigh- 
bourhood. The axis of 2, being taken in the 
direction parallel to the general trend of the 
river and that of y in the direction perpendi- 

6U vee) Apr 
Oy" 
comes into our question. The origin of 2, y 
is of course taken at the point of the torsion 
balance. 

An elementary water mass pdv situated at 
a distance 7 from the torsion balance in the 
horizontal direction m with reference to the 
x axis will produce a deflection of the torsion 
balance which is given by 


cular to it, the term ( 


’ SO) al 
-5 Fo cos Se. dv. 


Hence we have 
3 cle 2 
d0=—5 s pip | 20820 dv. 
By putting dv=rdpdrdh, we get 


s0=—% x0) cos apap | = "| an. 
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Actoal and calculated deflections of 
the torison balance and the height of the 
water level. 


By the customary graphical method, the values 
of 6@ can be calculated if the profiles of the 
river are given. 

The topographical features in the neigh- 
bourhood of the observation point is shown 
in Fig. 1, from which the profiles of the 
sections of the river were obtained. The 
expected deflections of the torsion balance 
were calculated in terms of the height of the 
water level of. the Kisogawa river, the beeans 
being shown in Fig. 3. 

The deflections of the torsion balance 
observed at the times of the discharge flows 
of the Kisogawa river are compared with the 
corresponding heights of the water level of 
the river in Fig. 2a,.b, and c. The curves 
drawn with full lines in these figures show 
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the actual deflections of the torsion balance. 
The heights of the water level are shown by 
broken lines while the calculated deflections 
by dotted lines. The difference between the 
calculated and observed deflections is also 
shown by thick dotted lines in the same 
figure. 


Deflection 
mm 
/0 | 
8+ +—-+ | | oz | 1 
6L 
; 
a 2 
2+ a! | 
0 = 
ra) 7 ee Fo 5) Os 
Height of Water level 


Fig. 3. The relation between the deflections 
of the torsion balance and the height 
of water level of the Kisogawa river. 


The difference between the observed and 
calculated deflections is at its maximum at 
the time when the height of the water level 
reaches its maximum and tends to zero as 
the water level recovers its original height. 
This fact may suggest that about the time 
when the discharge flow of the river is at its 
maximum, the water content in the soil near 
‘the observing point increases considerably. 

0’?6U 


In addition, the influence of the term - aay 
may not be negligible either. 

It must be mentioned that the height of the 
. water level of the Kisogawa river was mea- 
sured from February to April 1950 at Unuma, 
about 1km down along the river from the 
torsion balance station, and from May to 
September of the same year at Inuyama-bashi, 
about 200m down from the station. 

The records show fluctuations in deflections 
of the torsion balance which are more or less 
periodic. The fluctuations were studied with 
reference to the data obtained from February 
to May 1950. 

Periodogram analyses have shown that, of 
the variations with various assumed periods, 


the diurnal and semi-diurnal ones are of the 
largest amplitudes. The variations with 
shorter and longer periods are smaller in 
amplitude and fluctuating over large ranges. 
The hourly values of the deflections of the 
torsion balance are subjected to Fourier 
analysis. The results obtained for the assumed 
period of T=24h is 


= —0.002,-+0.186 cos ee (¢—3) 


4.0169 Cos 2 


24 (t+0.97) 


27 
+0.016 cos 34 (t—2.5), 


whereas the result obtained for the assumed 
period of T=24,84h, i.e., lunar diurnal period, 
is 


— e an — 
y=0.123+0.080s cos 94 gq (¢—19-1) 


4+-0.096 cos 2-57" -(#4+2.3)+... 
The values of y are given in Eotvos, i.e., 
LO RING ors 

As to the cause of the diurnal and semi- 
diurnal variations of the deflections of the 
torsion balance, however, nothing can be said, 
except that the influence of the tidal force 
alone cannot explain the above mentioned 
periodic variations. 

Since June 1950, the secular changes of the 
mean values of the deflections of the torsion 
balance became rather notable. It was noticed, 
however, that these secular changes may be 
due to the gradual bending of the suspended 
beam by the weight attached to its both ends. 
This partial defect in the structure of the 
instrument is now being repaired and the 
observations are expected to be revived soon. 

In conclusion of this note, the writers wish 
to express their thanks to the Government 
Authority of the Ministry of Education for 
the grant of the Science Research Fund by 
the aid of which the present observations 
were made possible. 
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or less irregular fluctuations. 


involving winter. 


In the present paper are given the results 
of analysis of the materials concerning the 
vertical movements of the surface of the 
ground in Marunouchi and Hibiya, in the 
central part of the Gty of Tokyo. The 
vertical displacements were deduced from the 
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Notes on the Vertical Deformations of the Ground Surface 
in the Central Part of the City of Tokyo. 


Abstract 

1. The vertics] displacements of the ground surface in the central part of the city 

A Tokyo, which have been detected by levelling surveys repeated recently, consist of 

two components, namely, the sinking due to the shrinkage of the soft superficial soil 

‘layer and the displacement resulting from tilt of the underlying layer of harder 
material. The former is progressive as 2 whole, while the latter is quasi-periodic with 

respect to time. The vertical displacements measured for longer intervals of time are 

proceeting in one direction, although those messured for shorter intervals show more 


2 As 2 whole, the rate of sinking has been evidently decreasing. 
3 The rate of sinking during the intervals involving summer is larger than that 


fixed to the wall of several selected buildings 
in the zrea, 28 well as to several other. ordi- 
nary benchmarks. The locations of the 
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Aa Part 


Fig. 1b. Reference to the area under 
consideration. 


40 Naomi MryABE and Akio KATAOKA. 


Table I. Vertical Displacements of Bench-Marks in Hibiya 
and Marunouchi Area (in mm.). 


en Uime May . | April : P Ov. ay 
Interval) 1938 | 1938 | 1939 1939 | 1939 | 1940 | 1940 | 1942} 1942 | 1943 | 1943 | 1944 
B.M. \ Nev. Feb. Al ah Nee eee oct. Apel Sedt. Feb. | Kae a nev. 
| 1938 | 1939} 19389] 1939| 1940] 1940| 1942| 1942| 1948| 1943 | 1944 | 1944 
Miya 3 — 3.6) +38.5 0 +2.1]| — 1.7 
Huku 3 14.8] —.7.6| — 7.8 +0.9 | -2.2 | +8.7| — 2.2 
KW — 1.9) — 0.4) +1.1] — 3.6 4.8) — 6.2 2.0; — 2.6) +3.3 | +0.4 |} +0.3] — 1.5 
KE — 1.2) — 1.2) +1.7] — 3.9 4.6) — 5.8 2.0| — 2.7) +3.4 | +0.4 | =-0.3 | — 1-6 
Miya 2 — 12:4) 0.84 6ri) St Bade ees, 
Miya 1 — 4.7) +3.9 | -1.2 | 41.0] — 0.1 
Kozi 2 — 2.2) — 3.9| +2.8| — 4.4) + 3.3) — 8.9} + 3.3) — 5.5) +5.6 | —0.7 | —1.5] —19.5 
M21W — 6.2} — 6.3) +1.7| — 7.7) + 2.7) — 9.6] — 2.9) — 6.6) +3.0 | —0.6 | 4+1.0°| — 2.6 
M21E = 3.8) = 6.1) +2.4 | — 8.5) + 1.6/ —10.0) = 2.9) = 5.3) 41.8) —022 | +140) = 384 
M21M — 6.2) — 6.7) +2.3 * 7.4) + 2.0) — 9.5) — 2.1) — 6.0) +2.2 | —0.1 | 41.1) — 2.7 
Kozi 3 | —22.9) — 8.4/°+0.2 | —15.3) + 5.4) —21.4) — 3.2) .—12.4) +7.7 | —3.2 | +6.4] 4+ 0.1 
A —27.5) —11.9) —0.3 | —20.3} + 3.6) —24.6 9.7) —14.5) +6.0 | —7.1 | +6.4] 4 0.6 
HE — 9.0) — 2.6) +3.4 | — 8.3) + 6.3) —12.3) + 6.7) — 6.4] +6.4 | —1.1 | 44.3) — 2.0 
HW = 8.9) = '2°6) 423.8 | — 8.5)\-- 6.7) 1227) 2 7 a= Goal -EGea Osa ear Guam 
Tetu i = 9:8) — 4.9) — 3.3} £5.0 |) —O°3 "E16 te 
Tetu 2 —12.8) + 5.3) = 6.5) te) — Ort Ok a eaters 
Tetu 3 —lf.2) 4 82 968) 28-5 1504) eeeken andes 
Tetu 3/ —14.4) + 5.3) = 8.2) +7.5 | —0.9 |. 4122) =e 
Tetu 4 —13.1) + 5.9) —97-3) 7.0 | —0.7 1210 eae 
Tetu 5 — 9.9) + 6.0| — 5.5) 45.5 | =—0-8-| 41.8) 2 aa 
OE — 9.4/ — 2.9) +1.3 | —18.8] +11.0) —15.9} + 8.1) — 9.7] 49.9] —9.8 4+1.9] — 1.8 
OW —138.2) — 2.0) +1.0 | —18.7) +10.5) —17.3) + 4.9) +10.0]) +8.8 | —2.1 +1.9 | — 0.9 
OM —13.2) — 2.1) +1.8 | —16.0} +10.6) —17.2) + 6.1] — 9.4] +8.4| -1.7 +2.6 | — 2.2 
SE — 6.8) — 2.5) +4.3) — 9.8) 410.3) —12.0] + 6.8) — 6.7 +670 |) — 059) | =e 24a one 
SW, — 8.2} — 8.8 +4.5 | — 9.6] + 9.5] -10.5] + 6.81 — 4.5) +3.9| -0.9| +2.6| — 4.9 
SW, - Go) — 8.8) +8.7 | —. 8,2) 469.8! —10.8).4 6:1! — 6.4 +6.2 | —0.9 | +2.85)— 4:5 
Pp’ —40.0) —14.6) —8.1 | —37.6) — 9.8) —33.9] —31.9] —19.7 +0.3 | —8.3 | +4.0] — 3.8 
py —27.3) — 9.8| —6.8 | —47.6] + 2.5) -81.7| -17.5] -15.3| +3.3| —7.4] — s 
Kozi 4 — 0.3) — 2.9) +2.0 | = 2.1) +.2.9) = 4.8) + 4.9). 27 +155") —0°6 190541 = a4 
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bench-marks in the area under consideration 
are shown in Fig. 1. The precise levels have 
been repeated over the area twelve times 
since May 1938. 

The vertical displacements of the bench- 
marks found are given in Table I, the first 
row of which shows the time interval during 
‘which they occurred. 

In a previous paper (Miyazse and Huxkunisat, 
1941), one of the present authors has already 
discussed the vertical movements of the 
surface of the ground in this same area which 
took place during 1938-1940. In the present 
paper, the study will be extended referring 
to the data given In Table I. 

At a glance on Table I, we notice the 
following facts. 

(1) Some of the time intervals are charac- 
teristic in that the vertical displacements of 
bench-marks which occurred during them are 
mostly positive. 

(2) The natural sinking of the buildings 
are very small. This is especially true for 
the buildings of Toho Theatre (now Earny 
Pyle Theatre), Shisei Kaikan and Tokyo 
Marine Insurance Co., to which the bench- 
marks HE, HW, SE, SW:i, SW:, KW and KE 
are respectively fixed. 

The foundations of most large buildings in 
this area are constructed in the manner as 
described below. After the ground surface 
is first excavated down to the depth of several 
meters, piles are driven into the depth of 
about 20 meters. The piles are usually con- 
crete rods or pine logs, about 20 meters long 
and 16-20cm in diameter. The base floor of 
“most buildings is constructed on so piled 
foundation. . 

Since the soft alluvium soil layer in the 
area under consideration is about 20 meters 
in thickness, it is hardly practicable to set 
up the bench-marks directly on the hard and 
stable material which underlie the soil layer. 


The bench-marks situated in this area are | 


mostly set up on the surface of the soft soil 
layer, except BM. Kozi-4, which is believed 
to lie directly on the harder material, for, at 
this point, -the thickness of the soft surface 
soil layer is almost zero. The piles of the 


foundations of the buildings, on the other hand, 
probably reach the layer of the harder material 
underlying the surface soil layer. The build- 
ings, the foundations of which are constructed 
as described above, are believed to be con- 
nected directly to the stable and harder layer. 

The vertical displacements measured as 
have occurred during a certain time interval 
consist of two components, namely, a) the 
sinking of the surface due to the shrinkage 
of the soft soil, of which the overlying layer 
is composed, and b) the displacement resulting 
from tilt of the underlying layer. While the 
latter is quasi-periodic with respect to time, 
as will be shown later, the former is progres- 
sive as a whole, although the vertical displace- 
ments measured for relatively short intervals 
of time may occasionaly be positive. The 
relative up- and downward movements of 
buildings, on the other hand, are likely to 
have almost nothing to do with the shrinkage 
of the soft soil layer. 

In order to obtain the magnitude and azi- 
muth of tilts of the ground surface, the 
values of G are plotted against @ for each 
bench-mark. Here G is the relative vertical 
displacement of a bench-mark with respect 
to BM. Kozi-4 divided by the distance between 
these two bench-marks. G will give therefore 
the component of tilt of the ground in the 
direction connecting the two bench-marks. On 
the other hand, @ is the azimuth of the line 
linking BM. K6zi-4 with the bench-mark for 
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_ 2a. G-@ diagram for displacements 
during June 1939—Nov. 1939. 
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which G is calculated. @ is measured in 
reference to N. Two example of the plots 
are given in Fig. 2. 


G 
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Fig. 2b. G-6@ diagram for displacements 
during Oct. 1940—Apr. 1943. 


If the soft surface soil does not shrink and 
the earth’s surface is as a whole tilted in a 
certain direction, a trigonometrical relation 

G=G,, COS (9—6 mn) 
must hold, where G,, is the magnitude of tilt 
and @,, is its direction. But in the case in 
which the vertical displacement is mainly 
produced by the sinking of the surface of the 
soft soil layer due to its shrinkage, the above 
relation may. fail. 

If the vertical displacements of the buildings 
of Toho Theatre, Shisei Kaikan and Tokyo 
Marine Insurance Co. are due to the tilting 
movements of the layer of material on which 
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8: 42 Apr— 42 Sept. 
9° 42 Sept— 43 Feb. 
10: 43 Feb.— 43 Nov 
71: 43 Nov- 44 May 8 YY 
72: 44 May— 44 Nov 


Fig. 3. Tilt vectors for separate time interval. 


they are built, a» sinusoidal. curve can be 
drawn in Fig. 2 through the points that corres- 
pond to the bench-marks fixed to these build- 
ings. The sinusoidal curves in Fig. 2, which 
are drawn in this way, give the values of 
G,, and 0,, of the tilts during the correspond- 
ing time intervals. The values of G,, and 
6,, for other time intervals were-also deter- 
mined similarly. Tilt vectors thus obtained 
corresponding to each time interval are shown 
by arrows in Fig. 3, while the successively 
integrated tilt vectors are shown by arrows 
in Fig: 4. 


© 
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Pe ee? a 
@ /938-May @ 
@ 38-Nov © 
@ 39-Feb. (©) 
@ 39-June @ 
@ 39-Nov @) 
© 40-Apr @ 
RY 
Fig. 4. Integrated tilt vector. 


In Fig. 3, we notice that, — 

1) the tilts during winter time, i.e., the 
time interval from autumn to next spring, is 
in general directed to NW or away from the 
neighbouring sea, whereas tilts during summer 
time, i.e., the time interval from spring to 
autumn, is directed to SE, or towards the sea, 

and that 

2) the vectors are very small in magnitude 
if their directions are deviated from NW-SE. 

Such. a behavior of tilts like this suggests 
that they may be due to the load of the 
annual tides of the sea water in the neigh- 
bouring Tokyo Bay. 

Naturally the load of the sea water due to 
the diurnal and semi-diurnal tides has greater 
effect on the tilt of the ground surface than 
that due to the annual tide. The effect how- 
ever is disregarded in the present study owing 
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to the reason that the measurements of the 
heights of the bench-marks in question by 
precise levellings take more than a week, so 
that the influence of the diural and semi- 
diurnal tides may be much smoothed out. 
Displacement components of a point in an 
elastic body caused by a distribution of load 
applied on its surface are given by Boustnrse 
in the following form; 


Sa 0? vipa 1 0 
Arp alps ~\ 4n(a+p) ) Ox Jloee+nar, 
1 @ 1 P 
v= ——— 
Arp al  An(A+ py u) a 8 (2+7) dP, 
ae 2 | 2443 | ‘dP 
4npn Oz? 4np(Atp)) 7 


with customary notations. 

Differentiating the vertical displacement w 
at z=0 with respect to x, we have 2-compo- 
nent of tilt of the surface of the elastic 
medium on which the load is distributed. 
Thus we have 


OWo_ —At+2u 8 | dP 
Ox 4np(A+p) Ox J 7 
Putting dP=pghdS, we have 
OWo_ — A+2p a | ad 
Ox 4n(A+p) Oo” Oa | r 
AZ cos 
Deneeant ogh Hes dédr 


Vin +1 


Sin O'n, m+1) log ———, 
Tm 


At2p 
Arp (A+p 
density of the sea water, h the amplitude of 
the annual tide, and 7 the distance of an 
elementary load from the point under consi- 
deration. We have a similar expression also 


=C 3'(Sin On, m+1— 


Phere C stands for ce ogh, p is the 


For the constants in the above expression, 
the following numerical values were taken: 


MeO, 
g=9.8 x 107 cm/sec? 
h=10cm_ | 
C.2.8- 
20 


og he 15.8 x 10® c.g-s. 


o=0.38, 


The values of « and o employed were taken 
from experimental results of IsHimoro and 
Iipa (1941). The value of h is the round: 
number of the amplitude of the annual tide 
observed by the Hydrographic Department at 
several mareographic stations along the coast 


of Tokyo Bay. To obtain the numerical 


OW 


value of tilt Oe the extreme value of 7 


should be assumed. By the extreme value of 
y is meant such a distance that variations of 
sea water mass beyond it have no effect on 
the tilt of the earth’s surface at the point in 
question. 

If the extreme value of 7 is taken as 40km, 
we have 


OW: 5.06x 10-5 

Ox 
Ow 34.16x 10-5, 

Oy 

whereas if it is taken as 25km., we have 

OW -6 

Ag 3:8 10- 

OW 3 6x 1)"8. 

Oy 


The resultant vectors of these calculated tilts 
are shown respectively by dotted arrows, Cal. 
I and II, in Fig. 3. For still larger values 
of 7, the calculated directions of tilt are much 
deviated from those observed. Whichever 
of the above two assumptions is taken as to the 
value of 7, the calculated magnitude and azi- 
muth of tilts agree well with those deduced 
from the observed data. Hence, it is likely 
that the actual vertical displacements mea- 
sured by the levelling surveys are more’ or 
less modified by the elastic deformation of 
the earth’s crust caused by the load of the 
annual tide. It is also suggested from the 
fact mentioned above that the annual tide of 
sea water within the distance of about 25-35 
km from a point alone has effects on the tilt 
of the earth’s surface at the point. 

On subtracting from the total vertical dis- 
placements those due to the tidal loads, the 
residuals will represent the actual sinking of 
the earth’s surface which may be due to the 
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Table II. 
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Revised Vertical Displacements of Bench-Marks in 
Hibiya and Marunouchi Area (in mm.). 


NS ime ay z April April | Sept. May 
\ Interval] 1988 | 1938; 1939 | 1989 | 1939| 1940) 1940) 1942| 1942 | 1943| 1943 | 1944 
Bat Nev. Feb. Fis Hey. April oct, Abril Sedt. Feb: ||. Nov. bits nee 

“| 1938} 1939] 1939] 1939] 1940| 1940| 1942 1942| 1943) 1943| 1944| 1944 

Miya 3 =.2,.7 40.8") 1.8 30. ented | 
Huku 3 —11.8]-~10.9| = 6.8) —3.79) 28.3) Sie) eee 
KW = 1.2] =, 9.8|'= 0.41 — 1.6 —.1-3) a4] Lp 21.2 21.84 Stell) Se eee 
KE 20.8) 18.6) 9.8 = qa ad) — 2.1) neal = 1.0) ee en ~8.1| - 1.0 
Miya 2 = 9.8] —6.2 | —6.8 } +0.9 | — 3:0 
Miya 1 =" 2.9) 2026 | 1.77 <2.18) seed 
Kozi 2 | +130) + 7.3) = 1.2) — 2.6) = 8.1) — 0.7) — 4.8) 16) SR | 1.64 
M21W — 3.0 = 9.4] — 2.2) — 5.5] — 8.2] — 1.6| —10.7) — 2.7| —5.6 | —-1.3| =b.0| — 1.6 
M21E — 0.6) —'8.4)-— 1.7) — 6.2} — 9.9 — 1.6) —11.0) — 1.3] —7.1 | -0.9| 24.6) = 2:3 
M21M — 2.9 — 9:9) — 1.8) — 5.2; — 9.8 = 1.1) ~10.2] = 2.0] ~6.7 | —0.8| —5.2| + 1.6 
Kozi,, 3 | —18.6) -11-1]) — 4.5] ~18.4| — 7.0| 111.4] -19.31 — 7.65 -1.8| -3.38] —2.8] + 1.2 
A —28.1] —14.5| — 5.0] 18.4) — 8.7} —14.7| —18.8| — 9.7 =3.4 | —7.2 |.-0.6 | + 1.8 
HE — 4.4) — 5.1] — 1.5] — 1.4 — 6.2) — 1.1] = 2.6 — 1.6 29.9 |= ee 6 Beas 
HW — 4.4 — 6.0, — 1.0) — 1.7°— 6.1] — 2.7] — 1.8) 2:6) 23.7} LOT) sea) Sib 
Tetu 1 + 1.5 — 5.2 1.6 —6.5 | 0.2) 6 ees 
Tetu 2 + 1.4 — 5.0) — 0.9) 3.4) 20/2) 25 90 oes 
Tety 8 | - 5.4) = 6.9 — 4.0) 9.2} —6.7)). 26.6 |.— 0.1 
Tetu: ‘3! —' 2:3) ..6.5] — 8.8) — 3.40) 20.6 ee © 
Tetu 4 — 0.8), — 5.1), — 18)4.0 | 0.8) ee) ong 
Tetu 5 + 3.1) — 5.J} + 0.7 —6.7 | + 0x8 |) aad enoee 
OE — 8.6 — 5.0| — 4.5) -12.1) — 3.5} — 3.5] — 2.8] — 3.6] -0.9 | -1.7] —6.5| — 0.9 
Ow — 7.6) = 4.0) — 4.6) -12.1) —14,8| — 5.8] = 6.6 = 4.0) ~1.7 | -1.6 | 6-2 be oe 
OM — 7.4) — 4,2) — 4.2) 14.8) — 8.8) — 4.9] — 4.7] — 3.21 -1.2| -1.1| —5.8| + 0.5 
SE - L8 - 4.2) — 0.1) — 8.5] - 0.9 — 2.5) 1.6 = 2.09.71) 05) 28 ues 
SW, rh ded) — Bed] + 0.2) = 8.3 = 158 = 1.8) 2.8 0 | 429) 095425 gee 
SW, + 1 = 5.4) = 0.6 = 7.1) =F} — 1.6) 2.0) 2.3) 20.5) 20 glee eee 
Pp — 85.2, - 16.2) —12.7) —36.4) -21.3, —23.9| —40.7| -14.3| -8.5 | —7.8| -2.8| — 9.9 
py — 25.9, —11.4) —11.2} -46.4) — 8.4) —22.4| —25.7| -10.7| -4.8| -6.9| — = 
Kozi 4 | + 0.9) — 4.2) + 0.5) — 1.7) — 1.3) = 1.81 — 18] & 6.8] te 


PT kIT 
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shrinkage of the surface layer of soft soil. 
These residual vertical displacements are 
given in Table II. 

As may be seen in the above table, the 
values are mostly negative, showing that most 
bench-marks sunk during each of successive 
time intervals. 


In Table II, the sinking of the bench- 
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marks, P’, P’’”, A and K6zi-3 are especially 
remarkable. To see the general tendency of 
the sinking of the earth’s surface in the area 
under consideration, the average values of 
monthly sinking rates of the four bench- 
marks mentioned above are calculated and 
are given in Table III. 


Table HI. Average Sinking of Bench-Marks P’, P’’, A and Kozi-3. 


ay | Nov. | Feb. | June 2 
1938 | 19388! 1939 | 1939 | 19389 
Period ) 


(G06 a” | Rh OEE as ae 
Nov. | Feb. | June | Noy. | April 


mm/month 4.3 | 4.4 | 21 | 6-7 | 2.8 


1938 | 1989 | 1989 | 1939 | 1940 | 


| | Feb. | | F ay 
1940 | 1940 | 1942) 1942 | 1943 |} 1943 | 1944 
t u eal Some! Be | 
Oct. | April! Sept. | Feb. | Nov. | May | Nov. 
1940 | 1942 1942 1943 1943 | 1944 | 1944 


Apri et. | Apri Sept. 


| | 
8.0 14 2 | 0.9 | 0.7 | 0.8 |0%4 


We notice in Table III that, 

1) the monthly rate of sinking has been 
gradually decreasing, 

and that 

2) in the interval involving summer and 


autumn, the rates are larger than in the 
intervals involving winter and spring. 


In conclusion, the writers wish to express 
their sincere thanks to the Municipal Authority 


of Tokyo who has supplied the valuable data 
of levellings and placed them at our disposal. 
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Magnitude-Frequency Relation for Earthquakes in and near Japan. 


By 
Chuji Tsusot. 
Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


By means of the formulas determined by the present author previously, the magni- 
tudes of 735 shallow earthquakes which took place in and near Japan during 1931-1950 
were determined. The mean annual numbers of earthquakes having the magnitude M 


were found to be expressed by 


log N= —1.60+1.06(8—M), for the area A, 
log N= —1.57+-0.72 (8—M), for the area B, 
log N= —1.61+0.66 (8— M), for the area C, 


log N= —1.46+1.04(8— M), 


for the area (A+B), 


log N= —1.33+1.01(8—M), for the area (A+B+C), 


the class interval of M being taken as 0.1. It is interesting to note that while the 


first constant does not differ much from area to area, the second does so. 


The second 


constant for an area is nearly proportional to the total number of earthquakes which 


occurred in that area. 


In his previous paper (Tsupor, 1951), the 
present writer obtained formulas which can 
conveniently be used for determining the 
RICHTER-GUTENBERG’S instrumental magnitude 
M of an earthquake, the maximum vibration 
amplitude A of the ground due to which is 
known at a known distance 4 from its epi- 
centre. The formulas obtained are: 


M=0.20 4+0.67 log A+3.80 (4<5), 
M=0.03,4+0.60 log A+5.00; (4 >5), 


where 4 is measured in 100km and A in 
micron. Log is the common logarithm. The 
former formula is applicable to the case when 
the epicentral distance is smaller than 500 
km, while the latter to the case when it is 
greater than 500km. In determining these 
formulas, only those earthquakes which had 
occurred in and near Japan and which had 
been observed in Japan were taken into 
consideration. It is not clear therefore whe- 
ther or not their validity can be claimed 
for earthquakes in other areas. 

During 20 years from 1931 to 1950 inclu- 
sive, 735 shallow earthquakes were reported 
by C.M.O. to have taken place in and near 


Japan with such magnitudes as to be felt 
beyond 200km from the respective epicentres. 
By means of the formulas given above, the 
magnitudes of all these earthquakes were 
determined, taking as materials the seismo- 
metrical observations obtained at the follow- 
ing six principal seismological stations in 
Japan : viz., Sapporo, Sendai, Tokyo, Nagoya, 
Kobe, and Fukuoka. The observational data 
at each of the six stations yield independently 
the magnitude for one earthquake. In most 
cases these six independent determinations 
result in a satisfactorily consistent magnitude 
of one and the same earthquake but in 
several other cases they do not so. In these 
latter cases, the mean of the six values is 
adopted. 

Among the 735 earthquakes, given in Table 
I, 38 were found to have had the magnitudes 
greater than 7.0. 

In the table, Mr is the magnitude as 
derermined by means of our formulas using 
seismometrical data obtained in Japan, while 
M, is that determined by GurmnsEeRG and 
Ricater and published in their book “‘Seismi- 
city of the Earth’’ (1949). 
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Table I. 38 Earthquakes with Magnitudes 1945 | Jan. 13 38 38 187.2 | 34.7 | 7.0 t 
reater than 7.0 Feb. 10°28) 58 | 142.1 | 4029s saiiice 
3 s Apr. 10 10 22/| 143.0 | 41.8 | 7.1 
ipidenthe 2423 36 | 139.54) S10 ie0 
Bee te Oct. 91823) Se | 147.05) 4s.8 3 aie 
ede Na a MG 1946 | Dec. 21 4 19 | 135.6 | 33.0 | 7.7 
x bs 1947 | Apr. 14 16 16 | 149.2 | 41-27) 7-0 
1981 | Feb. 20 14 34 | 188.0 | 44.0 | 7.0 Oct. 4 9 9/| 141.0 | 48.8] 7.1 
Mar: 9. 12-49. 141.9 | 40.6 7.2 |. 7.7 1948 | June 15 20 44) 185.5 | 33.8 | 7.2 
1932 | Nov. 13 13 48 ee a ae We0) ap 28 16 13) 136.2 | 36.1) 7-2 
1988 | Mar. 3 .2 31 | 144.7 | 39.1] 7.7 | 8. 21 | 148.8 | 46.0 | 7.4 
June 19 °6687'| 142.8 (88.6 | TAAR de ee 
1935 |} Sept. 11 23 4] 145.1 | 42.7 | 7.0 ie 
Oct WS 45) 14828) | 4024 720 | 73 2 
: 18 9 12] 143.8 | 40.2| 7.0 | 7.2 The frequencies of earthquakes are tabulat- 
19386 | Nov. 3.5 46 | 142.0 | 38.4 | 7.4 | 7.2 : : e 
1937| Feb. 21 16 3| 160.0 | 44.6 | 7.5 | 7.4 ed in Table i aCcOTdIaE to the magnitude 
drotlhss ae ZI Gaal) DONE || ea and the year in which they occurred. The 
VOSS aye elcome Oe IRL Sh al sock iO Gao Z 2 i rth | 
May. 28 16 18 || 141° |/86.7 | 7.27) 7.4 frequencies are given separately for ea 7 
Oct 129 9) 8b) 14430) 3928) 720 i ; 
Now Cee iy de ita te Ll ea ee quakes es — agit’ a the eS A, and 
5 19 50/141.7 | 37.2|7.0|7.7  B,andC as indicated in Fig. 1. This group- 
: x Me ee Ae Ue Oe ing of epicentral areas is not arbitrarily made. 
1939 ae ei Ne a jegie a oe i. On the contrary, it is suggested naturally by 
1940 Aug. 2 0 8)| 189.1| 44.3|7.3|7.7 the characteristically uneven distribution of 
1941.) Nov... 19 1 46s 18224 | 3203 0724 eres . = 
1943 | June 13 14 12| 142.7 | 40.9 | 7.0 | 7.4 epicentres among the areas. In the table, 
Sept 10 ie a nae oe . 7.4 the frequencies for (A+B) and (A+B+C) are 
1944] Dec. 7 18 36 | 136.2 | 33.7 | 7.4 also given. 
LO) 23 SOR 2420 ESS sO nie 


Table Il. Number of Earthquakes. 


a. Area A. 


: iM 
BaiVesalecelk s 
Year \ 


re 
kh 


Oo 
pelo hE 


OI 
oo 


HoH a 


DoH et bo po 
Peed) 

py 

ip 


wr ocUND | 
= 
ev) 


= 


PRNW HRW WO 
Hoo DoD He bo cob 


= 
lor) 


Py 
— 
OO NON WWHWe 
09 
—" 
ee =) 
re 
on 
ns 


we 

— 
rT 
ie.) 


NEw AHepw 
_ 
~] 


mee DO 
e bo 


[alll 5S cee 52 8 ell ell SS 


NRNRrw Oe RH WH REN Ne 
bo 09 0 
eb 
~~ 
— [eel SOM SM ct | 
oe 
i 
e 
= 
to 
oO 


SOON ATE WNROUOMON ATR WD 


WHNHHA wrote 


02 CO et DD et DO PO 
09 0 HY eB bo 


NwWNWNre 


it 
1 
3 


bo ce 


NNR 


ao 


1 1] 


l= 


en oe 


| 1 3] 7 18| 28 2 39 


35 43, 43 36) 41 36 27 21) 34/10/1831 8/101 3 


ey Magnitude-Frequency Relation for Harthquakes in and near J apan. 49 
y a ey b. Area B. 
95.0 a a c Ali 

2 2\ 1) 2 1 1 9 
2 1 ee i 4 a ADE 
3 1 Lied 1 Zi il rh 
4 a Li 3 
5 2 iti 3 
6 1) a 1 2 5 
‘ii ee a 3 it 1 8 
8 1 1 iz, 
9 | 1 i 
40 1 2 ob) al 5 
1 We a 3} 1 i 9 
2 1 1 A il 3 8 
S 1 at hi 4 
4 vl| eal Bua tient Piet il 9 
5 Ayah eal ae 1 it a yA 11 
6 il eles all ok a ie ae 
7 SU ad Cpe Wee Le etin to 19 
8 1 pL US ena or Bt) 1 A ay eG 
9 | LY 2) lly 3 i 1 i 
50 | ae 2| il 1 tft 1 6 

1| 1 2 3 4 6/10 9 


Go) AreaG: 
1 ) = i 
5.0, <1) 2.8) .4l 6 6 .7| 8 96.0 3 .2| .3) .4| .5] 6) 7] .8| .9/7.0| .1) .2| .3] .4| 5] .6|.7 Tota 
hear | | | | | | 

2 1 1 at 1 1 5 
3 1 et 2 5 
4 3 1 3 
5 1 il sei mal 1 q 
6 1 1 2 
7 1 leet 1 4 
i 8 1 1 4) 4) ah they 7 
eat 1 1 1 1 1 6 
40 al 1 a Al 1 5 
1 1 1| 1 1 1 6 
2 1 1 2 F 

3 2| 1) 2} 1 al al 2 i 
4 1} 1 1 3 
5 1; 1 slay, gah ib ’ 5 
6 1 1 2 
7 1] 2 Heh ky 6 
8 tite tp 1 1 6 
9 Ap 1 1| 2 6 
50 aed {ht ee. 

| eee eae ac Sat mle) alesh eat al cl 

“Total | 1| | 1 2 4| 4 3 5 512/11] a 7 ol 6 a 2 5 1] 3 2 2 2 | 104 


P ane eae 


Chuji TsuBot. 


50 


d. Area (A+B). 
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Fig. 1. The areas A, B and C. 


In Fig. 2, the logarithm of the mean annual 
number N of earthquakes is plotted against 
M. The value is at its maximum at about 
M=6.0 and decreases toward its both sides. 
While it is only natural that the number 
decreases toward larger M, that it does so 
also toward smaller M is due to the circum- 
stance that not all earthquakes whose magni- 
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Fig. 2a. Logarithm of the mean annual 
number of earthquakes N against . 


magnetude M (Area A). 


Fig. 2b. (Area B). 


Fig. 2c. (Area C). 
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Fig. 2d: (Area A+B). 


Fig. 2e. (Area A+B+C) 
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tudes are less than 6.0 are reported and there 
are therefore many of them which fall out 
of the statistics. The number of earthquakes 
having magnitudes greater than 6.0 only can 
be trusted upon. 

As to the relation between N and M, GuTEn- 
BERG and RicaTer proposed a formula of the 
form 

log N=a+b (8—M) 
in which log N is linear with respect to M. 
Here it is important to note that N is the 
mean annual number of earthquakes having 
the magnitudes between (M— 4M) and (M+ 
4M), where 4M is half the class interval of 
M for which JN is counted. 

The numerical values of the constants @ 
and } in the above formula were determined 
by the method of least squares so that the 
mean annual number of earthquakes in the 
areas A, B, C, (A+B), and (A+B+C) can be 
expressed, taking those earthquake which 
have greater magnitudes than 6.0. The 
results obtained corresponding to 4M=0.05 
are: 


Area a b 
A —1.60 1.06 
B —1.57 O22 
C —1.61 0.66 
(A+B) —1.46 1.04 
(A+B-+C) =e ae | 1.01 


It is interesting to note that while the con- 
stant b differs considerably from area to area, 
the constant a@ does not so very much. As 
is clear from the expression : 

log N=a+6b(8—M), 
the constant a@ is nothing but the logarithm 
of the mean annual number of earthquakes 
having the magnitude M=8.0. The approxi- 
mate constancy of the value of a for different 
areas therefore amounts to mean that the 
mean annual number of such large earth- 
quakes does not differ much from area to 
area, although there are considerable dif- 
ferences among the total numbers of earth- 
quakes. The total number of earthquakes is 
decidely the largest in the area A, and is 
followed by B and C as may be seen in Table 


lll. The constant b is also the largest in A. 
It was found that 5 is almost proportional 
to the total earthquakes number in the respec- 
tive areas as is shown in Fig. 3. This would 
indicate that if the total number of earth- 
quakes. in a certain area is twice as much as 
that in another, for instance, it does not 
mean that earthquakes of various magnitudes 
are equally double. Earthquakes with smaller 
magnitudes increase more rapidly than those 
with greater magnitudes do. 


Table III. 


Number of earthquakes of various 
magnitudes in the areas A, B, and 
C, (1931-1951). 


p eoagegty | B | Cc | Total 
6.0 43 15. “h) ak ieee 
6.1 36 6c | Se 
6.2 41 14: b= 9% 1 | 2 
6.3 36 11 9 56 
6.4 7 | 8 Gi tee Fae 
C55 Ve oy Sea Sy ieee 
ee f ae eae 2 | 46 
6272-1) 61, ea eer 14 
6.8 --| Is | peasy! 5 23 
UE ee ee 1 11 
7.0 10 3 2 15 
TA 3 0 2 5 
7.2 3 1 2 6 
7.3 1 0 3 4 
7.4 3 2 0 5 
7.5 2 0 0 2 
7.6 
7.7 1 1 0 2 
Total | 292 90 | 64 | 446 


On the other band, basing on their own 
materials, GUTENBERG and RicHTER gave for 
Japan 

log N= —0.90+0.80 (8—M), 


corresponding to 1/4 units of 4M. From the 
data we have been using, the following rela- 
tion is obtained : 


log N= —0.85+1.22 (8—M), 


corresponding also to 1/4 units of 4M. It is 
not surprising that while the values of a@ agree 
rather well, those of b differ sensibly from 
each other. The difference is presumably due 
to the circumstance that relatively smaller 
earthquakes with the magnitude a little over 
6.0 are apt to be partly dropped out of the 
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Fig. 3. Dependence of the constant b on 
the total number N of earthquakes. 


statistics of GuTENBERG and RicuTeErR based on . 


the distance earthquake observations. 

As was stated before, GurenBERG and RICHER 
proposed a formula of the form 

log N=a+b(8—M) 

for expressing N in terms of M. Strictly 
speaking, precautions must be paid before 
this expression is used, because it does not 
admit a linear superposition with respect to 
N. Suppose the constants a, 6 for a certain 
interval of time T, are a, and 0b;, and those 
for the same area but in another interval of 
time JT, are a, and b,. Then, 


log Ni=a,+ bi (8—M), 
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log No=a2-+b, (8—M), 
where a, and 6; may differ from a, and by 
respectively. If we take two statistics al- 


together, log (= T, Tp) comes into ques- 
2 


tion, but if the above two expressions are 
put into this one, the resultant expressions 
is no more linear with respect to M. In spite 
of this circumstance, it is always possible, 
by the very nature of the method of least 
squares, to determine the values of the con- 


si NiTi +N2T2 
stants a and 0 for expressing log ( TAT, ) 
in the form 
log (AN )= a+ 6-m™). 


The constants a and } may differ again from 
a, a, and b,, b:, respectively. With such 
variance, the constants have little significant 
meanings. da, a2, a, and bi, b:, 6 can be 
equal or at least approximately equal respec- 
tively with each other only when the time 
interval taken for the statistics is of such a 
sufficient length that the frequency distribu- 
tion of earthquakes with respect to M in 
that interval may be regarded as truely 
representing the average tendency peculiar to 
the area under consideration. If we are to 
attribute to the constants a, 6 any physically 


Table IV. 


The constants a and } determined for various time intervals. 


Interval a b Interval | a | b 


| 
AS 
| 


1931 


—0.40 0.52 

1931—1932 —0.97 0.79 
1931—1933 —0.67 0.60 
1931—1934 —0.80 0.65 
1931—1935 —0.86 0.69 
1931—1936 —1.05 0.79 
1931—1937 | —1.25 0.90 
1931—1938 —1.10 0.85 
1931—1939 | —1.15 0.86 
1931—1940 —1.23 0.90 
1931—1941 —1.25 0.91 
1931—1942 —1.30 0.95 
» 1931—1943 —1.35 0.98 
1931—1944 —1.36 0.99 
1931—1945 —1.44 1.04 
1931—1946 | —1.39 1.00 
1931—1947 | —1.41 1.01 
1931—1948 | —1.48 1.03 
1931—1949 | —1.46 | 1.04 
=—1.46 1.04 


1931—1950 | 


1950 = (ple) Oezil 
1950—1949 = 05083 0.36 
1950—1948 —0.90 0.64 
1950—1947 = 126 0.88 
1950—1946 Sesh | 0-91 
1950—1945 1h, | 0.90 
1950—1944 —1.30 | 0.91 
1950—1943 Sal cie | 0.96 
1950—1942 —1.46 | 02 
1950—1941 —1.49 | 1.08 
1950—1940 —1.56 1.07 
1950—1939 ae 56 1.07 
1950—1938 ile 1.08 
1950—1937 —1.46 | 1.08 
1950—1936 = hl | 1.04 
1950—1935 —1.48 1.05 
1950—1934 = 1,52 107 
1950—1933 =e 1.01 
1950-1932 ails | 1.08 
1950—1931 ~1.46 | 1.04 


Se 
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—@—e (starting from 1931) 


.--x---x (starting from 1950) 


= 10 20 


Number of years 
Fig. 4. Dependence of the constants a 
and 6 on the number of years 
taken in the statistics. 


significant meaning related to the nature of 
mechanism of earthquake occurrence, much 
care must be paid therefore as to the length 
of the time interval to be taken for the 
statistics. In this connection, it will be in- 
teresting to investigate how the constants a 
and 6 tend to their respective final values as 
we lengthen the time interval of statistics 
uccessively. Taking the area (A+B), the 


constants a@ and b were determined for re- 
presenting the mean annual earthquake 
frequencies in 1931, (1931-1932), (1931-1933), 
.....(1931-1950), adding one year successively 
each time. The same thing was done starting 
from 1950 and going back the other way, the 
interval being taken as 1950, (1950-1949), 
(1950-1948), .... (1950-1931). The results of 
calculations are given in Table IV and in 
Fig. 4. ; 

As is clear from Fig. 4, both a and 6 change 
conspicuously at first, but gradually tend to 
their respective asymptotic values when more 
and more years are added into the statistics. 
Beyond about »=15, the values of @ and b 
may be regarded as approximately stationary. 
How many years are necessary in order to 
get the stable value of a and b depends on 
the fluctuation in the mode of occurrence of 
earthquakes with particular reference to magni- 
tude. This presents an interesting problem 
and the writer expects to publish the results 


‘of his study of this problem in a near future. 
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Retarded Photographic Recording of Earthquake Motions. 


By 
Chuji Tsusor and Yoshibumi Tomopa. 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


An apparatus is described by which earthquake motions are recorded photographically 
with a certain time lag. The initial parts of the motions do not fail to be recorded 
and the speed of recording can be increased as much as desired. 


Earthquakes happen very rarely. It is a 
rather paradoxical fact that a seismograph is 
almost always indicating by its own record 
that there is nothing worth to be recorded. 
This is the very circumstance that makes the 
speedy running of the recording paper or film 
impracticable, if not impossible, for high costs 
and frequent attendances to the instrument 
indispensable for that purpose will be pro- 
hibitive. It is true that an ordinary starting 
device or a starter, as it is commonly called, 
which is appropriately connected to a seis- 
mograph drum and which makes it begin to 
rotate only when earthquake waves arrive 
the instrument is often useful for overcoming 
the difficulties above stated. But, since more 
or less sensible vibrations are needed for a 
starter to be excited to operate, initial small 
vibrations of earthquake motions can intrinsi- 
cally not be recorded and this is often fatal 
for many seismometrical purposes. 

An arrangement which was _ originally 
described by P. G. Ganz, H. J. Loair and J. H. 

_StTEPHEN (1949), which enables a_ retarded 
recording of earthquake motions is extremely 
interesting in this connection. The essential 
feature of this arrangement is that the mo- 
tions are temporarily ‘‘remembered’’ by some 
or other arrangement before they are ulti- 
mately recorded. In other words, the record- 
ing of the motions is made some time, say 7, 
after they actually take place. The initial 
portions of earthquake motions, however 
small, can in this way never fail to be 

-recorded, only if the final recording is started 
within time ¢ after the motions. For making 


the instrument temporarily remember the 
vibrations, GANE, Locim and SrepHen used a 
magnetic recording wire. 

The present writers have constructed an 
apparatus which works essentially on the 
same principle as that described by the above 
three authors, but which is purely mechano- 
optical and is therefore more dependable in 
routine observations. 

The working principle of the new apparatus 
may be understood in Fig. 1. Ais the last 
stage lever of a seismograph which makes 


Schematic sketch of the apparatus. 


Rigel. 


rotatory oscillations around the axis O at the 
time of earthquakes. Along the lever is 
attached a glass capillary tube G, the one 
end of which is dipped down into a small 
vessel I containing commercial Chinese black 
ink and the other rests lightly on a rotating 
cylinder D made of transparent glass. 
After flowing through the tube, ink comes 
out from its nozzle N and consequently a thin 
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black line is drawn on the cylinder. The 
cylinder is about 20cm in diameter and 10cm 
in breadth and is rotated at the rate of 1 
cm/sec. It is rotated continuously at the same 
position, no lateral shift being given as is the 
case in the ordinary seismograph drum. As 
the glass cylinder continues to rotate, the 
black line drawn on its surface is erased by 
the erasers E;, E;, and E;, which are all 
cloths wound cylindrically about respective 
axes and are rotating in the directions oppo- 
site to that of the drum. The first eraser Ei 
is fed water from a spray W. ‘The wet 
surface of the cylinder is wiped by E2 and 
then suitably heated by a small electric 
heater H before it meets the last eraser E;. 
In this way, a clean dry surface of the 
cylinder always comes up to the top to receive 
ink from the nozzle N. 

When earthquake vibrations arrive, the 
lamp L is lighted and the drum of a camera 
C which carries photographic film around it 
begins to rotate, by means of an ordinary 
starting mechanism S. By means of a lens, 
the image of the seismographic record 
drawn on the glass cylinder is focussed on 
the plane of a horizontal slit placed in front 
of the film. According to the rotation of the 
glass cylinder, the image moves across the 
slit, and since the film is also moving, the 
record drawn on the cylinder is photographed 


on the film. The film runs at a speed of 0.6 
cm/sec in our present apparatus, but the speed 
may be much increased if desired. As was 
stated before, the small initial portion of 
earthquake vibrations never fails to be 
recorded, if the starter comes into action 
before that very portion is in the position 
just to be photographed. This time of allow- 
ance corresponds to a quarter of the period 
of revolution of the glass drum and is 15 sec 
in our apparatus. The speedy running of the 
photographic film makes minute detail of the 
vibrations easily detectable. Five minutes 
after the beginning, everything stops auto- 
matically and comes back as it was at the 
initial state, so as to be ready to receive 
vibrations of a next earthquake. 

As may have been understood by the above 
descriptions, the merit of the present ap- 
paratus lies in the fact that there is prac- 
tically no loss of film and the speed of its 
running can be raised as much as desired. 
The apparatus has proved itself to be durable 
and dependable. Fig. 2 is a reproduction 
of a part of the record actually obtained at 
Tokyo on July 12, 1951. 
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Fig. 2, An example of the records obtained. Actual size. Numerals in sec. 
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On a Self-Exciting Process in Magneto-Hydrodynamics (ID). 


By 
Hitoshi Taxkrucni and Yasuo Saimazu. 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


The present paper is concerned with the problem of the self-exciting process in 
magneto-hydrodynamies studied in a previous paper (referred to as Paper I). In §2 
and § 3, by solving a differential equation numerically, the results in Paper I are checked 
in a conclusive way. The existence of a self-exciting dynamo having thus been shown, 
the detailed structure of this dynamo is studied in §4and §5. In §6, a new approximate 
method of solving the problem is proposed and applied to the problem in Paper I. This 
method of approximation is found to be fairly useful. By using this method, it is 
shown in §7 that the self-exciting dynamo is possible by the fluid motions of consider- 
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able varieties. The orders of magnitudes 
be the same as that in Paper I. 


§1. In a previous paper (which will be re- 
ferred to hereafter as Paper I (1952)), the 
present writers showed that a self-exciting 
dynamo is possible by which the earth’s main 
magnetic field may be produced and main- 
tained. The self-exciting dynamo considered 
was supposed to have the structure similar 
to that supposed by W.M. Etsasser (1946, 
1947) and E.C. Butuarp (1948, 1949) in their 
studies of the secular variations of the earth’s 
magnetic field. In the present paper, using 
the results obtained in Paper I, further studies 
will be made on this self-exciting dynamo. 

§2. In our self-exciting dynamo, the velocity 
field (relative to the earth’s mantle) of the 
conducting fluid in the earth’s core is assumed 
to be represented by arbitrary linear combi- 
nations of the following. two velocity fields. 

Ti°-type velocity field 


i= COS0 2 


1 r 
ViK7)=7V ul 1—-(Z) 
S.2°-type velocity field 
—6rV 2°°(7r) COs 2 Py = 
dV 9”° »¢\0(Cos 26 P»*) 
(2 Eso ieee amie 
O( ae 


\—( J ysin 006’ 


of the fluid velocities required are found to 


=3 sin? 6 


In (2.1), the (7, 6, ¢) components of the 
velocities refer to the polar axis which coin- 
cides with the rotation axis of the model 
earth, and a is the radius of the earth’s core. 
Vi,o and V2, are certain constants having 
the dimensions of velocity. Velocities given 

(2.1) become zero at r=a. As was point- 
ed out by BuLuarp (1949), the T1°-type velo- 
city in the earth’s core is considered to be 
much more predominant than the S,.°*-type 
velocity. So we shall, hereafter, confine our- 
selves to the study of the case in which 
V 1,0 


>. 2c 
In order to study this case, we must solve 


the following eigen-value problem. (see §5, 
Paper I). 


at. 0 


mov Ae Oa 


_2aT. 28 dT, 36x12 , 


de®. ee dé & de on ame 35.4 pl Se ag 
x[ga-pP +¢—up Tee |=0, @2) 
under the boundary conditions 
ive Re cea ies Se eS (2.3) 
& T° dT» pa 
and any=£ ae =. 198 Fig ee gee Taran 
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at €=1. In @.2) 


F= S and y=4zKaV 2,» , (2.5) 


« being the conductivity of the gy. core. 
In §5 in Paper I, the eigen-value Y = ay was 
found to be at about Y=100. In order to 
obtain the exact value for Y, the determinant 
A in (5.11) in Paper I is evaluated for several 
Y’s. The results obtained are shown in 
Table I and in Fig. 1. The exact eigen-value 
is thus shown to be at Y=130.6 (y*=2400). 


atc 
making adequate linear combinations of the 


A 


100 200 


Fig. 1. Y—J Curve. 

§3. In Table II (y7=2400) and Table III 
(y?=0 and 4890), are shown the solutions of 
(2.2) for which @ (&) in (2.4) becomes zero 
These solutions are obtained by 


solutions I and II (defined by (5.8) in Paper I) 
for each y*, and satisfy one of the two im- 
posed boundary conditions. 


Solution for y*=2400. 


Table I. Values of 4 for Y’s. 
yf VE A 
: 0 0 +18 
1200 Jaca +12 
1658 108.5 +8.2 
2400 130.6 —0.3 
3600 160.0 — 23 
4800 184.8 —63 
Table II. 
g T? TT! | 
0 ie | 1.281836 » | 
0.1 0.8706906 — 1.334939 
0.2 0.7167810 —1.870458 
(On8! 0.4687736 — 3.204420 
0.4 0.070183831 | — 4.677167 
0.5 —0.4210779 — 4.803987 
0.6 —0.8152594 — 2.742478 
Out: — 0.9294623 0.487666 
0.8 —0.7486155 | 2.986053 
0.9 —(0.3851389 3.919917 
1.0 0.00522213 3 


3.801665 | 


(0.00562) (8.805) 


Tor? Torr Torr 
0 0 0 
1.940291 — 49.48336 — 667.0607 
— 9.452420 — 88.86093 87.49544 
—16.11707 — 22.7512 1301.714 
—10.57754 139.9516 1663 .017 
9.721807 237.7424 — 8.738218 
29 .60671 125.7603 — 2018.81 
31.41339 — 84.74552 —1747.611 
17.13300 — 169.7589 75. 10816 
2.659509 — 105.6746 924.9715 
— 3.761189 — 30.81803 475.1291 
(— 3.876) (— 26.60) 


Table III. Solutions for y?=4800 and y?=0. 


ee 


é T° To’ 
0 1 — 15.2798 
0.1 — 0.528148 — 15.2677 
0.2 — 2.00586 ~13.5976 
0.3 — 2.97298 — 3.73610 
0.4 —2.33719 17.92768 
0.5 0.561066 37.5650 
0.6 4.29736 31.7832 
0.7 6.03172 0.54075 
0.8 4.45712 ~ 29.6135 
0.9 0.749519 ~ 41.4841 | 
1.0 ~3.410270 ~ 40.4950 
(—8.4101) (—40.417) 


Solution for y2= 0 is T.9=1 for all é. 


To" Tor Torn 
0 / 0 | 0 

1.70375 | 107.389 | 3893.37 
43.9570 828.173 | 9052.79 
162.5292 | 1351.721 | — 1999.88 

246 .823 | — 6.1808 | — 25773.3 

99.4471 | —2781.986 —10417.4 

— 216.182 — 2834.92 | 20351.0 

— 356.045 229.517 | 31746 .4 
— 215.848 2067.85 2618.63 
— 35.3183 1288.14 | ~12791.16 

39.0966 337.942 ~43761.8 

(36.230) (348.73) 
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As is seen from Table II, for the exact 
eigen-value 1?=2400, the other boundary con- 
dition (2.3) is satisfied simultaneously. For 
the solution given in Table III, this is not 
the case. T,%)’s (n=0, 1, 2, 3, 4) in Table 
II and III are connected with each other by 
the equation (2.2). A rough numerical check 
for the results in Table II and III is made 
as follows. Executing, for example, the 


integration | "7 T.°’’” (©) d& numerically, we 


0 


have 7.” (E) for each &. The values of 


T.”’’ (£) thus obtained should agree with 
those obtained by direct numerical integration 


of (2.2). "TB" (B) dE, , 


ee. 0 
T.” (&) d& (Stmpson’s method of numerical 
0 


The values of 


integration is used) are shown in the brackets 
in Tables II and III. These values agree 
with 72” (€=1), ...., 72° (E=1) satisfac- 
torily well. Small discrepancies among these 
values are probably due to the roughness of 
the Simpson’s method of numerical integra- 
tion. In short, in the last and the present 
sections, we have shown conclusively the 
existence of the real eigen-value ¥’=2400 (Y 
130.6) for the equation (2.2) with (2.3) and 
(2.4). A self-exciting dynamo is thus found 
by which the earth’s main magnetic field 
may be produced and maintained. Having 
thus got the exact eigen-value, we shall now 
study the detailed structure of the self-exciting 
‘dynamo. 

§4. In studying the structure of the self- 
exciting dynamo, uses will be made of the 
following expressions for the (7, 9, ¢) com- 
ponents of the velocity and the magnetic 

fields. 
T,°-type velocity field 
0, 
i 


V bogey %E) sin @, 


V 


Mya cos 8, 


curl V 1 
eV ana E ra 


0, 
Vi%E)= ISé: 


(E°V 1°) sin , 


S»2**-type ace field 
: igure qEV ME) cos 2¢ sin? 0 


tp 3 
i : 
—3V,2o(E* si RRA, 


6V 2,00( ” ) sin 2¢ sin 0, 
0 


V 


Beek ae Mast (pp oe tee *)sin 2¢sin8, 
V2, 26 
ats Ss ( 7 
Ve lee 

T,°-type magnetic field 
0 
tel 0 
—6&°T.(E)X sin 20, 


)cos2¢sin20, 


24% 5 Ty (78 cos 26-41), 


curl H dT, 


one (e% 


p+3Ty ) sin 20, 


Si°-type magnetic field 
—2S1°E) cos 6, 


H} (c& aS!  +25:°) sin 6, 


2a 


T.°5-type magnetic field 
0, 
Hi —6#T275(E) cos 2¢ sin 8, 
3E°T.”’ sin 2¢ sin 20, 
aa 187 € T(E) sin 2¢ sin? 0 


curl HJ —3- 1 or +3ET: a sin 2¢ sin 20, 
saat a )cos 26 sin 6, 
a 
T@)= yf cape +0—upT], 
HE OLW. yas (4.5) 
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3 dial iations of velocity fields. 
The ‘values Of ME VEO) ant “EEO Gee 


| pOrmOm eee 28) ate - 8) 
the equations (4.1) and (4.2)) are shown in  ——s _ Sng 5 gent 
= ate 0 | 
Table IV. The values of 6£?T,°, 2S.°, (€S:° 0.1 | 0.36 0.54675 
+2Si°), and 67.25 (see the equations (4.3), a4 | oe Oe ecare 
(4.4) and (4.5)) are shown in Table V and in 0.4 | 0.96 | 0.9720 
Figs. 2, 3 and 4, 2S,°(€=1) being taken ee | ee eat 
2S1° (E=1)=1 tentatively. 0.7 0.84 0.42525 
, ae ; 0.8 0.64 0.2160 
Table V. Radial variations of magnetic 0.9 0.36 0.06075 
fields. 1.0 0 0 
; 6eT. : 28,9 ES19/ +28, 6e°T 78 

Om eas 0 a " = 0.403753 “2298 0140875 r 0 ee 

Ook 0..0°914169 0.480666 0.47931 0.773998 

0.2 0 .07301030 | 0.688403 1.12906 1.93851 

0.3 0.07442964 1.263152 2.20530 1.47394 

0.4 0. 07117817 1.695332 1.88497 — 2.32076 

0.5 —0.0110526 | 1.257983 —1.19524 — 7.06826 

0.6 —0.0308149 — 0.068727 —4.52108 —5.46677 

0.7 —0.0478179 —1.307977 — 4.25821 5.69262 

0.8 —0.0499679 —1.660091 —1.32671 17.93345 

0.9 | —0.0328391 —1.360353 0.44662 18.40224 

U@ | 0 | — 1.000000 0.50000 0 

°2 
*10 
rigone 
0 ! fs 
4 8 
mars | 
Fig. 2. 6é277,0 0 1 
t Fig. 4. 6227.28 


The maximum value of 6£?T,° (and also of 
6&T.**) is at about £=0.8, and the ratio of 
6X &T,°(the maximum absolute value through- 
out the earth’s core of the T.°-type magnetic 
field) to 2S,° (E=1) (the maximum value at 


abe the surface of the earth’s core of the So 
ig. 3. - 2S)9, --- &(S,°)/+2S,9 type magnetic field) is about 


On a Self-Exciting Process in Magneto-Hydrodynamies (II). 61 


6X Eis" = 0.05x- 


58 (4.6) 


Taking 2S:° (€=1)=3 gauss and taking X 
= 1000 (see Table II in Paper I) tentatively, 
we have in the earth’s core, the T.°-type 
magnetic field of about 150 gauss. This 
result is worthy of notice. In his study on 
the secular variations in the earth’s magnetic 
field, BuLLARD (1948) concludes that many of 
the observed facts can be accounted for by 
an induction current caused by a circular 
motion with a diameter of a few hundred 
kilometers situated near the surface of the 
core, but in order to get numerical agree- 
ment, an improbably strong (inducing) field 
at the core has to be assumed. This diffi- 
culty is, however, removed by taking the 
above self-exciting dynamo into consideration. 
This dynamo, besides giving us an explana- 
tion for the maintenance of the dipole magne- 
tic field of the earth, provides us also a way 
to have a large magnetic field in the core. 
The secular variations of the earth’s magne- 
tic field may be caused by an eddy motion 
of rather local type in the presence of this 
large magnetic field. 

§5. In this section, we shall examine the 


case when X=4rra- Y= 1000. Since the 


value X=10* is taken rather tentatively, 
great importance should not be attached to 
the results in this section. Furthermore, the 
results may be changed drastically by the 
future hydrodynamical studies on the fluid 


motion in the earth’s core. 


Vi,0 


[1] In the expression for X, = is the 


maximum value of the 7T1°-type velocity 
field. Taking a=3.5x10%cm. and «=3x10-° 


Vi,0 
e.m.u. aS are used usually, we have <e 


=(0.076 ae The value of Y for this case 


is obtained as follows. As is seen from Table 
Il in Paper I, the approximate value of Y 
for this case is nearly the same as X=©o. 
In short, X=10? is taken to be the same as 
X= co from the physical point of view. Thus 


it may be allowed to take Y =130 (see Table 
I) for X=1000. The corresponding maximum 


value of the S.?*-type velocity Lg be- 


© 
8 
So.  teaeen secs 
[2] The energy dissipation W by viscosity 
(the kinetic energy transformed into heat by 
friction) per unit time is 


comes Vaoo=Ot th pe —— 


w= al curl V/dv, (5.1) 
where y: is the coefficient of viscosity. The 
integration in (5.1) extends through the earth’s 
core. Since we are studying the case in 
which X->co, it is sufficient for our purpose 
to take the T;°-type velocity only in the V 
in (5.1). The value of W: thus calculated is 


Wi = PndHy,,.) 25.5 x x10 oe? 
SEC; 


(5.2) 
in which a and Vi, are taken to be the same 
as in [1]. 

[3] The energy W, dissipated as the JouLE’s 
heat per unit time is 


ire HH) dv. (5.3) 
Taking only the T."°-type magnetic field into 
consideration, we have 


(5.4) 


W2=0.075xa(aV,0)"=8.8% 10" 0 


The energy dissipation takes its maximum 
value at €=1, i.e., near the surface of the 
earth’s core. Since the mass of thc core is 
1.910" gr., the corresponding JouLn’s dis- 
sipation in each gram of the core is 1.1 10-" 
cal. 

sec.” 

[4] The ratio of the energy dissipation by 
the mechanical friction to that by the JouLn’s 
heat is given by (5.2) and (5.4). This will 
become equal to 1 for the viscosity “ of the 
order of 10” c.g.s.. If the viscosity of the 
earth’s core is smaller than this critical value, 
the energy dissipation by viscosity becomes 
negligible compared with that by the JouLE’s 


heat. 
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[5] The energy dissipation W, of the field 
is compensated by the work done on the 
field by the fluid motion. This will be un- 
derstood as follows. The pondermotive force 
F of the field upon the current is given by 
the well-known FLEmMine’s law 


F=IxH=;- curl Hx H. (5.5) 
Taking the negative of the work done on the 
fluid by this force per unit time, we have 


Ws=—\V-F d= j,\V- xcurl H)dv, 


(5.6) 


the integral extending through the earth’s 
core. This is nothing but the work done on 
the field by the fluid motion. Using the 
results obtained in the last section, we have 


(5.7) 


W528 8x 1002 =e 


‘ = 
SEC. 


Thus it is shown that the energy flow in the 
course of hydromagnetic interchange is from 
the fluid motion to the field. The work W; 
takes its maximum value at about €= s In 
short, the field takes its energy in the deeper 
part of the earth’s core and consumes it near 
the surface of the core. 

[6] If the fluid motion and the field are to 
be stationary, there must be a source of 
energy which compensates the energy dis- 
sipation by the mechanical friction and the 
Joutr’s heat. The energy transformation in 
our self-exciting dynamo may thus be shown 
schematically asin Fig. 5. Taking the results 
in [3] and [4] into consideration, and taking 
u=10" c.g.s. tentatively, the required energy 
is estimated to be of the order of 2x10-" 

cal. 
gr. sec.” 
Panertu (1942) were found to contain radio- 
active substances, and to generate Abs lig) 3 

cal. 
gr. sec.” 
large an amount as that required above. 

[7] The current density J at the surface 
of the earth’s core is calculated by the rela- 


Six meteorites examined by F. A. 


This is about twenty times as 


_ ENERGY SOURCE 


FLUID MOTION 


Wa 


Fig. 5. Schematic representation of the energy 
transformation in the present model of the self- 
exciting dynamo. 


tion J= pcurl H. Taking only the 7.°-type 


magnetic field into account, we have the 
maximum value of J of the order of 2.7 10-7 
e.m.u. Two-dimensional current sheet of this 
current density and the skin depth of about 
10 km. will give values for the current den- 
sities near the surface of the core, which are 
appropriate to explain the secular variations 
in the earth’s magnetic field. 


§6. In the last section, we have examined 
the case when X=1000. No absurdity was 
met with in the course of this examination. 
In. order to complete the study of the self- 
exciting dynamo, however, there remains a 
difficult problem to solve. That is the pro- 
blem on the fluid motion. Is it possible that 
such a fluid motion takes place in the earth’s 
core? Is the fluid motion appropriate to 
maintain the self-exciting dynamo supposed 
in the present paper ? 

Putting aside these problems for a while, 
we shall now show that the self-exciting 
dynamo is possible by the fluid motions of 
considerable varieties. In (2.1), we put 
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eo Mat-(O0 


va Yfs-(2)) 


tentatively, for the radial parts of the velo- 
city fields. There is, however, no reason why 
the radial functions should be of these types. 
In order to make the matter clearer, studies 
will be made whether or not the self-exciting 
dynamo is possible by the fluid motions of 
the following types 


Vir(r)= yee), Va(E)=EI-EF, 


(6.1) 


V(r) ="ViE), Vi) =E"1-£), 


F 


= oa (6.2) 


where o and r are certain integers. By the 

same reason as is stated in §2, we shall confine 

ourselves to the study of the case in which 

Ee Thus, 
2,2¢ 


problem is reduced to solve 


as 11-$5.in Paper J, the 


aaceye. oh Silty 
TINE)= TEs HE) 
d 4 207 ad 20\ET. 
x| geEV Te) + SAEs ET], 63) 
@S 4 dSi° we 2 
a EVENT, 
Eee ge 5 Pe 
Of eae (6.4) 
ET.) 6 dT. 2 dV 
dee de ~3E dp ie 
under the boundary conditions 
ie” (6.6) 
dSv 
at f=: 


- It is, however, almost hopeless to solve the 
eigen-value problem for y by the trial and 
error method. We shall solve the problem 
by am approximate way which is similar to 
that adopted: in §4 of Paper I. - At first, 
taking (6.6) into consideration, we shall put 


T(E)=Ao(1—£)+ ArE(1—€)+ «.- 
eA nek ee) ets (6.8) 


fom T(E Ans 7=02 tee... ) in (6.8) are 
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certain constants to be determined. The radial 
function 


TE, n)=E"(1—&), (6.9) 
under A, will be called the »-th coordinate 
function. Inserting (6.9) in (6.3) and (6.4) and 
taking (6.7) into consideration, we get 


S.(6)=" 7p SD AnSiME, 7), 


(6.10) 

in which S,° (&, ) is 

1 (WT An +33 Gn +57 

S\N 440 N45 aN. 

4! +43. W- #12) Wan sae 
N48) t(NaaNGTE 


ORS Sil 
dn’ +33 ong 16m +57 5 
rT +FON GO 


SE, 2)=— 


~ (N+5)\N +8) 
4n’ +43 (EN 41 

~(N +7)(N +10) 
n ee lZ a +8) 

(N +8)(N +11) 
In (6.11), 2’ and N are 
=n+2o, N=n'—t=n+20—-1, (6.12) 
respectively. The constants A,’s will be 
determined by the following orthogonality 

conditions (see §3 and §4 of Paper I) 


+ 


(6.11) 


@T? , 6dT? 2 avs 
\er 0(E, mG ie lag 
air See ), (6.13) 


in which T.° and S,° are to be taken as those 
in (6.8) and (6.10). If we take (6.8) into con- 
sideration, we see at once that (6.13) is noth- 
ing but the equation from which (6.5) itself 
is obtained. Executing the integral in (6.13), 
we get the equations of the same types as 
those in Table I of Paper I. Equating the 
determinant formed by the coefficients of 
A,’s to zero, we have the equation by which 
the eigen-value y=4rKaV2,.. is determined. 
By the present method, better approximation 
for y is obtained with less numerical compu- 
tations than by that in Paper I. In fact, 
putting ¢=r=0 in (6.2) and using the method 
in the present paper, we have the value of 
y=sy as is shown in Table VI. In Table 


VI, the second approximation, for example, 
means the value of Y which is obtained by 
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using two coordinates functions and solving 
the matrix equation of rank 2. From this 
table, we see that our approximate values 
for Y agree with exact eigen-value satisfac- 
torily. Looking at the result from a different 
angle, we may say that the result in §2 and 
§3 of the present paper is checked in an in- 
dependent way. 


Table VI. The eigen-value for the case 
o=t=0 obtained by the new approxi- 
mate method. 


approximation Ve 
al tials Nea 
2 129.3 
3 130.2 
exact 130.6 


§7. Returning to our main subject, we shall 
now study the case in which o=r(=<0). For 
the larger value of o=r, the depth of the 
fluid motion becomes the shallower. In fact, 
the value of & for which the fluid velocity 
becomes maximum is given by 


E ==((250 say) i= = pe (7.1) 


The maximum value of the radial velocity 
of S.°*-type velocity field is given by 


Max. vel.=18&°*(1—Em)*V2,26, (7.2) 


where V2,2- is the constant in (6.2). Defining 
Y by the following equation 


Y =4r«a-(Max. vel.), (7.3) 


the value of Y is calculated for several og=r 
up to the second approximation. The first 
approximation for Y differs but little from 
the second approximation which is shown in 
Table VII. From Table VII, we see that the 
self-exciting dynamo is possible for all o=r 
studied, and the values of Y obtained do not 
differ substantially from that for g=r=0 in 
Paper I. In addition to the results in Table 
VII, the values of Y are calculated for several 
values of o and t (<a). Only the first ap- 
proximation for Y is calculated, and the 
results of calculations are shown in Table 
VIII. 


Table VII Table VIII 


The eigen-values for various types of 
fluid motions. 


o=t te o t Me 
} 
——~ =< c aa } 
0 129.3 2 0 56.2 
2 98.2 4 0 42.2 
4 44.5 8 Tle! 3200 
8 28.9 0 eet 154.0 
16 25.6 0 Ser 90.8 


In view of the results obtained in this sec- 
tion, we may say that the self-exciting dynamo 
is possible by the fluid motions of consider- 
able varieties. The only condition required 
for the existence of the self-exciting dynamo 
is that the radial velocity is of the order of 
magnitude for which Y in (7.3) becomes 


Y ~10°. (7.4) 
Putting z=3.5x108cm. and «=3x10-* e.m.u., 


the corresponding maximum radial velocity 
is calculated to be 


Max. vel. =0.01 soe (7.5) 


Is it possible that the fluid motion, of this 
order of magnitude, takes place in the earth’s 
core? Under what condition is the fluid 
motion possible? It, thus, becomes desirable 
to study these problems rather in detail. This 
will be done in a forthcoming paper. 
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On Self-Exciting Processes 
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In a recent important paper, Takrucni and 
Sximazu (1952) have investigated the possibility 
of a simple type of motion in a fluid conduct- 
ing sphere acting as a dynamo. At the end 
of their paper they thank the present authors 
for discussion and add, ‘‘It is, however, fair 
to remark here that Professor BuLLAaRp and 
Mr. GeLuMAN are of the opinion that the 
self-exciting process is impossible by the 
process considered here.’’ In fact, there is 
no serious disagreement between our results 
and those of Takrucni and Sarmmazu and the 
purpose of this note is to explain what we 
regard as proved and what as_ uncertain. 
We use their notation and refer to their 
paper TS. 

We have used equations exactly similar to 
their (3.9), but have worked mainly on the 
case where only an S.** motion is present (i.e. 
a=0 in TS). We have taken V,2¢ (7) as 


2 V.2¢ r“(1—r/a)? T(r) 


where z is a positive integer or zero and f(r) 
“has no zeros in0<r<a. This is more general 
than TS (3.10) which corresponds to u=0, 
f(n=1. With a=0, TS figure 1 and (3.9) 
fall into two pairs; (A) equations for S;° and 
T.*s and (B) equations for T,° and T.**. We 
have proved by the WBK method that (B) 
posesses two types of solutions for all w>0. 
In one type T," and T.”* have the same sign 
at rv=0; in the other they have opposite 
signs. We have investigated the first type 
numerically for f(r)=1, w=1 and w=14. The 
second type probably has a higher eigen 
value than the first, but has not yet been 
studied in detail. For w=0, f=1, the second 
type of solution exists but the first does not. 
We have proved by the same methods that 
equations (A) possess solutions for w>0. We 


in Magneto-Hydrodynamics. 


By 


and H. GeLLMAn. 


ry and University of Toronto. 


have not yet completed the study of these 
equations for “=0. 

We therefore agree with TS that solutions 
of TS (3.9) exist and have analytical proofs 
of existence not depending on the convergence 
of iterative numerical processes. We have 
done a considerable amount of numerical 
work by iterative methods and have obtained 
results converging about as well as TS Table 
II. We find it hard to judge whether such 
processes are converging and have one ex- 
ample where a solution was known to exist, 
but where convergence of the numerical pro- 
cess was still in doubt after 20 iterations. We 
suspect that the first line of TS Table I does 
not really converge. 

However this may be, there is agreement 
that there are systems for which TS (3.9) 
possesses real eigen values. The existence of 
such eigen values does not establish the ex- 
istence of dynamo solutions of Maxwell’s 
equations, since in TS figure 1 all harmonics 
beyond the second have been omitted. We 
have derived equations containing all har- 
monics and have calculated the coefficients 
for harmonics up to the fourth. Some work 
has been done on the solution of these equa- 
tions, and there is no reason to suppose that 
solutions do not exist. The existence of solu- 
tions does not however follow necessarily from 
anything that TS or we have done; though 
naturally the existence of solutions of the 
simpler equations makes it more likely that 
the complete set also has solutions. 

15th July, 1952. 
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Physical States of the Earth’s Core. 


By 
Haruo Mixt. 
Geophysical Institute, Faculty of Science, Kyoto University, Kyoto. 


Abstract 


1. The question, ‘Is the earth’s core gas or liquid?’? has been discussed from the 
standpoint of statistical mechanics. It has been found that for each element there is 
a critical temperature, below which the liquid state is stable in the eore and above 
which the gaseous state is stable and this temperature depends on the atomic number 
of the element. This temperature is named the critical temperature for that element. 

2. The molecular viscosity coefficient of the core has been computed. In case the 
liquid core is stable, the coefficient is 10-?~10-1 poise, while it is 10-? poise in case 
the gaseous core is stable. Moreover, it has been shown that we can estimate the 
temperature of the liquid core, if its composition is known, or, conversely, we can 
estimate the composition, if its temperature is known. 

3. Using the theory of visco-elasticity, the attenuation coefficient of the core for 
seismic waves has been obtained. The result agrees with observations of seismic waves. 

4. The electric and thermal conductivities have been obtained, using formulas of 
Following values have been obtained: 
2~4~x10% ohm-! em.-!, 
0.05~0.5 cal. em.-! sec.-! deg.-1, 


statistical mechanics. 
electric conductivity: 
thermal conductivity: 


67 


§1. Gas or Liquid ? 

In order to investigate the problem, <‘Is 
the earth’s core gas or liquid ?’’, we must 
first examine its Gipps’s energy, 


G=E+tPYV— IS, 


-for each case, where 


E: Energy 
iP Pressure 
V: Volume 
T: Temperature 
S: Entropy. 
The phase that gives smaller G is the one 
that is stable. In this connection, the follow- 
ing experimental and theoretical results in 
high pressure physics must be remembered. 
a) Theoretically, the density of hydrogen 
increases rapidly from 0.35 to 0.77 at 0.8 x 10° 
atm., making a transition to the metallic 
phase. Experimentally, transitions at high 
pressures have been found for arsenic, tin 
and phosphorous. Taking these facts into 
considerations, W. H. Ramsey (1949) attributed 


the rapid increases of density at the boundaries 
of the earth’s core and of its inner core to 
the first and second phase transitions, respec- 
tively. 

b) Although the earth is one of the planets 
belonging to the solar system, it does not 
differ very much in physical properties from 
fixed stars in that it has a great mass. Much 
suggestions can therefore be obtained, if we 
consider that the earth is now at one parti- 
cular stage of revolution of a star. The ab- 
normally great densities of white dwarfs 
were estimated by W.S. Apams in 1925. 
R.H. Fowter explained these great densities 
by assuming that all the electrons surround- 
ing atomic nuclei in such stars are split out 
of them owing to their so small mutual dis- 
stances and that the electrons form a dege- 
nerated electron gas. In studying the earth’s 
core, Ramsey assumed that only a small 
fraction of atomic nuclei is split of electrons. 
Thus the state in white dwarfs may be con- 
sidered as the extreme of what was assumed 
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by Ramsey for the earth. 

Let us consider that the earth’s core con- 
sists of free eleetrons and ions ionized by 
high pressure and that its structure resembles 
to the metallic bond structure if it is liquid 
or solid. Following O. K. Rice (1940), let us 
put the potential energy as follows: 


U=U44+U ez, 


Us= SAAN: \ 


V 
3h? 3. \?!3 
Ue= (0m) Nave), 
where A: A constant such that (Ai/7?)=2.0 
~2.4, 
ry: The number of free electron per 
one atom, 
h: PLANCK’s constant, 
me: Mass of electron, 
: The number of free electron per 
one mol. 
Then, 
EF s5=U4+UR+3RT, 


E,=UstU n+ SRT+RT ASE 


Ey =E electron + SRT, 


where E ¢ectron: Energy of free electron 
Vy : Free volume of liquid. 
In case of ErnstEIn’s model 
iets \a 


2nmiv? 


In case of DeEByr-Born’s model 


: -( fr rel. fx 3/2 
Jaks er ea aa, : 


Therefore, in either case, 


Alogv; _ 3 
IS eye 


Suffixes s, / and g represent solid, 


liquid and gas respectively. 
So, we get Bs=Ey=U4t+Uer+3RT. 


Furthermore, 


Opn 4 
TS;=TS7=3RT ) Slog 
. roe +s} 


; 5 2nmikT \3!2 V 
TS,=>-RT+RTI A 
eee Bes on|(“ he ) ne 


Drsye’s characteristic tempera- 
ture, 

k: BoLTzMANN’s constant, 

m;: Mass of ion. 
Neglecting PV and TS due to electrons be- 
cause these are common to the three phases, 
we have 

Gs=G =U 4+Ur+3RT—3RT 


x {log 4 3 ai 


G gE electron a NSB 


ciog{ 2A") 
= h Nz). 


Since G; and G, are equal, we cannot say 
which of the solid or liquid phases is more 
stable. We shall compare the stability of 
the gaseous phase with that of the liquid 
phase. 

Now, following Ricr, if we assume E prectyon 
=Ur, we get 

Gi.—G,=U4+RT 


ce 60 2 
i 


and if G:—G,>0, the gaseous phase is stable 
and if <0, the liquid phase is stable. 
As was shown before (Mrx1, 1952) we have 


h 


where Q@Qp: 


(1) 


Ov=7¥m, 
‘ Seas 4 i x 
4 ev 1 at 
De DS 
_ _ Olog ym 
Te ~“G log V’ 


where y¢ is GRUNEISEN constant. According 
to these relations, we get the following 
values of @p and vy» at 2900 km. by using 
Jerrreys’ density distribution. 


Atomic Number | Op | Yon 
10 | 1290°K 2.700 x 1013 
20 | 1033°K 2.152 x 1013 


30 877°K |; 1.825. 1038 


Also, the values of @p and »» at both sides 
across the core boundary are related as fol- 
lows: 
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Od, = = (22) Ob, ’ 
Oc 


Ye 
Ymse = 4% Yinys >» 


in which suffixes c and s indicate the core 
and mantle sides, respectively. The assump- 
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tion introduced by Ramsry that the composing 
element does not changed across the core 
boundary was adopted also in the above 
equations. 

Since p;=6.02 and p,=8.92, we get the 
following numerical values : 


Atomic | Ye=1.0 Ye=1.5 
Number ae ' ay 7 } 
Op, c Vm, ec Op, c Vimy e 
| Sete 
10 | 1910°K 3.98 x 1018 2320°K 4.84 x 1018 
20 | 1530°K 3.19 x10 1860°K 3.87 x 1038 
30 1300°K 2.70 x 1018 1580°K 3.29 x 1018 


Using these values, we can evaluate the second term in equation (1). 


The result is shown 


in Table I. 
Table I. The value of RT log|(“ arm ae On’ fat 
( x10” ¢.g.s. unit) 
2 Atomic ; 

Pe taichor 10 | 20 30 

TOK = Fs 7 a ie _ 4 bate a <— a fa 
= Ye=1.0 %e=1-5 | te= 1.0 fered We) Neue | Ye=15 
3000 1.98 2.12 2.25 2.40 2.42 2.58 
5000 2.97 3.22 3.41 3.66 3.72 3.96 
8000 4.28 4.65 4.96 5.38 5.46 5.85 
11000 | 5.44 5.98 6.40 6.94 | 7.08 7.61 
15000 6.88 | 7.58 8.18 8.90 9.10 | 9.86 
20000 8.40 9.40 10.2 | a Oye | 11.4 | 12.4 


- By comparing the above values with attrac- 
tive potential energy of liquid, we can examine 
the stabilities of the liquid or gaseous core. 
Before doing this, we must know the number 
of free electrons per one atom. ‘This can be 
done as follows. 

According to GRUNEISEN’s equation of state, 


_dU , 3RroT os 

dV Vale 
Furthermore, since we have assumed that 
the core consists of electrons and ions, we get 


CLL hee NH RE celia fe hepa 
dV 3 ( us 2071¢ i Vv. 

We may evaluate the number of free 
electrons per one atom assuming the relation 


between V and V, properly. Provided the 


P= U=U4tU e. 


volume per one atom is represented by 
-=e(2a)? f 
a: Radius of ion, 
the values of c corresponding to several crystal 
structures are as follows; 
The body-centered cubic lattice : 


diet dd 
(V3) 
The face-centered cubic lattice : 
1 ® 
(Ci / Re 
; eae ren oe ib, 
The diamond lattice: c= (3). 
The close-packed hexagonal lattice : 
1 
C= io ik 
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Let us take ae because the close- 


packed lattice structure is the most probable 
under high pressures in the core. The space 
for free electrons per one atom is therefore 
An 
3 
if an atom is represented by a sphere. 

If we denote the volume of space for free 
electrons per one mol by Vz, then 


a@=1,47a*, 


MIKI. 


V 
Vv 2230/5), 

Since P =1.39 x 10° atm. at the core bound- 
ary, we can find the relation between tem- 
perature and the number of free electrons 
per one atom from equation (2) and also the 
relation between V and Vz. 

The results of calculations are shown in 
Table II. 


Table II. Relation between temperature and the number of free 
electrons per one atom. 
Pe Rumer 10 ee | 80 

To, os Tee LO Ye=1.5 Y¢=1.0 Ye=1.5 | Y¢=1.0 geet: 
3000 0.1576 0.1420 | 0.351 0.335 0.612 | 0.600 

5000 0.1364 | 0.1063 0.333 0.305 0.590 | 0.561 

8000 0.1000 0.0300 0.298 0.255 0.560 | 0.515 

11000 0.0529 | — 0.266 0.195 0.527 0.465 

15000 — | — 0.218 0.094 0.480 | 0.385 

20000 — | — | 0.142 | — 0.422 0.281 


We can find the relation between tempera- 
ture and U4 from Table II and furthermore, 
examine the stability of gaseous or liquid core 


by equation (1). The attractive potential 
energy is given in Table III and Figure 1 as 
a function of temperature. 


Table III. Relation between temperature and attractive potential energy, —U4. 
(102 ¢.¢-s- unit) 
Atomic ‘ 
Sa 10 20 30 
SRG , ea ; Bis * ha : 2 4 or 
Ya¢=1.0 Yg=1.56 | Yg=1.0 ¥e=L5 Yg=1.0 Vaso 
3000 0.445~0.543 | 0.360~0.432 | 1.76~2.10 | 1.60~1.92 | 4.60~5.52 | 4.42~5.30 
5000 0.334~0.400 | 0.2083~0.243 | 1.58~1.90 1.39160." 228 —5 12 | 3.87~4.65 
8000 0.179~0.215 | 0.016~0.021 | 1.27~1.52 | 0.98~1.12 | 3.85~4.61 3.25~3.90 
11000 0.050~0.060| — ““— | T.01~t.21 | 0ls4~0l65 |) 341-210) Bene 18 
| | | 

15000 — —-) = — | 0.68~0.81 | 0.18~0.15 2.838~3.40 1.82~2.19 
20000 ~ — — | 0.29~0.85 | — — 2.19~2.62 | 0.97~1.16 


In the figure, we can see that each element 
has its own critical temperature corresponding 
to the atomic number that determines which 
of the gaseous or liquid phases can be stable 
in the earth’s core. Above this temperature, 
the gaseous phase is stable and below it, the 
liquid phase is stable. For instance, if the 
atomic number of the composing element of 
the earth’s core is 30, the critical temperature 


is 5000~7000°K and is below 2000°K, if the 
atomic number does not exceed 20. There- 
fore, if the earth’s core is of the metallic 
phase of SiO, (mean atomic number 10), as 
suggested by Ramsry, the earth’s core must 
be gas at some thousand degrees of tempera- 
ture as usually adopted. This result may 
answer the question why the idea of phase 
transition at the depth of 2900 km. can explain 
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the steep increase in density, but it cannot 
explain why the solid phase-changes suddenly 
to the non-solid one. 

The critical temperature decreases with the 
decreasing atomic number as we see in the 
above computations. 


§2. Viscosity 

No decisive conclusion was obtained in the 
preceding chapter, whether the earth’s core 
is in the gaseous or liquid phase. Gas and 
liquid are very much different in their struc- 
tures. Therefore, we shall divide our dis- 
cussions into two sections, corresponding to 
the gaseous core and liquid core, in the fol- 
lowing study of viscosity of the core. 

a) Viscosity of the liquid core 

In the theory of visco-elasticity, MAxWELL’s 
fundamental equation is as follows: 


aS= Gd dt. 


T 


where S: Stress 
¢: Strain 
G: Rigidity 


t: Relaxation time. 
Since dS=0 under a constant surface force, 


s=6r%., 

On the other-hand, 
_ ae 
P az ? 


z: Direction normal to the flow. 


Since the velocity of flow is dz/df=v in the 
definition of viscosity coefficient S=vxdv/dz, 
we see 


7 tae (33) 


Recent experimental results with high fre- 
quency technics have shown that the rigidity 
of liquid is about one or two figures smaller 
than that of solid. We can evaluate the 
value of viscosity coefficient if we know re- 
laxation time. 

H. Eyrine (1941) has shown that the visco- 
sity coefficient of liquid can be represented 
as follows, 


pa One ch Te lP?, (4) 


where 


2791, 1 
a ; ty se-U/ RP | 
/ ABA 


G’ is the work necessary to make holes of 
one mol in. the liquid, that is, G’=G,—G, 
and we evaluated this quantity in §1. Then, 
we can evaluate G, ct and y respectively from 
equations (3) and (4). 
If we adopt the Drsyr-Born’s model, 
1 


res e-U|RP 
oD ¢ 


and we can get the minimum relaxation time, 
Tm, aS U<0. These values are shown in 
Table IV. 


Table IV. Relaxation time for 
atomic number 30. 


[°K T2100 eet is 
1.40 x 10-13-~ 1.20 x10=18~ 
3000 Bilan 10228 4.02 x 10-22 
7.95 x10-M~ 5.65 x 10a 
oe 6.30 x 10-13 3.59 x 10-18 
ie 5.17 10-14 4.25 %10-14 


We can evaluate the minimum rigidity, 
Gn, as G’>0. These values of G are shown 
in Tabl V. 


Table V. Rigidity for atomic number 30. 
Gm: Minimum rigidity 


T°K jal) Gaasleo Gin 

4.87x10U~ | 4.26x100~ 

3000 | "7.02x104 | — 6.08x 100 

A268 100 4.49 x 101 
6.00 x10 | TIER 


2.05 x 1012 


5000 


3.41 x10" 


We can evaluate the viscosity coefficient 
from the equation (3). These values of 7 
are shown in Table VI. 


Table VI. Viscosity coefficient of the 
liquid core for atomic number 30. 


ILENE Ye=1-0 Ye— 0 


3000 6.81 x 10-2~3.64 5.11 x10-2~2.45 


3.41 x 10-2~ eke eh eo 
5000 Ginga tol 1.61 x 10 


i EEE EEE EEIE EERE 
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Minimum viscosity coefficient of the liquid 
core for atomic number 30. 


T°K 1@=1.0 vant 
3000 1.06 x 10-2 8.70 x 10-3 
5000 1.76 x 10-2 1.45 x 10-2 


It is natural that the rigidity increases with 
decreasing temperature as we see in Table V. 
The values of rigidity and viscosity seem to 
be near minimum values, because the experi- 
mental study shows that the value of G’ is about 
a~ ts of the theoretical. 
the bottom of the earth’s mantle is known 
to be 3.21x10". We have seen that the coef- 
ficient of rigidity of liquid is one or two figures 
smaller than that of the corresponding solid. 
Therefore the rigidity given in Table V is of 
reasonable magnitude and the temperatures 
3000~5000° are likely to be a good approxi- 
mation. From Fig. 1, these temperatures are 
seen not far from the critical temperature. 


The rigidity at 


x10"Ces. unit 


6 
4 
2 
Ss 5 8 " 15x10 PALS 
Fig. 1. Stability condition for the element of 


atomie number 30. 


The above result that the temperature of 
the earth’s core is in the neighbourhood of 
the critical temperature was obtained for 
atomic number 30. This can be said also 
for all elements of which the liquid core can 
be considered to be composed. Therefore, 
the viscosity of the liquid core is about 10-2 
~10-! poise. 

That the temperature of the liquid core is 
in the neighbourhood of critical temperature 


is an important result. From this, we can 
estimate the temperature of the liquid core, 
if the composing material is known, or, con- 
versely, we can know the composition if the 
temperature is known. 

b) Viscosity of the gaseous core 

The viscosity coefficient of gas is given by 


0 


Prey 
= 3 ym, 


by statistical mechanics, where 
N: Number of atoms in volume V 
M: Mass 
7: Mean free path 
v: Mean velocity. 

We have assumed that the earth’s core is 
composed of electrons and ions in §1. We 
shall evaluate the viscosity coefficient of mixed 
gas of electrons and ions. 

The mean free paths for the mixed gas are 


C— —— Us oA . 
4raZ2Ne+r(ae+ai)?Ni ‘ 

Se eee 

al 7(de+ai)?Ne+4zaz2N; " 


1 


where a: Radius of electron or ion. 
Suffixes e and z represent electron and ion, 
respectively. Since a.~0, the above equations 
are transformed to 


i Peal 
alas ke aay gs 
=> rai? yy + sae 
Consequently, 
1 Nim; 1 N.mev 
=71+7-=— ——_— = oe 
at 3 rairzNe+4ra;27Mi 3 za?N; 


On the other hand, we have for the pressure 


TES NEL INGNES ONG 
P= 25) ie seis 
el 7 ) ( Ve ) i V La 


The first and second terms on the right hand 
side are due to electrons and ions, respectively. 
We assume that ions are pluged closely in 
the gaseous core as in the liquid core, so 
that, V/V.=3.85. Then the number of free 
electrons per one atom, r=N,/N;, can be 
easily obtained taking P=1.39x10® Atm. The 
calculated numbers are given in Table VII. 
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Table VII. Number of free electrons per 


one atom in the gaseous core. 
AS TSE Cole Ss LS 


Atomic ] 
Number | 
ae i 10 20 .| 30 
: res EESEE eee aw 
3000 0.170 | 0.340 | — 
5000 0.164 | 0.335 — 
8000 0.155 | 0.325 0.550 
11000 0.145 0.318 0.546 
15000 0.132 0.305 0.535 
20000 0.114 0.290 0.520 


The mean velocities of electrons and ions 
can be determined from the following equa- 


tions : 
E me MeV” —— 38h? € N. 2/3 
2 407m, “48 \ee ; 
ine oy 
vg ive? 3 pr 


from the assumption. We can evaluate the 
viscosity coefficient of the gaseous core using 
the above numerical values of N., Ni, ve, vi 
and a;. The results of calculations are given 
in Table VIII. 


Table VIII. Viscosity coefficient of the gaseous core. 
= Atomic | 
Reh Bye 10 20 | 30 
i nace i ™| | 
3000 | 2.16 x10-?+3.40x10-5 1.85 x 10-2 +4.35 x 10-5 | — _ 
5000 | 2.80-x 10-8+ 3.22 x 10-5 2.40 x10-2+ 4.30 x 10-5 — —_ 
8000 3.14 x 10-3+3.00 x 10-5 3.04 x10-3+4.05 x 10-5 2.65 x10-?+ 5.08 x 10-5 
11000 | 4.15 x 10-3+-2.72.x 10-5 3.55 x 10-3 4+ 3.92 x 10-5 3.12 x 10-34. 4.90 x 10-5 
15000 | 4.88x10-%+2.40x 10-5 4.16 x10-3+3.71 x 10-5 3.65 x 10-2 +4.80 x 10-5 
20000 | 5.63 x 10-3+1.98 x 10-5 4.81 x 10-?+3.46 x 10-° 4.20 x 10-344.60 x 10-5 


The first and second terms in Table VIII 
are due to ions and electrons, respectively. 

In §1, we tried to see whether the earth’s 
core is gas or liquid by comparing Grpps’s 
energy for each case. But we could not 
' get a decisive conclusion. In this section we 
have evaluated the viscosity coefficient, the 
quality which is sensitive to the difference 
in state, and we have obtained the value of 
10--~10-! poise for the liquid state and 10-° 
poise for the gaseous state. But we shall 
not be able to distinguish the viscosity coef- 
ficients by means of observations in near 
future. Therefore, whether the earth’s core 
is gas or liquid cannot be known from vis- 
cosity either. 


§3. Attenuation of Seismic Waves 

The attenuation coefficient for longitudinal 
Seismic waves per one wave length depends 
on the frequency and on the medium through 
which the waves are propagated. In case 


that seismic waves travel through the earth’s 
core, the attenuation coefficient can be given 
as follows: 


‘i Kk 
fen NI) 


where «: Static compressibility, 
xs: Adiabatic compressibility. 


: : : 8x 
Since «=x; in our case, a; is equal to “3 KON. 
If we put the value of y evaluated in the 


preceding séction into the above equation, 
we find 
a2 x10-“4~2 x 10-38. 

Therefore, the amplitude of longitudinal waves 
decreases by the factor e-10-™* _o-10-™ per one 
wave length when the seismic waves travel 
through the earth’s core. 

The velocity of transverse waves in the 
core is about 10-1cm./sec. from the following 
equation 
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na 
0 

Since the attenuation coefficient of trans- 
verse waves per one wave length is 47, the 
amplitude decreases by the factor e~** per 
one wave length. 

From the above results we can see that 
the longitudinal waves travelling through the 
core are scarcely attenuated. 


§4. Electric Conductivity 

There is no direct method of measuring 
the eleetric conductivity in the central part 
of the earth. The depths to which the cur- 
rents induced by the variation of geomagnetic 
field penetrate depend on its period of varia- 
tion. By analysing many types of variation 
of geomagnetic field, T. Rixirake (1950,.1951) 
found the electric conductivity of the part 
shallower than 1500 km. His result is that the 
electric conductivity of the portion shallower 
than 400 km. is 10-° ohm-!-cm-! while that 
deeper than 400 km..is 10-? ohm-! cm-?. 

Recently W. M..Eusasser and E. C. BuLtuarp 
have explained the geomagnetic field by assum- 
ing the liquid motion in the earth’s core. In 
their theories, the electric conductivity of the 
earth’s core was taken 4.5x10° ohm-! cm-!. 
This value is very different from that of 
Rikxirake. Let us evalute the electric con- 
ductivity in the earth’s core along the line 
of this paper. 

The electric conductivity is given by 

NE etl 
18 ye TDi cot 

from the statistical mechanics, where v, is 
the velocity of electrons at the top of the 
filled cells and is given by 


; ={ hh? (S Ne 2/3) 1/2 
i Ame \ 7 7) 3 


From these equations, we can calculate the 
electric conductivity. by inserting the values 
obtained in the preceding sections, if we 
assume 

Ma y A 
ta? N; 

The results of numerical calculations are 

as follows; 
Atomic Number 10 
2.5 x 10'~1.9 x 103 ohm-! cm.-} 


Atomic Number 20 
3.1 x 10?~2.8 x 10 ohm-! cm.-? 
Atomic Number 30 
4.0 x 10?~3.6 x 10? ohm-! cm.-1 
The thermal conductivity can be obtained 
easily from the law of WIEDEMANN and FRANZ, 


3z {fk 
h= 3 (5) Te. 


The results are as follows; 
Atomic Number 10 
0.0b5~0.25teala deg.—* cmigeisecs: 
Atomic Number 20 
0:06~0.37 cal. deg.—+ cm.-! sec.-? 
Atomic Number 30 
0.08~0.47 cal. deg.-t cm.-! sec.-? 
In the present paper, it has been attempted 
to investigate the physical state of the earth’s 
core from the stand point of quantum statis- 
tical mechanics. It is admitted that there are 
not a few points which need further consider- 
ation.. It is hoped that these will be revised 
by the writer or others in the future. 


Acknowledgment 
The writer wishes to express his hearty 
thanks to Prof. K. Sassa, Prof. E. Nisaimura 
and Prof. S. Mryamoro for their instructions, 


encouragements and many valuable discus- 


sions. Many thanks are also due to Mr. J. 
NisHimuraA and Mr. M. Nisar for their as- 
sistance in this work. 


References 
EYRNG, H., etc.: 
1941 
Hill Co., New York, 480-485. 
MipOi dae 
1952 ‘‘GRUNEISEN’s Parameter and the Depth 
of Isostatie Compensation’’. Journ. Phys. 
Harth. teelo: 


RAMSEY, W.H.: 


1949 ‘On the Nature of the Earth’s Core. 
M.N.R.A.S., Geophys. Suppl., 5, 409. 
RICH, OKs 
1940 «Electronic Structure and Chemical Bind- 


ing’’. 
376. 
RIKITAKE, T.: 
1950 ‘Electromagnetic Induction within the 
Earth and its Relation to the Electrical 
State of the Earth’s Interior’. Bull. 
Earthq. Res. Inst., 28, 45; 219. 


“Electromagnetic Induction within the 
Earth and its Relation to the Electrical 
State of the Harth’s Interior’. Bull. 
Earthq. Res. Inst., 29, 263. 


McGraw-Hill Co., New York, 370- 


“The Theory of Rate Processes’?. McGraw- 


JOURNAL OF PHYSICS OF THE EARTH, VOL. 1, No. 2, 1952. 5 


Characteristic Tilt of the Ground that Preceded the Occurrence of 
the Strong Earthquake of March 7, 1952. 
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. 


Abstract 


At one of the tiltmeter stations operated by us, an unusually large tilt of the 
ground was observed during 3 months preceding the occurrence of the strong earth- 
quake of March 7, 1952, of which the epicentre happened to be not very far from the 


station. 


After careful investigations of the reliability of the observations, it has been con- 
cluded that the ground tilt was really in an intimate connection with the earthquake 


occurrence. 


It is emphatically stated that the observations of the ground tilt are one of 


the promising methods useful for predicting the occurrences of destructive earthquakes. 


§1. Conspicuous upheaval and subsidence 
of the ground have often been observed in 
epicentral areas of destructive earthquakes 
in connection with their occurrences. Also 
the surveys of precise levelling and _ trian- 
gulation have shown that these deformations 
of the ground take place not only at the time 
of earthquake occurrences, but also during 
several years or months preceding them. The 
deformations of the ground before and after 
the occurrences of earthquakes are in the 
opposite sense in some cases, while they are 
in the same sense in others. In the former 
-cases, the rebound motion of the ground may 
be presumed. Continuous observations of 
the deformations of the ground have been 
believed to be one of the promising means 
for predicting the occurrences of destructive 
earthquakes. 

For this purpose, dense nets and frequent 
occupations of the surveys of levelling and 
triangulation in the areas in question are 
undoubtedly highly desirable. But these will 
both mean prohibitably expensive and _ labor- 


ious works. 
As one of the supplementations of these 


works, continuous recordings of tilts of the 
ground by highly sensitive tiltmeters at as 
many stations as practicable may be recom- 
mended. In fact, since about 1930, such tilt- 


meter observations have been continued by 
us at several stations in Japan. In recent 
years, observations with extensometer, gravity 
meter, magnetometer and others were added 
at the stations, with a hope to get further 
informations, if possible, about the deforma- 
tions and state changes of the earth’s crust 
which may take place concurrently with 
seismic activities. 

Concerning with tiltmeter observations, 
characteristic tilts of the ground have been 
observed shortly before. the occurrences of 
several destructive earthquakes. These were 
studied by K. Sassa and E. Nisaimura (1951) 
and were believed to be directly connected 
with the seismic phenomena. The tilts are 
very small in amount, but peculiar in form 
and, in most cases, are likely to become not- 
able several hours before the occurrences of 
earthquakes. Since they have their own 
characteristic mode, the tilts are easily dis- 
tinguishable from those due to disturbances 
that have no connection with the occurrences 
of earthquakes. It does not seem right, how- 
ever, at the present stage of our knowledge 
of these matters, to insist that every destruc- 
tive earthquake is necessarily preceded by 
such a forerunning phenomenon and that, if 
it occurs, it can be observed at a station 
situated in any direction with respect to the 
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epicentre. Long and patient observations are 
needed before this important problem can be 
solved. . 

§2. Considerable difficulties are inherent in 
the observations of tilts of the ground, owing 
to the circumstances to be stated below. So 
far as our experiences are concerned, the tilts 
of the ground which were observed several 
hours before the occurrences of destructive 
earthquakes were about 0’’.1 in amount at 
observation stations lying within a distance 
of 100km. from the respective epicentres. 
Such minute tilts can be weil observed only 
when tilts due to meteorological and other 
disturbances are eliminated to a_ sufficient 
degree. This is only possible at a deep seated 
observation room which is suitably selected 
and proved stable. Moreover, for obtaining 
a clear and unmistakable record of that 
characteristic tilt, the observations must be 
made within a distance of 100 km. at least, 
or preferably at a much less distance, from 
the epicentre. The direction of the station 
as seen from the epicentre is also an import- 
ant factor to be considered, since the amount 
of tilt may differ considerably in different 
directions, or even may be zero in a certain 
direction, owing to the mechanism of occur- 
rence of particular earthquakes. If we are 
to establish tiltmeter stations in Japan so that 
any earthquake epicentre in her neighbour- 
hood shall lie within 50km. from any one of 
the stations, about 100 stations are necessary, 
while the distance of 100 km. is allowed, the 
number will be reduced to 40. In the last 
quarter of this century, about 30 destructive 
earthquakes have taken place here and there 
in Japan. At present, only about 10 tiltmeter 
stations are being operated by us. It is regret- 
fully admitted that under these circumstances, 
the study of the forerunning phenomena of 
earthquakes is very slow in progress and has 
much to depend upon mere lucky chance. 

It has been stated that characteristic tilts 
of the ground have been observed several 
hours before the occurrences of destructive 
earthquakes. Although these phenomena are 
of prime importance for studying the possible 
change in the state of the earth’s crust in 


respective hypocentral areas right before the 
occurrence of the earthquakes and also for 
studying the mode of release of seismic 
energies, they are too small in amount to be 
taken as an unmistakable sign of the comming 
catastrophe and if they occur, they do so too 
shortly before the earthquakes to be of any 
practical value for making the findings public 
in order to prevent disasters. 

For the public prediction of earthquakes to 
be of any reliable certainty and of any prac- 
tical effectiveness, such phenomena as are in 
an unmistakably direct connection with the 
occurrences of earthquakes and as will apear 
several days, several months, or even several 


years before the occurrences must be sought 
and continuously watched. 


§3. In this paper, the present writer is go- 
ing to study secular tilts of the ground be- 
fore the occurrences of destructive earth- 
quakes, especially from the standpoint of 
prediction. That gradual and long persistent 
tilt of the ground occurs before the occurences 
of destructive earthquakes has been establish- 
ed by means of levelling surveys. This is a 
very useful phenomenon for our purpose, but 
if we are to observe it by tiltmeters, a num- 
ber of difficult problems arise. At under- 
ground observation rooms, if they are built 
shallower than 20m., it is our experience 
that the tilts due to daily and annual varia- 
tions of temperature, heavy rainfall and to 
other unknown origins are so conspicuously 
recorded down there that those due to earth- 
quakes, if they exist, can hardly be recog- 
nized. If the room is several hundred meters 
below the earth’s surface, no such disturb- 
ances are present, but instead another diffi- 
culty arises. A large continuous tilt into a 
certain direction is usually observed. This is 
likely to be a result of the enormous earth- 
pressure working around the room, for the 
station is usually situated in an adit of a mine. 
Such tilts depend on geological structure, 
topographies, as well as other conditions of 
the mine. Delicate tiltmetric observations in 
a coal mine adit are almost impossible owing 
to this reason. In our case, a stable place in 
a metal mine has carefully been selected, 
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where the above mentioned secular deforma- 
tions are negligible. a 

On the other hand, in long observations of 
secular change ‘in tilt, cautions must be paid 
against apparent tilts caused by instrumental 
defects. It is ordinarily difficult to decide 
whether the observed secular tilts are real 
ones which are our object to watch or they 
are merely apparent ones due to secular 
change on the part of the tiltmeter used. 
Concurrent observations with two or more 
tiltmeters of the same type at the same place 
are highly desirable. In addition, observa- 
tions with tiltmeters of different types, at 
the same place, or at different adits are desir- 
able. These comparative studies are now 
being made and it is hoped that the results 
will shortly be published. 

§4. What will be described below is the 
results of our tiltmetric observations made at 
a deep room in a mine near which a very 
strong earthquake occurred on March 7, 1952. 

The epicentre of this earthquake has been 
located at a point 20km. offshore towards 
south-west from Noto Peninsula, Ishikawa 
Prefecture. The exact position of the epi- 
centre was at 136°.2 E and 36°.5 N. The 
energy of the earthquake has been estimated 
to be about 107% ergs. An earthquake with 
a nearly same magnitude took place at 136°.3E 
and 36°.1N, on June 28, 1948, not very far 
from the one in question. Since the 1948 
~ earthquake took place on land, it wrought 
tremendous damages to the city of Fukui and 
its. neighbourhood near its epicentre. 

Since the 1948 earthquake, tiltmetric obser- 
vations has been continuously being made at 
the depth of about 300m. in an adit of the 
Ogoya copper mine. The place is at 136°.5E 
and 36°.3N and is therefore 40 km. distant 
from the epicentre of the 1952 earthquake. 
(Fig. 1). The surrounding mountains are 
mainly Tertiary tuff in geological formation. 

The tiltmeter observations were started 
soon after the 1948 earthquake, together with 
those with extensometers and magnetometers, 
primarily for the purpose of studying post- 
seismic deformations of the earth’s crust. 

§5. Before going to study the tilts observed 
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quake and the locations of the tiltmetric 
observation stations. 
x Epicentre A Stations 
in connection with the 1952 earthquake, it 
will be appropriate to describe our instrument 
and the method of observations in some 
detail. 

From the day of beginning of observations 
until the end of 1950, the tiltmeters were set 
up directly on the floor of natural rock ex- 
posed in the adit. From September 1951, 
they were set up on a concrete foundation 
prepared at the same place. There was an 
interruption of observations of 10 months 
which intervened between the two periods. 
The effect of change in position of observa- 
tion will be discussed in detail later. 

On March 4, 1952, 3 days before the occur- 
rence of the March 7 earthquake, which we 
are going to study, another earthquake 
occurred at 144°.0E and 42°.0N, south of 
Hokkaido. This earthquake was a big one, 
having the energy 4 times as large as the 
March 7 earthquake. In spite of the great 
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distance of 900km. between the epicentre 
and Ogoya, which is our tiltmeter station, 
the tiltmeters were badly shaken by the 
earthquake waves, owing probably to reson- 
ance of the tiltmeter pendulums to the wave 
motions. 4 days were needed to repair and 
adjust the instruments, so that the tiltmetric 
observations were interrupted from March 4 
to 8. Unfortunately enough, the March 7 
earthquake happened to take place just in 
this interruption. It is more than regrettable 
that a chance was lost for completely study- 
ing the forerunning phenomena of a destruc- 
tive earthquake, although this was an un- 
avoidable circumstance. 

At our station at Ogoya, tiltmeters, extenso- 
meters and magnetometers are working. The 
tiltmeters are essentially Zollner-suspension 
type horizontal pendulums, almost wholly 
made of super-invar. The supporting frame 
of the pendulum is about 20cm. high. The 
mass of the pendulum which is 15gr. in 
weight is suspended by super-invar wire, 30 
in diameter. The equivalent pendulum length 
is 6cm. As may be seen in Fig. 2, the de- 
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Fig. 2. Tiltmeter of horizontal pendulum-type. 
The cylinder on the left is the cover. 

flection of the pendulum caused by tilt of the 
ground is photographically recorded by the 
lamp and mirror method at a distance of 2m. 
A drum around which a photographic paper 
is wound makes a revolution in one week. 
The velocity of running is 1.7mm. per hour. 

The sensitivity of the tiltmeters are deter- 
mined either by giving them a small tilt of a- 
known amount or by measuring the periods 
of their free oscillations. 

The sensitivities of our tiltmeters were kept 


a 1 +—$— 
B ieee 
CA ie A Spahr 
7 
Yo 
Fh TALE Le 
oat B 
r T i as = 4 40 
| A 
| | a! ie | Earthquake 60. 
| 
in Ss anes ska a —1__§— ed 
1948 1949 1950 1951 1952 | 
Fig. 38. Secular tilts of the ground at Ogoya. Arrows indicate the time of occurrence 


of the Earthquake of March 7, 1952. 


Tilt of Ground that Preceded the Earthquake of March 7, 1952. 


at 0.013’’"/1mm. deflection from July 1948- 
October 1950 and at 0.078’/1mm. deflection 
from September 1951. The periods of free 
oscillation are 30 sec. and 12 sec. respectively. 

§6. Two tiltmeters are used at our Ogoya 
station, the one to record the N-S component 
of the tilt and the other to record the E-W 
component. These will be called (A) and (B) 
respectively. In Fig. 3, the records obtained 
by (A) and (B) during the last 4 years are 
shown. If there were interruptions of obser- 
vations, the curves before and after them are 
connected reasonably by broken lines. 

Referring to the (A)-component, the ground 
tilted a little toward N soon after the obser- 
vations were started in July, 1948, and after 
2 months, the direction of tilt turned toward 
S and remained so for the following 2 years. 
Since July, 1950, a gradual tilt toward N be- 
came again apparent. Since September, 1951, 
when the tiltmeters were reset after an inter- 
ruption of 10 months as stated above, the 
tiltmeter had been showing a large tilt toward 
N, until the occurrence of the March 7 earth- 
quake. It is noteworthy that the tilt became 
smaller after the earthquake. 

As to the (B)-component, the tilt was very 
small during the first 2 years. After the 
tiltmeter was reset in September 1951, small 
fluctuations in tilt were observed for some time. 
But, since the beginning of December, 1951, 
an exceedingly large and characteristic tilt to- 
~ ward W became notable and continued up to 
the time of the earthquake occurrence. After 
the event, the tilt recovered its direction to- 
ward E and judging from its present rate, 
it is likely that the eastward tilt will before 
long cover the previous westward tilt. 

The sensitivities of the tiltmeters were fre- 
quently calibrated by measuring their free 
oscillation periods. The sensitivicy of a hori- 
zontal pendulum type tiltmeter depends on 
the small angle 7 between the directions the 
axis of its rotation and of the plumb line. 
The angle of rotation @ caused by a tilt ¢ 
in the direction perpendicular to the plane 
containing the two directions is given by @ 
=g/i. The free oscillation period T of pen- 
dulum is given by 
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T=2n/ Tg, 

where / is the length of the equivalent simple 
pendulum. Since the angle 7 is very small, 
it is easily affected by the tilt of the ground 
in the plane containing the rotation axis and 
the plumb line. In our case, the angle 7 is 
about 54’ or 340” corresponding to the period 
of free oscillation of 30 sec. or 12 sec. tre- 
spectively. The ground tilt of 10’ or 60” 
will cause the change in the period by 3.0 sec. 
or 1.5 sec. Frequent precise measurements 
of the free oscillation periods of a tiltmeter 
are useful not only for the purpose of its” 
calibration. A large tilt observed by a tilt- 
meter in one direction can be proved to be 
a real one by observing the change in free 
oscillation period of another which is set at 
right angle with it. That the large tilts ob- 


served by our tiltmeters were real is beyond 
doubt. 


These large tilts cannot be of local nature; 
for instance, they are due neither to defor- 
mations of the concrete block on which the 
tiltmeters are placed, nor to those of the 
adit due to earth-pressure. Tilts due tu these 
local deformations are likely to be gradual 
and one-directional. As the curves in Fig. 3 
indicate, the observed tilts are sometimes 
gradual and other times rapid and not infre- 
quently change their direction. Especially 
before and after earthquake occurrences, the 
observed tilts are very characteristic in ap- 
pearance. Strictly speaking, the above state- 
ments are not fully ascertained unless con- 
current observations with several tiltmeters 
are made at different adit of the Ogoya mine. 
But the characteristic tilts observed several 
months before and after the earthquake, especi- 
aily those in EW direction observed by the 
(B)-component are believed to be intimately 
related to the occurrence of the earthquake. 
As was stated before, the tiltmeters were 
reset in September 1951. After that, the (A)- 
component tiltmeter continued to record a 
large and monotonous tilt towards N, and the 
tilt amounted to nearly 65’’ just before the 
outbreak of the earthquake. After the event, 
the northward tilt decreased in its rate as 
shown in Fig. 3. But since February 10, 


80 Kennosuke HosoyAMA. 


1952, the tilt increased a Jittle in rate. As to 
the (B)-component, on the other hand, the 
rate of tilt was small during the first 3 months 
after resetting, but since the beginning of 
December, 1951, a rapid tilt toward west be- 
came notable. This westward tilt amounted 
to nearly 55’ in the succeeding 3 months up 
to the occurrence of the earthquake. Corres- 
ponding to the notable change in rate of the 
northward tilt on Feb. 10, a similar change 
is also found also on the (B)-component record. 
Especially after Feb. 25, the westward tilt 
reached 15’’ in 10 days. It is remarkable 
that this rapid tilt ceased soon after the 
earthquake and its rate at present is as it 
was before. Combining two components of 
the tilt, it was found that the ground at 
Ogoya began to tilt toward NW. 3 months 
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Fig. 4. Tiltgrams obtained at Ogoya during the epochs of (a) October 12-18, 1951 
February 18-25, 1952 and (ce) May 14-21, 1952. ; 


before the occurrence of the March 7 earth- 
quake. It is notable that the epicentre of 
this earthquake was also 40 km. NW from 
Ogoya. A more rapid NW tilt began 10 days 
before the earthquake, but its rate became 
suddenly small soon after it. At present, the 
tilt is as it was before. The problem of this 
recovery of tilt is a very interesting and im- 
portant one, which calls for its solution long 
and patient observations. The tilt records 
obtained during the said three periods are 
shown in Fig. 4. Corresponding to the 
change in rate of tilt observed at Ogoya 10 
days before the earthquake occurrence, one 
which is similar in shape but much smaller 
in amount was also observed at some other 
tiltmeter stations. These are Kamigamo, 
Ikuno and Kishu, which are respectively, 160, 


: (b) 
Lines indicated by A,B, M and 


O are the ground tilts to NS- and EW-directions, change of magnetic declination 
> 


and datum line respectively. 
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190 and 280km. distant from the epicentre 
as shown in Fig. 1. The tilt records obtained 
at these stations showed characteristic tilts 
which were about 0.1, 0’’.2 during 10 days 
before the earthquake occurrence but which 
ceased soon after it. This small amount of 
tilt is roughly one hundredth of that observed 
at Ogoya, and this ratio is likely to be reason- 
able considering the difference of distances 
of the stations from the epicentre. The 
secular tilts observed at these stations from 
Feb. 7 to March 13, 1952, are shown in Fig. 5. 

$7. Summary 

An example of characteristic tilts of the 
ground before and after the occurrence of 
an earthquake has been observed at a tilt- 
meter station which is not very far from the 
epicentre. This was really a uniquely for- 
tunate chance which has never been met with 
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Fig. 5. Ground tilts observed at Kamigamo, 
Ikuno and Kishu, during the period of Feb- 
ruary 7-March 13, 1952. Arrows indicate the 
time of occurrence of the earthquake of March 
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in our 20 year tiltmeter observations at several 
stations. The observations have been valuable 
in two points: first the station was only 40 
km. distant from the epicentre and secondly 
the observations have been made long enough 
so that characteristic tilts both before and 
after the earthquake were successfully ob- 
tained. In the beginning part of this paper, 
the necessity and difficulty of long tiltmeter 
observations were emphatically stated in con- 
nection with the problem of earthquake pre- 
diction. Then follows a description of de- 
tails of observations at Ogoya. After care- 
ful investigations of the record obtained at 
Ogoya from various points, the conclusion 
was reached that the large and characteristic 
tilts of the ground observed at Ogoya 3 
months before the March 7 earthquake, more 
particularly those observed 10 days before it 
are intimately related to the occurrence of 
the earthquake. It is concluded that tiltmetric 
observations of this kind must be continued 
in order to solve the important problem of 
earthquake prediction. 


Acknowledgements 

The writer wishes to express his cordial 
thanks to Professor E. Nisammura for the 
guidance throughout this study and also to 
Professor K. Sassa for his advice in inter- 
preting and analysing the observed records. 
It is a pleasure to the writer to record with 
thanks the good understanding and kind help 
of the staff members of the Ogoya Mine 
Cooperation, which the writer received during 
the observations and without which the ob- 
servations could not have been possible. The 
cost for the observations was covered by the 
Grant in Aid for Scientific Research of the 
Ministry of Education. 


Reference 


Sassa, K., and NISHIMURA, E.: 


1951 “‘On Phenomena Forerunning LHarth- 
quakes”. Trans. Amer. Geophys. Union, 
32, 1-6. 


‘ 


< a! 


Lae iret) al : 


AF) ies irae 


‘Seae 4 dé: eats 
eee}: ey aed d 
Magee! 3 


in ' Rie it P og N Wt ae Lelie 


We 


fen PF 


Parte 
» VIR ae 


@' 


JOURNAL OF PHYSICS OF THE EARTH, Vou. 1, No. 2, 1952 83 
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al la Efikoj de la Nankaido Tertremo, la 21-an de Decembro, 
1946, en Koti Prefektejo, Sikoku, Japanujo.* 
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Abstract 


Geographical Distribution of Earthquake Damages with Special Reference 
to the Effects of the Nankaido Earthquake, Dec. 21, 1946, in 
Kochi Prefecture, Shikoku, Japan. 


Concerning the relations between geomorphological ground conditions and earth- 
quake damages throughout the Prefecture of Kochi, Shikoku, Japan, due to the Great 
Nankaido Earthquake of December 21st, 1946, we have noticed several interesting facts: 

i) On the sandy ground of beach ridges, damages were found only in those parts, 
where the basement bed rock lies deep, or where some kind of soft mud layers lies 
beneath sand, or where the ground consists of very soft sand, which has been deposited 
in very recent years or accumulated by wind blow. On the contrary, there was scarcely 
any damage, where the bed rock lies near the surface, or where hard materials such 
as pebbles or stones lie under the sandy ground, or where the surface sand are con- 
solidated. Practically, of course, some of the above mentioned conditions coexist com- 
plicatedly combined in the same place. 

ii) There were severe damages on swampy low lands behind the beach ridge. 

iii) Even on alluvial plains, we found patches of rather slight damages among other- 
wise severe damages. Under those patches of slight damages stony layers were usually 
found as is the case on beach ridges. 

We have examined the hardness of the surface ground by a simple hand-boring 
stick of about 1m. length. 

Existence of the stony layers under the ground was conjectured by the distribution 
of Cl’ contents in underground water. Where there are stony layers under the ground, 
underground water flows constantly and the Cl contents must be small. On the con- 
trary, underground water cannot flow easily in clay, mud or fine sand layers, and the 
Cl contents must become great there. A typical example of the above two cases was 
encountered in the town of Nakamura, where the geographical distribution of earth- 
quake damages clearly agrees with the Cl’ contents distributions. (Fig. 10) 


* International language ESPERANTO. 
** Tokio Toritu Daigaku signifas Cefurba Universitato de Tokio at Metropolitan University of 
Tokyo. La junaj aittoroj apartenis gis 1949 al la Geografia Instituto, Naturscienca Fakultato, 


Universitato de Tokio. 


Setumi MiyAmura, Yatuka Saizyo kaj Sohei KAIZUKA. 


84 


‘ofnueder ‘nyoytg ‘ofeyyeferg Oy 
ep fozyoqun fey fosepia (nig ue owosrjA0y, opreyueN SEES Bl op q fouropsoy ap aes Rasa! ep peo I4SIp ea iy Pra) Aes 


(IVNVY 4 O¥-d avy) 


Wodd lunzisy} 


Beco, © O0l=d >0s@ OWNHNS 2 
PENGTN Va OS > d>02 @/LNIVMY HOW IS 9 
2 = aon = Og > d>g © VLV9O ¢ 

TEE G>d>/l @® VYNWVMVN pb 
arte 1>d=0 o JNON ¢ 

Od 0 VMVOVYIN 2 

WO 

Atoka d rogvwoa IM 

FQ CVLNII0¥d 


WHO? 


aie 


> 


OF 
dees 
SHA. ts fh 


wv 
OPE! ecel 


OLYVWONVT VSYWIONI 


Geografia Distribuo de Tertremaj Domagoj. 85 


De Ja granda tertremo de 21-a de decembro, 
1946, kiu enerale estas nomata Nankaido- 
au Nankai-Tertremo, la plej severaj domagoj 
estis suferataj en la suda parto de Sikoku 
Insulo, nome en la Prefektejo Koti. La 
epicentro de ¢i tiu tertremo estis proksimume 
Ge 135°E, 33°N lati provizora publikigo de la 
Centra Meteorologia Observatorio. Sed la 
domagoj de la tertremo montris tre kompli- 
kitan distribuon en Giuj partoj de la sud- 
okcidenta Japanujo lati la diversaj naturoj de 
grundoj kaj aliaj kondi¢oj, sendepende de 
epicentra distanco. Ekzemple, la _ urbeto 
Nakamura kaj gia Girkatiajo en la kantono 
Hata, Koti-Prefektejo, suferis la plej grandajn 
domagojn de ¢i tiu tertremo, kvankam la loko 
situacias Cirkatii 200 km malproksime de la 
epicentro. 

Pri la domagoj de Ja tertremo en Koti- 
Prefektejo, esploroj de diversaj vidpunktoj 
kaj diversaj metodoj estis jam faritaj de 
kelkaj kolegoj en nia instituto, nome, de d-roj 
K. Kanar, H. Kawasumi, T. Minaxami, S. 
SakuMA, Y. Sato kaj aliaj. La nunaj atitoroj 
esploris la temon precipe pri interrilatoj de 
tertremaj domagoj kaj naturoj de grundoj de 
la vidpunktoj de geomorfologio, hidrologio 
kaj geologio. Nu sekvos sube la generala 
aspekto de la rezulto de niaj esploroj pri Giu 
geomorfologia distrikto aparte, kiu estas 
montrata en Fig. 1. 


I. Distrikto de Muroto-duoninsulo. 


La tuta distrikto apartenas administre al la 
kantono Aki kaj gia tereno estas plejparte 
montaro kun mallargaj valoj. Sekve vilagoj 
kaj urbetoj tie situacias generale sur la 
marbordo kiel strat-vilaSoj, escepte:de kelkaj, 
kiuj trovigas sur la roka grundo de la mont- 
flanko de valoj en interna lando. Lokoj de 
urbetoj kaj vilaSoj estas klasifikataj lati la 
naturoj de grundo kiel sube: 

a) roka grundo, b) riverateraso, c) mar- 
borda teraso, d) malalta perono de marerozio, 
e) marborda digo, f) cirkonferenco de aluvio 
g) interna parto de aluvio, h) aluvio Girkat 
river-enfluejo. 

Tamen plejmulto da vilagoj kaj urbetoj 
situacias sur la grundoj a), d) kaj e). Surla 


grundoj a) kaj d) oni trovis neniajn domagojn 
de la tertremo, kaj sur e) trovigis konsider- 
indaj domagoj en kelkaj partoj. Sur la 
grundo f) estis malgrandaj domagoj kaj sur 
g) kaj h) trovigis iom grandaj domagoj. 

Tamen Cie la grado de la tertremaj domagoj 
estis tre komplikitaj lati la aliaj diversaj 
kondicoj de lokoj. 

Kiel supre skribite, riveroj en la distrikto 
havas tre mallargajn valojn kaj formas nur 
malgrandajn aluviajn ebenejojn ée siaj buSsoj 
al maro. Tial la dikecoj de tiuj aluvioj estas 
malgrandaj kaj, en konsekvenco, oni ne 
trovas tie tiel grandajn tertremajn domagojn, 
kiel oni generale povas trovi sur tiel-nomata 
aluvia ebenejo. Precipe, ju pli suden oni iras 
des pli granda oni trovas la lastatempan 
levigon de grundo, la terasojn altaj, kaj 
levigintan peronon de marerozio kreskanta, kaj 
konsekvence, kvankam la aluviaj ebenejoj ¢e 
river-enfluejoj estas sufice largaj, la aluvio 
estas tre maldika kaj sur riverfundo malko- 
vrigas roka grundo, ekzemple en la valo de 
Huttu, en la urbeto Muroto. Tertremaj 
domagoj estas sendube tre malgrandegaj tie. 
En la norda parto, la marbordo transigas al 
malleviga tipo de leviga tipo, kaj pli norde de 
cirkat’ Aki la riverbuSaj ebenejoj estas iomete 
pli larga kaj ankaw iliaj aluviaj tavoloj estas 
pli dikaj. Sekve trovigas iom da domagoj en 
vilagoj kaj urbetoj sur tieaj aluvioj. Latlonge 
de la marbordo kreskas sablejo Gie, krom tie, 
kie la marbordo estas borderita de leviginta 
perono de mar-erozio. Kreskado de la sablejo 
estas tre rimarkinda en la norda parto de la 
distrikto, kaj latilonge de la marbordo de 
riverbusa aluvio-ebenejo kreskas bekforma 
att digforma sablajo, malantati kiuj estas 
malaltaj kaj malsekaj terenoj. Marborda 
sabldigo formas du ati tri paralelajn vicojn 
respondante al siaj nombroj de fojoj de 
nekontinuaj leviéoj en la pasintaj jaro). 

Sub la sablejo, kiu borderas la terasojn 
proksimigantajn al la maro aii montojn, estas 
roka bazo malprofunda kaj la tertremaj 
domagoj tie estas ne tre severaj. Kontrate 
la sablejo, kiu borderas la riverbusgajn aluvio- 
ebenejojn, estas bekforma sablajo Ce golf- 
enirejo att sablodigo en bordmalproksimeco 
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antati la aluvioj estis enterigitaj kaj sekve la 
dikeco de sablejo estas granda kaj en kelkaj 
partoj sur la tiaj sablejoj trovigis sufiée 
grandaj domagoj de la tertremo. Nome, kiel 
oni renkontas en la urbetoj Aki, Kiragawa 
kaj None, severe suferintaj partoj en la 
urbetoj lokigéas generale aii sur sablotavolo 
pli nova ol aliaj partoj ati sur sablejoj, sub 
kiuj roka bazo estas malprofunda, kaj precipe 
tiaj partoj, kiaj konsistas el blovita tavolo de 
sablejo sur la malseka tereno malantati la 
marborda digo ati sablajo, suferis tre grandajn 
domagojn de la tertremo. Tio Sajnas kom- 
preneble ne katizita de tio, ke en tiuj partoj 
speciale amasigas malbonaj domoj. (Komparu 
kazojn de Aki kaj None.) 

Sube ni priskribos la rezultojn de la esploroj 
en kelkaj lokoj de la distrikto. 

(1) Urbeto Aki. La urbeto kugas latilonge 
sur du vicoj de sablaj digoj, el kiuj norda 
estas iomete pli alta, kaj la topografiaj 


konturlinioj en la urbeto montrigas en Fig. 
2. (Ankati vidu Fig. 1.) 


qOm 
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Fig. 2. Domagoj de.la Nankaido 'Tertremo en 


la urbeto Aki, Koti Prefektejo. 


De nordo al oriento fluas unu _ rivereto, 
kontratt kiun malaltigas la tereno kaj ankati 
en la nord-okcidento la tereno estas malalta, 
tusante al la malantatia terciara montaro. 
Tiu sablejo estis verSajne bord-malproksima 
sablajo. Gia meza parto estas la ple] mainova, 
dum la ambati ekstremoj estas pli novaj, kaj 
el la du vicoj la norda estas sendube pli 
malnova. Tamen, la suprajo de la norda 
sabl-digo Sajnas formita de blovado. 

Lati la diroj de putfosistoj, subtera tavolo 
en norda parto estas speco de maréa tero kaj 


sub la sablejo de la orienta parto kusas 
riverstonetoj, dum alie estas nur sufi¢e dika 
tavolo de sablo. 

Ni ekzamenis Ja distribuon de Cl’-enhavo 
de la subtera akvo en la urbeto. Generale la 
distribuo de Cl’-enhavo en subtera akvo 
proksima de marbordo estas paralela al la 
marborda linio, kiel montrate de la esploro de 
d-ro Sinki@i Yosurmura, en Simoda, Sizuoka 
Prefektejo, sed ¢i tie tiela rilato tute ne 
ekzistas. Krom marborda regiono Cl’-distribuo 
estas Generale superregata de influoj de sterko 
en la kampo kaj higiena stato en domaro, kaj 
tio estis videbla ankaii Gi tie. Plie, Gar ni 
scias, ke la Cl’-enhavo estas multe influata de 
movado de subterakvo, ni povas konjekti per 
tio, ke la akvoportanta tavolo sub la tereno 
estas maréa ati Stonetplena terajo. Generale 
en maréa tavolo la subterakvo staras senmova 
kaj kontratie i fluadas en Stonetplena tavolo. 
La distribuo de Cl’-enhavo en la urbeto Aki 
estas montrata en Fig. 3 kun nombroj de 
gram-valoloj por litro. Kvankam la influo de 
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Fig. 3. Cl’-enhavo en subterakvo kaj rezulto 
de borbastona ekzameno en la urbeto 
Aki, Koti Prefektejo. 

(x mola sablajo A malmola sablajo) 
kampo en norda ¢ikatiajo kaj de malpura 
distrikto konsistanta el Svarmo de domaéoj en 
nord-okcidenta parto estas ambaii iomete 
klara, la plej granda tendenco, kiu duigas la 
urbeton en partojn okcidentan kaj orientan, 
estas tio, ke en la orienta parto la subtera 
tavolo konsistas el riverStonetoj kaj sekve tie 
fluadas pura subterakvo kun malmulta Cl- 
enhavo, kaj en la okcidenta parto la subtera 
terajo estas Slimo, argilo, ati sablo havanta 
malpuran kaj Cl’-riéan subterakvon. Kom- 
preneble en la orienta parto trovigas ankati 
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Tabelo ye Domagoj en la urbeto Aki, Koti Prefektejo. 
Logdomoj Tenejoj, Staloj, keto pe 

Nombroj de | Pro- Nomb: ae ail aoe 
No. Suburbeto Detruitaj Tutaj _centaj| Detr ore = cane 
“Tute | Duone | 1") | Doma- “Tute | Duone  T%t@j | Doma: 

a met ee No Ne B. 80) P No | iNn Nv g0j P 
1 1-ku it SHH 182 2.9 0.7 
2 2-ku | 3 | 106 1.4 2 1 2.4 
3 3-ku celine: hie ieee BT a heey a | 0.0 
4 4-ku ee Sere FT os | eel 2 2.1 
5 5-ku i) Bt) | Soa a soe de | 1 2 2.4 
6 6-ku ane 0 2.6 0.0 
7 | 7-ku 66 0.0 2 1.5 
8 8-ku 10). 196 4.4 6 6 71 
9 9-ku | | 109 0.0 | | 0.0 
10 | 11-ku | 1 | 128 0.8 | 2 1.6 
11 12-ku | lates) 19% 1.2 1 1 1.2 
12 13-ku 3 BO) : 
| 104. 6.8 4 2 4.8 

13 | 14-ku 5 6 104 77 "33 10 26.9 
14 | 15-ku 2 10 TORO. 1K) 0.0 
15 16-ku 10 13 84 | 19.7 CW use| | 22.0 
16 17-ku 18 28 88 | 27.3 4 | 1-7 5.1 
17 18-ku 8 47 | 96 | 82.8 3 | 3.1 
18 | 19-ku bbvep uel’ 45 127 cae Ole ed lata ae 16 15.2 
19 20-ku 13 104 106. (ase 7) 10 8 deed 
20 21-ku 6 20 163 tA Hert ie 2 } Ste2 
Pak | 22-ku Dey AT 102 bet mh | 6 Sal: 
22 23-ku il | 45 Zila | e259) | 0.0 
23 24-ku | 140 0.0 | | 0.0 
24 25-ku 2 She MAL. Qed hice a ah 8.9 
25 Anauti B} Bis Nh Ope! | | 0.0 
| Sumo 83 406 3160 | Beata 80s) e SA | 4.0 


N.B. P={(Ng+Np12)/N} x 100%. 
ke N estas sama kun tio de loSdomoj. 


MD <3 O% 300M 


Fig. 4. Procentaj domagoj P de la 
Nankaido Tertremo en la urbeto 
Kiragawa. 


domaro de malpura stato! 

Por ekzameni solidecon de plej supra terajo, 
ni uzis borbastonon (kies longeco kaj diametro 
estas ¢. 15m kaj 2cm _ respektive) kun 
&talpinto ensovante gin en la_ teron. La 
rezultoj de borbaston-ekzameno en la urbeto 
montrigas en Fig. 3 kun la signoj x kaj 4, 
kiuj montras molajn kaj malmolajn lokojn 


P de tenejoj k.t.p. estas kalkulata, kun supozo 
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Fig. 5. Cl’enhavo en la subtera akro en la 
vilago de Kiragawa, Koti Prefektejo. 


ENHAVO a a 


DE CL’ 


respektive. 

Domaj domagoj de la lasta granda Nankai- 
tertremo en éiuj kvartaloj (suburbetoj) de la 
urbeto Aki estas montrata en Tabelo I kaj 
Fig. 2. La domagoj estas generale severaj 
sur nova kaj moja tereno, kaj en la orienta 
parto, sub kie kuSas_ riverstonetoj, trovigas 
nenia granda domago. En la orienta ekstremo 
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Tabelo II. Domagoj en la urbeto Kiragawa. 
| Logdomoj Tenejoj, Staloj, k.t.p. 
S een ~ _Nombroj de | | Pro- Nombroj de Pro: 
Ong) POUDUEVEtO | spetruiray . | centaj Detruitaj . | centaj 
Tute | Duone oe Doma-| Tute | Duone ae Doma- 
| Ng | Nn, 80) Dee Ng Ni, gO) Ie 
pan) , > see ae as | { 
1 | Higasi-mati 1(1)| 14(15) 54 | 14.9 8(11) | 14( 24) 27.8 
20 | Nici-mati 4(4)| 85(89)/ 66 | 34.8 | 4( 4)) 18( 23) 19.7 
3. Naka-mati 1(1)| 10711) 59 | 10.2 | 9 9)| 14( 28) 2771 
4 | Kami-mati-nisi | 12(16) | 45 13.4 4( 4)) 11¢ 15) 21.1 
5 | Kami-mati- 8( 8) 54 7.4 | 9( 12) 8.3 
 higasi | | 
6  Miya-no-uti 1 8) ATi obT5 1S SC SuIeeGe 6.4 
7 | Kita-mura | 10 0.0 2 2) 10.0 
8) | Linzi | 1¢ 1} 36 1.4 | 2 2) 2.8 
9  Naka-no-kawa | | | 63 0.0 6( 6) 4.8 
10 Turi-no-kuti | 26 0.0 | ete a) 1.9 
11 | Nada 369 0.0 |  2¢ 2)) 18( 19)| 3.0 
12 | Aliaj 414 0.0 — 0.0 
| Sumo 6(6)| 87(98)| 1243 4.0 | 29(83)| 97(130) | 6.2 


T oa aaa 


de la urbeto trovigas iom da domagoj videblaj 
kompare al en la centra parto. Tamen en la 
3-a kvartalo aperas jam influo de riverstonetoj 
kaj la domagoj estas iomete malpli severa. 
Domoj en la okcidenta parto estas sufi¢e mal- 
fortikaj, sed tamen tio verSajne estas katizo 
de domagoj tie, éar la kazo de la vilaSo None 
montras la rilaton tute inversan kiel ni 
priskribos pli poste. 

(2) Urbeto Kiragawa. Kiel montrate en 
Fig. 4, la vilaSo kuSas sur la montflanko 
proksimiganta al la maro kaj sur la sablajoj, 
kiuj de la mont-bazo etendigas al la direktoj 


orienta kaj okcikenta respektive. La rivero, 


kiu enfluas en la maron en la orienta ekstremo 
de Ja vilago, havas neniun aluvion. Tie roka 
bazo estas pli malprofunda kaj riverstonetaro 
pli granda ol en okcidenta limo. 

Tertremaj domagoj estis severaj sur okci- 
denta sablajo kaj la parto sur la montflanko 
neniom suferis. (Vidu Tablon II kaj Fig 4.) 
Cl’-enhavo distribuigas kiel en Fig. 5 kaj 
la influo de subkuSanta fluo de rivero estas 
granda en okcidenta parto. Sekve la rilato 
de Cl’-enhavo en la subterakvo al la tertremaj 
domagoj estas negativa, kaj la domagoj estas 
Gefe regataj nur de la profundeco de roka bazo 
subtera. Sajnas probable, ke la sablajo de 


Tabelo II]. Domagoj en la vilago None. 
Logdomoj Tenejoj, Staloj, k.t.p. 
é ae Nombroj de | Pro- Nombroj de Pro- 
No Subvilago Detruitaj | Tutas centaj|  Detruitaj © . | eentaj 
Tute | Duone aa ES Doma- ‘ute | Duone Tutaj | Doma- 
- ea ees | te) PN, leant NO poi 
1 | Ura-1-ku lpanccdis ise Al 993 | 30.7 7 32 23.7 
2-| Ura-2-ku Leto 68 g2. | 5s.) i “18 27 25.8 
3 | Ura-3-ku V7 |, Bt 98 | 58.7 27 53 51.0 
4 | Higasi-mati el hcl 315 | 25.4 | . 31 96 24.6 
5 | Ike-no-aima 8 23 . 49 | 39.8 17 38 25.5 
6 | Nakamura 35 45 83 | 72.0 55 83 93.7 
7 | Nakazima 3 11 44 | 19.3 2 8 10.9 
8 | Kuzukago 3 6 16 87.5 2 9 31.4 
9 | Osino | 1 8 idehag7 Ais GBI 2 1.2 
10 | Utida | 6 ST dees k 2 tee 
11 | Oto | 4 31 6.5 3 5.8 
12 | Kawaguti | 3 36 4.2 2 1.4 
13 | Masunase | | 3 82 4.7 1 0.9 
14 | Betuyaku ans. 32° | 3.2 1 0.8 
15 | Natome-gawa 4 | 46 | 108 | 25-0 8 21 7.5 
Sumo 99 | A821 080 ihe Stall Nee 15 alae re 20.2 
| | 
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la okcidenta parto de la vilago konsistis el 
blovita sablo de pli marproksima tereno, kaj 
konsekvence la severaj domagoj tie farigas 
klaraj. 

(3) Vilago None. Distribuo en Giuj sub- 
vilagoj estas montrata en Tabelo III. Sub- 
vilago Nakamura (N-ro 6), kiu kuSas meze de 
la aluvio, suferis la plej severajn domagojn 
kaj Ura-2, Ura-3, sur la marborda digo sekvis 


SIS/KUI 


Fig. 6. Indeksa Landkarto de la subvilagoj 
en la Vilago None, Koti Prefektejo. 


< 30% 
fm <20% 


Fig. 7. Domagoj de la Nankaido Tertremo en 
la vilago None, Koti Prefektejo. 


gin en domagoj. Estas malgrandaj domagoj 
ankatt en Ura-1, Higasimati, Ike-no-aima, 


Nakazima, kaj Kuzukago; la wunuaj du 
situacias sur la marborda digo de sabla terajo 
kai la lastaj tri estas Girkati kaj en la aluvia 
ebenejo. (Vidu Figurojn 6 kaj 7.) 

Kiel ¢i supre skribite, la domagoj en Naka- 
mura estas tre grandaj. Surfaca terajo tie 
konsistas el komplekso de argilo kaj stonetoj, 
kaj tio montras, ke gi estas riveramaso. Sed 
lai geomorfologia vidpunkto ni povas 
konjekti, ke subtere estas malgrandaj graj- 
netoj de argilo kaj sekve la tereno estis iam 
laguno malanatat: la marborda digo. 

Montflanko éirkat la aluvio estas tre kruta 
kaj tio montras la grandan dikecon de la 
aluvio. La lastaj konjektoj probable klarigas 
katizon de la severeco de la tertremaj domagoj 
tie. 

Sur la marbordo de la sud-okcidenta parto 
(Ura-2 kaj Ura 3) suferas severe de la 
tertremo, kaj en la nord-orienta parto (Ura-1 
kaj Higasimati) la domagoj estas malgrandaj. 
(Komparu la paragrafon pri Aki urbeto.) 

Cl’-enhavo de put-akvo en la vilago estas 
generale tre malmulta kaj influigas iomete de 
irigacita rizokampo en okcidento kaj svar- 
mado de domaéoj en nord-orienta parto. En 
nord-orienta parto ni povas supozi movigadon 
de subterakvo. (Vidu Fig. 8) 


CL'-ENHAVO 
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Fig. 8. Cl’-enhavo en subterakvo kaj rezulto 
de borbastona ekzameno en la vilago None, 


Koti Prefektejo. (x mola A malmola) 


Lati la borbastona ekzameno, surfaca terajo 
tie trovigis mola cirkonference de la marborda 
digo. Denaske la marborda digo kreskis de 
la nord-oriento kie] bekforma sablajo lati nia 
konjekto, tial la distribuo de domagoj klarigas 
de noveco de naskigo de la sablajo. Ankatt 
sub la nord-orienta parto estas roka bazo 
malprofunda. 
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II. Riverregiono de la Rivero Mononobe kaj 
la Girkatiajo de la Urbo Koti. 


La valo de la supra parto de la rivero 
Mononobe estas relative mallarga sur si kaj 
en diversaj lokoj trovigas riverterasoj kun 
maldika amaso. Tertremaj domagoj estas tie 
malgrandaj. 

La orienta parto de Koti-ebenejo estas 
okupata de ventumilforma tereno de malsupra 
parto de R. Mononobe kaj trovigas tie 
nemalmultaj urbetoj kaj vilaSoj. Tertremaj 
domagoj sendube estas malgrandaj, ¢ar 
generale la ventumilforma tereno estas stonet- 
plena kaj rezistas forte al tertremo. Kontrat 
tio, la okcidenta parto, ¢irkatiajo de la urbo 
Koti konsistas el aluviaj ebenejoj, kiuj 
malaltiéas inter dislokigo0j, kurantaj de oriento 
al okcidento, tereno tie estas malalta kaj 
malseka kaj oni trovas sufi¢e severajn 
tertremajn domagojn. 

Partoj de la urbo Koti, kiuj estas mal- 
fortaj al tertremo, estis plene forbruligita de 
usona aeratako kaj tial procentaj domagoj de 
la tertremo estas tre malgrandaj malgrat 
supozo. Sekve, por tiu ¢@i distrikto estas 
necesa singardo. 

Marborda distrikto de Koti-ebenejo kon- 
sistas el marborda digo kaj malalta kaj 
malseka tereno malantati gi. Ekzemple en la 
urbetoj Kisimoto, Miwa kaj aliaj trovigas 
severaj domagoj sur ¢i tiuj malaltaj terenoj 
malantat la digoj. 


Ill. 


La supra parto de ‘‘antecedanta’’ rivero 
Yosino ne havas largan valon, kaj sekve 
trovigas malmultaj vilagoj en la regiono. 
Tertremaj domagoj ne ekzistas tie, katize de 
solida grund-naturo. 


Riverregiono de la R. Yosino. 


IV. Riverregiono de Ja R. Niyodo. 

En la supra parto trovigas malgrandaj 
domagoj de la tertremo, sed la malsupra 
parto suferas iom severe kaj estas pli dangere 
tie ol la riverregiono de la R. Mononobe. 


V. Regiono de la 
Simando. 

Meza parto de R. Simando meandras 

enmetate konsekvence de lasta ter-levigo kaj 


supra parto de R. 


havas riverterason, sur kiu estas dislokitaj 
urbetoj kaj vilagoj. Ne trovigas tre severaj 
tertremaj domagoj, krom en urbeto Kubo- 
kawa, kie procentaj domagoj superas 1%. 
Tamen la lasta urbeto situacias sur iom larga 
inunda aluvio, kiu naskigis pro tio, ke nur tiu 
parto mallevigis en plej nova geologia epoko 
kontrati pli malsupra parto, kiu_ levigis 
samtempe. 


VI. Riverregiono de la malsupra parto de 
R. Simando, Nakasuzi Terfosajo k.a. 

En Ja centro de la regiono situacias la 
urbeto Nakamura, kiu suferis la plej severajn 
domagojn de Nankaido-tertremo kiel jam 
skribite. La malsupra parto dela R. Simando 
nomigas ankati R. Watari, kaj gi fluas tie 
ankati kun iom kruta klinigo. Kontratie $iaj 
enfluantoj Nakasuzi kaj Usiro havas tre 
malkrutajn riverujojn. Precipe la lito de R. 
Nakasuzi estas pli malalta ol tiu de R. Watari 
(Gefa fluo) Ge la enfluejo, kaj gis 8.5km supra 
punkto gia fluosuprajo restas horizonta. ‘Tial 
la supra riverregiono de R. Watari estas tre 
vasta kaj okaze de inundo ekfluadas grandega 
kvanto da akvo, parto de akvo enfluas tiam 
kontratiflue en la fluejon de R. Nakasuzi kaj 
R. Usiro ofte okazigante superakvegon sur 
la Girkatiajo. 


Fig. 9. Geomorfologia Historio éirkaij 
Nakamura, kantono de Hata, Koti 
Prefektejo. 


La terfosajo de Nakasuzi (germane Naka- 
suzi-graben) naskigis post la kreta periodo, 
sed la diluvia tavolo tie suferas nenian agiti- 
gon, kaj tio montras, ke la dislokigoj ne 
movigas post la kvarta periodo. Tial oni ne 
povas rilatigi la dislokigojn kun nuna ter- 
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tremo, malgrati asertoj de kelkaj atitoroj. 

Ebenejo en tiu distrikto estas tre dika 
aluvio, kies amaso konsistas el §limo el mal- 
grandegaj eroj, kaj e¢ tre proksime de monto 
roka bazo estas tiel profunda, kiel 50m 
borado ne povas atingi gin. La geomorfologia 
historio de la regiono povas doni bonan klarigon 
al tio. Nome, kiel montrate en Fig. 9, unue 
naskigis largaj peronoj katizitaj de mar- ai 
river-erozio, due ili levigis kaj formis terasojn, 
ekzemple la malaltan terason norde de urbeto 
Simoda, trie ili plie levigis intermitante kaj 
formas kelkajn Stupojn de terasoj (ekzemple 
marterasoj kiuj subportas subvilagon Hirano, 
urbeton Simoda, ati subvilagon Iya, vilagon 
Taguti kaj riverterason, sur kiuj restas 
subvilago Taguro, vilago Gudo.) Levigoj de 
€i tiu periodo estas sufi¢e grandaj, la erozio 
datiras longe, kaj formas profundajn valojn, 
kun tre krutaj flankoj. Kvare, tiam la 
tereno malievigis kaj maro invadis en la 
valojn, kie akumligis dika aluvia tavolo @is 
la marsuprajo. Post tiam R. Simando multigis 
sian transportan forton, katizis levigon en sia 
mez-flua regiono, kaj portis stonetojn, kiuj 
kovris maldike la suprajon de la aluvio. (Vidu 
priskribon pri la urbeto Nakamura.) 

La baza roko en tiu distrikto estas alternaj 
tavoloj de sablostonoj kaj skistargilo. Mal- 
graui ke la kialo estas kompreneble nekonataj 
partoj de skistargilo detruigas tre facile pro 
efloresko. Tio alportas tertremajn domagojn 
pro terfalajo en la regiono de Nakasuzi kaj 
ankatit havus iom da influo al generala 
severeco de domagoj en la tuta distrikto. 

Proksime al montaro trovigas tie ankat 
konsiderindaj domagoj de la tertremo malgratt 
generalaj ekzemploj. Krom la supre citita 
naturo de roka bazo, tio devenas de artefarita 
amasigo, kiu estas populara en tiu regiono aul 
por akiri largan konstruan terenon ¢e la 
kruta montflanko ati por eviti superakvon de 
oftaj inundoj. 

Kruteco de valformo estas klara per sisma 
ter-esploro de d-ro K. Kanai kaj ekzameno 
de putoj proksime de montflanko. 

(1) Urbeto Nakamura. Pri la distribuo de 
tertremaj domagoj en la urbeto d-ro Kanai 
jam aliloke skribis.» Nu sekvos Gi-sube kelkaj 


CL-ENHAVO 
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Fig. 10. Cl’-enhayo en subterakvo sub la 
urbeto Nakamura, Koti Prefektejo. 


rimarkoj pri interrilatoj de domagoj kaj la 
naturoj de grundo. 

En la orienta duono de Ja urbeto suferas 
pli severe de la tertremo kaj tio bone 
respondas al la distribuo de riéa Cl’-enhavo 
en subtera akvo kiel montrate en Fig. 10. 

La subtera akvo en la okcidenta parto 
enhayas pli malmulte da klor-jono ol en la 
orienta parto, Gar la unua fluadas tra la 
stonetplena akvoportanta tavolo amasigita de 
R. Simando kaj la lasta restas, kontrate, 
en la Slima subtera tavolo amasigita de R. 
Usiro, kaj tiu malsameco inter du _ partoj 
donas sendube supre priskribitan distribuon 
de tertremaj domagoj. Kiel aliaj atitoroj 
skribis, Giuj domoj en la urbetoj konstruigis sur 
amasigataj terenoj por bari la oftajn super- 
akvegojn. 


VIII. Marborda distrikto de sud-okcidenta 


parto de la kantono Hata. 


La distrikto enhavas marbordon de ¢irkata- 
jo de Suzaki, preter promontoro Asizuri 
gis la urbeto Sukumo. Kvankam en la 
norda parto ne trovigas marborda teraso, oni 
trovas gin de ¢irkatt Saga kaj, malsame de 
kazo de Muroto duoninsulo, vilagoj restas sur 
&i. Kompreneble estas ankati vilagoj sur 
marborda digo kaj malsupro de monto pro- 
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ksimigas al la maro. Sube sekvos rezultoj de 
esploroj en la vilagoj Ogata, Simokawaguti 
kaj la urbeto Sukumo. 

(1) Vilago Ogata. Nur en la subvilago 
Irino ni faris niajn esplorojn. La tereno estas 
malalta sutplenigita sablejo inter marbordo 
malproksima sablajo kaj teraso kaj suferas 
treege de la tertremo. Kiel montrate en 
Fig. 11, domagoj estis severaj en Hayasaki (1), 
Siba (6), Mangyo (5) kaj Mati(3). Aliflanke la 
domagoj estis iom malseveraj en Honmura (4) 


PROCENTA 
. DOMAGO 


BESS < 50% 
KZA < 30% 
(ml < /0% 


Bigs J. 
Irino, vilago Ogata, Koti Prefektejo. 


Domagoj de la Nankaido Tertremo en 


CL~ENHAVO 


MPC’ CO« 50MG/L 
°c MM <IOOMG/L 
on EQ >/OOMG/L 


Fig. 12. Cl’-enhavo en subterakvo en la 
subvilago Irino, vilago Ogata, Koti 
Prefektejo. 


proksima al malsupro de teraso kaj en 
Hamanomiya (2), kiu situacias sur la sablejo de 


marborda digo. Cl’-enhavo en la subtera akvo 
distribuigas tie kiel en Fig. 12, kie klare 
montrigas la parto kun staranta subterakvo. 
Granda Cl’-enhavo en Mangyo devenas ankat 
de malpura stato de domaro. 

(2) VilaZo Simokawaguti. Tertremaj do- 
magoj estas montrataj en Fig. 18, Severaj 
domagoj distribuigas en la malalta tereno 
maiantaii sablajo kaj en la malantatia flanko 
de la sablajo. Kontratie sur la sablajo mem 
precipe en la proksima parto de roka bazo 
sin trovas neniu domago. 


PAZ, 
Waitin, 


CAE egies 


Fig. 13. Domagoj de la Nankaido Tertremo 
en la vilago Simokawaguti, Koti 
Prefektejo. (x tute A duone) 


Ciuj supraj respondas al la distribuo de Cl’- 
enhavo en la subterakvo, montrata en Fig. 14. 
Oni ne povas rimarki subteran fluon influatan 
de proksima rivero. La valoroj de Cl’-enhavo 
mem estas tre grandaj, sed la kialo estas 
neklara. 


ENHAVO:: DE-CL’-(MG/L)-. \\. 
ZZA<50 wa<200 wi 
150 


ES > 200 


Fig. 14: Cl’-enhavo en subtera akvo en la 
vilago de Simokawaguti, Koti 
Prefektejo. 
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BE FORBRULO 


eae 


Fig. 15. Domagoj de la Nankaido Tertremo 
en la urbeto Sukumo, Koti Prefektejo. 
(x tuta detruigo A duona detruigo) 


(3) Urbeto Sukumo. Sukumo estas inter- 
landa urbeto kaj ne marborda. Kvankam 
la subterakva Cl’-enhavo klare dividas la 
urbeton en du partojn tute same kiel en la 
urbeto Nakamura, tamen estas ne tiel klara 
rilato al la tertremaj domagoj en tiu ¢i okazo. 
(Vidu Fig. 15) La domagoj estas regataj Gefe 
de la dikeco de aluvio, kiu farigas iom post 
iom maldika proksimigante al la monto. 


Fine la atitoroj deziras esprimi siajn korajn 
dankojn al la logantoj kaj la attoritatularo 


EZZAZ < 2OMG/L 


LE, V 
CL™-ENHAVO =} > SOMG/L 


Fig. 16. Cl’-enhavo en subterakvo en la 
urbeto Sukumo, Koti Prefektejo. 


éie en Koti Prefektejo, kiuj bonvole helpis 
iliajn esplorojn. Ankati dankon al d-ro K. 
Kanat de la Instituto -por Esploroj de 
Tertremo por lia bonkora rilatigo al la estraro 
de la prefektejo. S-rinoj K. Noumoto kaj S. 
Konpo helpis la atitorojn en preparado de 
l’figuroj. 
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The First and Second Vertical Derivatives of Gravity. 


By 


Chuji Tsusor1 
Geophysical Institute, Faculty of Science, Tokyo University, Tokyo, 


and 


Motohiko Ka'ro 
Teikoku Oil Company, Tokyo. 


Abstract 


A method has been shown to calculate the distribution of 0g/0z from that of g. 
The second derivative method has been critically reviewed. 


§1. A paper by H.M. Evsen (1936) was 
likely to be the first to direct our interest 
anew to the problems of the vertical deriva- 
tives of gravity. More recently, the second 
derivatives have become the subject of many 
studies. Particularly THomas A. Exxins (1951) 
has given a detailed description of the second 
derivative method of gravity interpretation 
and by several examples, has shown how 
this method is effective for detecting small 
irregularities in gravity anomalies and is 
therefore useful for figuring out minute 
underground mass distribution which are 
liable to be overlooked by the ordinary me- 
thod. There can be no objection to ELxrins’ 
statements expressed in his paper. It is 
interesting to note that he wrote ‘‘its (the 
method) use is justified only on data of high 
accuracy’’. 


§2. The second derivative method consists 
essentially of calculating oe and Be from a 
given distribution of g and of summing up 
the two. It is evident from the LapLacr’s 


equation, that this sum will give with 


the algebraic sign reversed. In the LapLacr’s 


equation 

AOE ON ag OMe OM 0 

Oat ak Oy* = Oz,’ 
where U is the gravity potential, the last 
term on the left side is nothing but 0g/0z, 
while the first two have no direct connection 


with g and can therefore not be calculated 
from g at once. Here lies the reason why 
Evxkins differentiated each term of the Lap- 
LACE’S equation with respect to z and obtained 
the relation 


This is a very important relation in that a 
possibility is thus opened to calculate a ver- 
tical derivative (second) of g from its hori- 
zontal derivative (second). But at the same 
time, this was made possible at the cost of 
accuracy, because the second derivatives are 
too much sensitive to small errors in gravity 
measurements as Exxins admitted. In our 
opinion, what was aimed at by the second 
derivative method has been not so much to 
obtain the very quantity 0°g/0z" as to get one 
which has a definite physical meaning and is 
reasonably sensitive to small irregularities in 
gravity distributions and which can be cal- 
culated directly from them. 


§3. Let us consider a thin sheet of mass 
at the depth d from the earth’s surface and 
the mass along the sheet be expressed by 

M (ay) = X1d1 Cn sin 2X sin 2Y- 
The gravitational potential due to this sheet 
is given by 


—V 2 +n? (z+) 


U (xyz) 
cos mx cos ny e - 


9 + Cinn 
=2nh > ye pin sin 
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where z is taken vertically upward with the 
origin at the earth’s surface. k® is the uni-- 
versal constant of gravitation. The potential 
is so taken that its derivative in a certain 
direction will at once give the force in that 


g(xy)=2nk* Seu Cin ate NX oe 


direction. From the expression of the poten- 
tial, it is evident that the horizontal distribu- 
tions of gravity and its vertical derivatives 
are given by 


-y, eS oe 
ny e m2+n2 a 


—V m+n? a 


28 ary) = — Ark DS Vm mn Sn My € ; 


24 sy) =2nk* SD (m2 n®) Cnn 8 8 
es 


COS yyy 4° COS ny e —Vmi+n? a 


If the observed horizontal distribution of gravity is expressed by 


g(xy) = Sy Bmn ses mx es ny ; 


fen 29 Smite Buin 2% mx 28 my , 
ond eos S (m+n?) Bina 8 mx ny . 


Since these expressions do not involve d, they 
can be extended to the case of any three 
dimensional underground mass distribution. 
There are two points to be noticed. The 
one is that the first vertical derivative of 
gravity can also be calculated from the hori- 
zontal distribution of gravity, if arithmetical 
Jabours in doing so are not spared. If the 
horizontal distribution of gravity is expressed 
by a double Fourtrr series and if each coef- 
ficient is multiplied by the corresponding 
Vm?+n? and again synthesized, the results 
will at once give the horizontal distribution 
of the first vertical derivative of gravity. 
The computational labour is not so great as 
it at first appears to be. With well designed 
stencils for harmonic analyses and a calculat- 
ing machine, the computations can be carried 
out in a relatively short time. In a forth- 
coming paper, one of the authors is going 
to describe his results of calculation of 0g/0z 
from an actual distribution of gravity. 

The second point to be noticed is that by the 
expressions given above, we can see to what 
degree 0g/0z and 6’9/0z" are sensitive against 
local horizontal gravity variations. If these 
variations are regarded to consist of various 
components with different wave-lengths, the 


former is sensitive in the first power of the 
reciprocal of the wave-length, while the 
latter is so in its second power. This is also: 
the very reason that 0?g/0z* is sensitive to 
observation errors. 

After all, if the object of the second deri- 
vative method is really to see the distribution. 
of 0°9/0z*, the method of ELxins is to be 
adopted. But it is to be remembered that 
09/0z can also be calculated from gravity 
distribution. If the object of the second 
derivative method is merely to get a quantity 
which is reasonably sensitive to small irre- 
gularities in gravity distribution, the first 
power quantities seem to be preferable. An 


example of these is / (GE) +5 if and this 


quantitity can be calculated easily from a 
given gravity distribution. 
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Anomalies of the Vertical Gradient of Gravity 


Associated with Anomalies of Gravity. 


By 


Chuji Tsusotr. 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


The gravity anomaly is necessarily associated with the anomaly of the vertical 


gradient of gravity. 
invest gated. 


The relation between the two kinds of anomaly has been 
The anomaly of the vertical gradient is sometimes too large to be 


neglected in the free air reduction of gravity. 


§1. In a previous paper, (Tsusor and 
Karo, 1952), the writer has shown with M. 
Karo that if the distribution of gravity anoma- 
ly along the earth’s surface is expressed by 


Agiz)=> Bu mea , 
this is necessarily associated with the anomaly 
of the vertical gradient of gravity which can 
be written by 


429 = —S mBp we me , 
Oz : 


G5), Gm Se 71a).. 


This expressions gives the anomaly per one 
radian. In order to obtain its value per unit 
length, the right hand side of the above equa- 
tion must be divided by Z/2z, where Z is the 
length which is taken as 2x in the Fourier 


expression. Thus we have 


2xzm 
A ae ip Bun Sema : 


Since Z/y is nothing but the wave-length 1», 
of the 72th harmonics, we have 


0 
4 7 _> Grice 012 ? 


2r 
=— > Bngime , 


Or Gn =— — 

m 
It is therefore seen that the amplitude, Gyn, of 
AOg/0z is 27/4 times that of the correspond- 
ing gravity anomaiy, B,. The first vertical 


gradient is usually taken as 0.3086 mgal/m, 
that is 3086 E.U. But cautions must be taken 
for the fact that, in many cases, its anomaly 
is by no means negligible as the numerical 
values in Table I show. 

Table I. Amplitudes of the Anomaly of 0g/0z 


in Eotvés Unit for Various Gravity Anomaly 
Amplitudes and Wave-Lengths. 


SS Baal paises 1 Oseal | 100meal ; 
Am ~~. cet pete ee 
km | 
1 See eS i4! |, or628 6283 
10 | 6 63 | 628 
100 | al | 6 63 


The gravity anomaly of 10 mgal having the 
wave-length of several km. is not uncommon, 
especially where detailed gravimetric survey 
is to be made. At such places, it must always 
be kept in mind that the normal value of 
0.3086 may easily vary by more than its 10 
per cent. 

It is true, however, that 40g/0z itself 
changes quickly in the vertical direction, so 
that in order to make an accurate free air 
reduction, it is sometimes necessary to go 
beyond the first derivative. 

§2. The gravity potential due to a plane 
mass which is situated at the depth d is given 
by 

Ue2=> Am Ug Hip chee as 
m 


where the positive direction of z is taken 
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vertically upward. The potential is taken in 
such a way that if it is differentiated in a 
certain direction, it will at once give the com- 
ponent of force in that direction. From this 
expression for the potential, it follows that 
the gravity at the height z is given by 


Ag(x2z)= —L= > Am° mxe-GFD ‘ 


But since 
Aga, 0)= Xi Amo mae-™4 , 
=>) Bn(0)28 mz , 
we can write 
Ag(az)= > Bri@ysma= > BnO)gnmne-™? . 


Each component of 4g(x, 0) and that of 4g(x, z) 
are connected through the relation 


Bn(Z)= BnO)e-™ 4 
or BO) => By(ze™ yi . 
If the first order only is taken, then, 


B»(O)= Byr(z)(1+ mz) 

= Bm(z)+m2zBm(zZ) 

=Bn(z)+mzBn(0) 

=Bnrlz)—zG nl(0) , 
where we have put B,(0)=B,(z) for the 
second term. So long as this is a good 
approximation, the neglect of the second and 
higher order derivatives are justified. But 
this is not always the case. Generally we 
have to compare (1+mz) and e”*. Besides, 
Biz) and B»(0) are not exactly equal. Since 
z is equal to 2xh/Z, where h is the corres- 
ponding actual height, we are to compare 
(1+2xzmh/L) and exp(2zmh/L), or the same 
thing, (1+2zh/An) and exp(2zh/Am). These 
will depend on the ratio h/2{, and we get the 


values given in Table II. The values in 


Table II. Comparison of exp(2ch/d,) and 

(1+2rh/4m) for Various Values of h/A,.- 
han exp(2rh/Am) | 1+2nh/Am 

0.00001 1.0000628 | 1. 0000628 i 
0.0001 1.000628 1.000628 
0.001 1.00629 | 1.00629 
0.01 1.0648 1.0628 
0.1 1.874 1.628 


Table II show that if the height is more than 
one per cent of the wave-length, the linear 
correction is already insufficient. 


§3. In order to make the circumstances — 
clear, let us take a numerical example. 
Assume an area with the horizontal extent 
100 km. and the gravity values at a constant 
height of 1km. be expressed by 


ara + 10 cos eee. 


g(ah)=const.+ 10 cos 100 100.’ 


where « is measured in km. and g in milligal. 
The effect of topography is considered to have 
been already corrected. If the uniform correc- 
tion of 0.3086 1000 is applied, we shall get 


5-272x 
100 


go—const.+10 cos ae 


+0.3086 x 1000 . 


If we take into account the anomalies of 0g/0z, 
we shall get 


A) +10 cos 


2r2x 5-2nax 
1 Tas, 
100 +10 cos 


+0.3086 x 1000 


B) g =const.+10 cos 


+00 x 10 x cos —— + 


100 
5-272 
il SS [Sci : 
x 10 x cos 100 |x 000 


100 


Table III. Values of {g)—const.—(0.3086 x 1000)} 
calculated according to the Methods A), B) 
and C), in mgal. 


ckm | Aye fl By} C) 
0 20 27 25 
5 10 11 10 

10 = 2 rah dae ks 

15 6 qT 6 

20 138 19) =| 17 

25 0 0 0 

30 —13 —19 —17 

35 — 6 - 4% = % 

40 2 7 5 

45 —10 —11 —10 

50 — 20 —27 —25 

55 ~10 ad: i et0 

60 2 7 5 

65 — 6 — 7 6 

70 -13 —19 Seley 

75 0 ie 0 

80 13 19 17 

85 6 7 6 

90 2 — 7 — 5 

95 10 11 10 

100 20. | 27 25 


If we take the exponential expression, we 
shall get 
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C) go=const. +0.3086 x 1000 
+] 10 cos. 775 exp( oe ) 


100 100 
5-272 5-27 
10.x : BELLE 
Ft a 9 exp( 100 if 


The calculated values {g,—const.—(0.3086 
x 1000)} for A), B), and C) are compared in 
Table III. 

As the values given in Table III show, the 
calculated free air reductions considerably 
differ according to the method of reduction. 
The difference amounts even to 7 mgal in 
. our present example. 

In actual case, the gravity stations are not 
situated at the same altitude, so that the above 
discussions cannot at once be applicable. The 
gravity anomaly along the earth’s surface 
cannot be known until the anomaly of the 
vertical gradient is known and the anomaly 
of the vertical gradient cannot be known until 
the gravity anomaly is known. There can 
therefore be no simple method by which to 
make a rigorous free air reduction. The only 
way left to us is likely to resort to the 
successive approximation. First apply the 
normal reduction (0.3086xh) and get the 
approximate distribution of gravity anomalies. 
From that distribution, caiculate 4dg/dz and 
making use of it, obtain the corrections to 
be made to the normal reduction of (0.3086 x h). 
Repeat this process until no further correc- 
tions will be needed. This is evidently a 
laborious work to do, but at least the first 
corrections are necessary in most cases. 

In the above discussions, the one dimen- 
sional case has been dealt with. The two 
dimensional case can be investigated quite 
similarly. If we use a double Fourter series 
instead of a single one and use 7/742+ 72 in- 
stead of 22, the above discussions will readily 
be applied. 

§ 4. In the case of a spherical earth, things 
are a little more complicated. Let us assume 
that the mass distribution on a_ spherical 
surface at the depth d is given by 


Mi HO Y n(ue) - 


The gravity potential outside this surface is 


é —d)m+l 
Tie ye (a oo AmY m(ug?) ? 


yme+l 
and that inside it is 


oe WB a 
(a—d)™ 
Since the potential is continuous at 7=(a—d), 
A,, and A’, must be equal. We can write 


U'= 


Alin Y n(u@) . 


therefore 
(a—d)™+t P 
U=>; ymel 7 Am Y m( uP) , 
and 
. ym 
= Ya—aym An mln?) 4 


Fig. 1. 


For r>(a—d), we have 


: Hii ) — m+1 

Ag(r)= — Ol == vis “ = a AmYm(uP). 
Or ym? 
On the other hand, putting r=a, we get 
= m+1 
4g(a)=>; EINE = a AmY m(uG) 
qm 
=>Bn Y n( uP) - 


We can therefore write 


M+2 


Ag 7) = SB -Y m( ue?) , 


m+2 


and 
2) m 
a ey = —> tages Z -Y m( 2@) : 
OP) r= a 


We see thus the amplitude of the mth 
harmonics of the anomaly of the vertical 
gradient of gravity is (7+2)/a times that of 
the corresponding harmonics of the anomaly 


100 


of gravity. 

The gravity values at the poles are 983 and 
that at the equator is 978, so that the varia- 
tion is 2.5 in amplitude and m=2 in order. 
Thus the difference of the vertical gradients 
of gravity at the poles and at the equator is 


_4x5 
6.4 x 108 


This is never a quantity negligible against the 
normal vaiue of 3086 10-°. 

If m=400, the corresponding wave-length 
is roughly 100 km, and if the amplitude of 
the gravity anomaly is 100 mgal, the difference 


=230%10-%. 
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of the largest and smallest vertical gradients 
is 
402 «0.1 
70 X10 . 
6.4 x 108 i 
This will cause a difference of 20 mgal in the 
free air reduction for an altitude of 3000 m. 
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